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Abstract 


A general characteristic of the atmosphere is its streakiness. It manifests itself in many pheno- 
mena such as the “cloudstreets’’ of polar outbreaks and tradewinds, the radar bands of hurricanes 
and certain precipitation patterns, the squall lines and the jet stream cirrus. 

The banded structure can be traced back to convection in heated air flows with a curved 
velocity profile of rather uniform direction. The typical vertical curvature was found to be 10-7 
(cm sec-1). A preliminary note deals with some theoretical considerations. 

Synoptic cloud analysis by weather satellites will greatly profit from the wide spread existence 


of cloud bands. 


I. Introduction 


The advent of the jet and rocket age pre- 
sents the meteorological observer with a reveal- 
ing bird’s eye view of the atmosphere. One of 
the most striking features of tropospheric 
weather as seen by a stratospheric observer is 
the streakiness of the cloud systems. Almost 
invariably a considerable fraction of cloud 
layers spreading under the observer is organized 
in bands running in the general direction of the 
flow (Figs. 3 to 5). With some practice one can 
estimate the wind direction at cloud level 
purely from looking at the cloud structure (a 
possibility to be exploited by meteorological 
satellites). The band structure of cloud systems 
is also apparent on many of the composite 
pictures taken from high flying balloons or 
rockets (Figs. 6 to 8). (BERGSTRALH 1947, Hu- 
BERT and BERG 1955). 

That the meteorological literature contains 
so little on this widespread phenomenon of the 
“cloudstreet” may be in part due to the worm’s 


eye view of the ground observer, in part it 
may be a consequence of the regrettable 
neglect, direct cloud observation in general 
and nephanalysis in particular have suffered in 
the last decade. A survey on the occurrence 
of cloudstreets indicates that they are frequent 
in all latitudes and, in some areas even represent 
the prevailing type of atmospheric convection 
clouds. Much information on the structure of 
cloudstreets has been supplied by glider pilots 
who are fond of this phenomenon and who 
introduced the term “wind thermals” to express 
their observation that a combination of thermal 
convection and strong winds is required to 
produce cloudstreets (KUETTNER 1949). How- 
ever, the priority of their discovery that the 
tedious circling technique in columnar thermals 
could be replaced by a comfortable tailwind 
flight beneath the cloudstreets must be given 
to the sea gulls (Woopcock 1942) who seem 
to have enjoyed this technique for millions of 
years. 

A more systematic study of the banded 


1 Present Affiliation: U.S. Army Ballistic Missile Agency 
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structure in the atmosphere is suggested by 
certain similarities with laboratory experi- 
ments. While a considerable literature exists 
on the artificial formation of “longitudinal 
rolls” or “strips” in flowing media heated 
from below, and on the analogy with certain 
cloud patterns in the altocumulus level (Mar 
1931, AVSEC 1939, BRUNT 1951) no satisfactory 
theoretical interpretation has been offered so 
far. As the equations governing this flow 
embrace the effects of viscosity, buoyancy and 
vorticity, the silence in this matter may well 
be a reflection of some discouraging efforts of 
workers in this field. 

A first step towards the solution of this 
problem is a more thorough understanding of 
the conditions under which bands form in the 
atmosphere. A general description of the phe- 
nomenon and the results of a special aerological 
investigation will be given in Part 1 of this 
study. Part II, to be published separately, will 
discuss the theory of organized convection in 
the atmosphere. 


2. General Description and Definition 
of »Cloudstreets» 


The most common type of cloudstreets may 
be defined as a system of long rows of cumulus 
clouds stretching in the general wind direction. 
The pictures, 1 to 8 give examples of convection 
bands as they occur in the atmosphere. From 
Fig. 1 and 2 it can be seen that the cloudstreets 
sometimes extend in parallel bands across 
the whole visible sky. Most of them show 
indications of individual cumuli lined up like 
pearls on a string. But the most spectacular 
cases of parallel cloudstreets resemble con- 
tinuous roll clouds (Fig. 3). Such clouds form 
during a similar diurnal period as normal 
cumulus clouds do. Although individual cumu- 
lus streets sometimes originate over and down- 
wind from a heat source, such as an industrial 
plant, a fire, a heated island, a warm lake or 
along a mountain range, the phenomenon 
discussed here is, in general, independent of 
orography. In fact, orography tends to disturb 
rather than to create cloudstreets and probably 
for that reason they are best developed over 
oceans. Interesting in this respect is the observa- 
tion that cloudstreets are the prevailing form 
of thermal convection during the Arctic sum- 
mer over Scandinavia (KUETTNER 1937). Here 
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the terrain consists mostly of tundra, lakes 
and hills, but none of the features of the 
topography is reflected in the cloud picture. 
Even over the Arctic pack-ice convection in 
the lower levels forms regular cloudstreets as 
far north as 80° during the almost uninterrupted 
solar heating of the summer months.? Simi- 
larly, observations in the tradewinds show that 
cloudstreets are a characteristic form of con- 
vective activity over the tropical oceans. 

Frequently the cumulus tops of cloudstreets 
have a definite ceiling, indicating the existence 
of a temperature inversion which restricts ver- 
tical motion to a well defined layer over the 
ground. For instance, the cloudstreets of North- 
ern Scandinavia form preferably in fresh out- 
breaks of Arctic air which is heated quickly 
from the ground during the long Arctic sum- 
mer days. An aerological study of this process 
(KUETTNER loc. cit.) shows how an adiabatic 
layer establishes itself through a depth of 2 to 3 
km, topped by a well defined inversion and 
maintained for as long as 15 hours daily. Most 
of these characteristics were already found in 
cloudstreet studies by the German Soaring 
Research Institute, DFS. (MALETZKE, 1953.) 

Some estimates of the length and spacing of 
cloud streets may be appropriate. Typical values 
are: so km in length and $ to 10 km spacing 
from axis to axis. However, cloudstreets of 
several hundred kilometers have been seen. 
With a characteristic layer depth of 1% to 3 
km the spacing of the cloudstreets lies fre- 
quently between 2 and 3 times the height of 
the convective layer. Occasionally isolated 
cloudstreets are observed. 

Patient observations reveal that cloudstreets 
have a cycle with a lifetime of the order of 
one hour. They appear to build up and dissipate 
on both ends, sometimes in such a fashion that 
the street as a whole drifts slower than the ele- 
ments in the street. This indicates that the street 
propagates upwind with respect to the mean 
wind and, in rare cases, even with respect 
to the ground. The dissipation is often seen to 
start through interference with other growing 
cloudstreets. 

From glider observations there is sufficient 
evidence that vigorous updrafts exist under cer- 
tain sections of the cloud base while down- 


? Observations by the writer on weather reconnais- 
sance from Alaska poleward during June 1957. 
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THE BAND STRUCTURE OF THE ATMOSPHERE 


Fig. 6. Cumulus in bands 
over New Mexico, photo- 
graphed from Aerobee rock- 
sm Oct, 1954. (Photo: 
Naval Research Lab.) 


drafts occupy essentially the space between the 
bands. Compared to the flow in columnar 
thermals the updrafts in cloudstreets are gener- 
ally estimated to be stronger under otherwise 
similar thermal conditions. 

More information about these details is 
available from the use of whole-sky or wide 
angle time lapse pictures, as are now being 
taken at the Harvard Blue Hill Observatory 
near Boston under contract with the Air Force 
Cambridge Research Center. Another source 
of information lies in pilot’s reports. A vivid 
glider flight account is given in the Appendix. 
Here Dr. Steinhoff describes his surprising ex- 
perience when he encountered and recorded 
this interesting phenomenon for the first time 
in flight history. In the light of subsequent ex- 
perience his observations on the distribution of 
drafts and turbulent zones are of special in- 
terest. 


Fig. 4. Early morning cumulus cloud streets, oriented 

in the direction of low level winds, as seen from 

appr. 10 km height over Tampa, Florida, 12 August 

1957. (Courtesy: V. Plank, Airforce Cambridge Re- 
search Center.) 


Fig. 5. Cloudstreets lined up in wind direction near 

Tampa, Florida, 30 August 1957, as seen from about 

5 km altitude. (Courtesy: V. Plank, Airforce Cambridge 
Research Center.) 


Fig. 7. Cumulus bands as seen from 150 km altitude, 
5 Oct 1954. (Aerobee rocket photo by Naval Re- 
search Lab.) 


Fig. 8. Composite picture of more than one million 

square kilometers of Northern Mexico and South- 

western United States. (Photos taken from V2 rocket 

by Naval Research Lab.) Date: 7 March 1947. Altitude: 

160 km. Cloudbands in lower right and center follow 
surface flow converging over Arizona. 
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3. Aerological and Synoptic Analysis 
of Extra-Tropical Cloud-Streets 


This study was stimulated by the following 
thoughts. The reports by glider pilots that 
only the fast drifting “wind thermals” trans- 
form into streets and the corresponding evi- 
dence from laboratory experiments suggested 
an extension of Rayleigh’s theory of cellular 
convection to flowing media. No effects of 
consequence are, however, obtained from an 
introduction of a uniform basic flow into the 
equations. Apparently additional characteristic 
parameters are required. Some authors, (see 
Brunt, 1951) have stated that wind shear is a 
decisive factor. But to the writer’s knowledge 
a systematic investigation of the wind condi- 
tions connected with cloudstreets has never 
been made. 

Under these circumstances an acrological 
and synoptical study was undertaken in the 
fall of 1953 in the Boston area which has covered 
over 50 cloudstreet days so far. Routine ob- 
servations of wind and temperature were used 
in this preliminary investigation. Data included 
pibal soundings of the Weather Bureau Air- 
port Station at Boston and, if available, radio- 
sonde and rawins taken by the Air Force at 
Bedford, a suburb of Boston, and by the 
Army Signal Corps at Ft. Banks in Boston’s 
harbor. When required use was made of other 
New England upper air stations such as Nan- 
tucket (ACK) 90 miles SE, Portland (PMW) 
90 miles NNE, and Mitchell Field (HEM) 175 
miles SW of Boston. 

Consistently a striking feature was found—a 
strong curvature of the wind profile containing 
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Fig. 9. Vertical wind profiles on days with cloudstreet development as observed near Boston, Mass. Windspeed 
in knots (running scale) vs. altitude in feet. Small numbers along curves denote wind direction in tenths of degrees. 
CL. B. = Cloud Base. E = estimated. Cloud photo, Fig. 2, corresponds to sounding at lower right (17 Dec. 


1953). 


in most cases a definite wind maximum within 
the convective layer. There is very little indi- 
cation of a linear wind shear or uniform flow, 
Fig. 9 to Fig. 11. The illustrations also show 
that variations of wind direction with height 
are quite small, in fact, they permit replacing 
the vector wind profile with the scalar wind 
profile without appreciable error. This calls 
for discussion in a later section. It is evident 
from the figures that the mean variation of 
the vertical wind shear with height throughout 
the convective layer is of the order of —r1o-? 
(cm! sec-!) corresponding to a negative cur- 
vature of about 10 m/sec per km?. However, 
in extreme cases, values as high as the threefold 
and as low as one-third of this figure are 
observed. 

With rare exceptions the wind speed in the 
lower atmospheric layer is considerably higher 


than normal. The average scalar surface wind 
on days with cloudstreets in the Boston area 
was about 70 % above the 60 year scalar cli- 
matological mean. (15 knots as compared to an 
annual mean of 9.0 knots and a maximal 
monthly mean of ro.2 knots in March accord- 
ing to the Weather Bureau’s Climatic Sum- 
mary of the U.S.). This situation seems to 
justify the use of the term “wind thermals” 
by glider pilots. 

In some cases with low level cloudstreets 
we find a jet-like wind maximum close to 
the ground (Fig. 10 upper part) which will 
be discussed in connection with the synoptic 
interpretation of the wind profile. Our findings 
on the wind have been summarized in Fig. 25 
(left curve), by averaging height and wind 
velocity for the 3 levels: ground, wind maxi- 
mum and upper wind minimum. (In this 
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Fig. 10. Same as Fig. 9. Upper series shows cases of “low level jets’’. 


figure the forementioned shallow “ground 
jets” have not been included.) 

The earlier observations that cloudstreets 
are topped by a stable layer or inversion are 
borne out by most of our data. On the average, 
this lid was found at about 2 km above ground, 
quite close to the level of the upper wind min- 
imum. See Figs. 9 and 10. 

In a few cases it has been possible to take 
some rough measurements of the cloudstreet 
dimensions by checking cloud shadows with 
land marks in the lowlands around the Blue 
Hills Observatory. For example, on Dec. 16, 
1953, the three nearest cloudstreets had a 
mutual spacing of 7.2 and 8.5 km respectively 
(from axis to axis) while the mean width of 
each cioudstreet was 2.3 km and their length 
more than 40 km. With a cloud top around 
2.5 km the spacing-depth ratio was about 3. 
Evaluation of Fig. 3 indicates an average 
spacing of 8 to 10 km at a layer depth of 2.5 
to 3 km and a length of 50 to 80 km. 
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The actual flow pattern in cloudstreets is 
not sufficiently known. Its measurement repre- 
sents a major problem and the method used 
to explore a tradewind cloudstreet in flight 
(described in Section 6) may give some 
suggestions how it should be done. The quasi- 
steady state of cloudstreets facilitates somewhat 
the study of the airflow in comparison with 
that of the isolated columnar convection for 
which the instantaneous flow pattern is largely 
unknown. The flight report of Steinhoff (Ap- 
pendix) gives some evidence that under the 
base of the cloudstreet the draft distribution is 
not symmetrical with respect to the axis and 
further that there are definite zones of strong 
turbulence imbedded in this draft structure. 
As to the circulation near the cloud top an 
interesting observation by aurM Kamps, 
WEICKMANN and KELty (1953) indicates that 
there is horizontal divergence towards both 
sides of the cloudstreet. These authors write 


(p- 22) 
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Fig. 11. Same as Fig. 9. Observations over Eastern United States. Panoramic cloud photo, Fig. 3, corresponds to 
sounding in upper right (ro April 1957). 


“ 


. we had some dry ice left which we dispersed at a 
rate of about three pounds per mile for eight miles along 
the middle ridge of a cloud line consisting of supercooled 
water clouds. After the seeding was finished nothing 
happened immediately, but nine minutes later ice crystals 
appeared left and right of the cloud line but none under- 
neath. The cloud line itself was not modified. This ex- 
periment indicated that the circulation which formed the 
cloud line consisted of an up draft under the ridge, diverg- 
ing horizontal flow near the top of the cloud line and 
downdrafts at some distance left and right of the cloud 
edges. The velocity in the updraft must have been higher 
than the fall velocity of the generated crystals which was 
about 1% meter per second. Thus, they did not stay in 
the cloud and were carried away horizontally.”’ 


Time lapse pictures of cloudstreets do not 
give a clear impression about the circulation 


as horizontal and vertical motions cannot 
easily be separated visually. Stereoscopic mo- 
tion pictures will be required. 

A study of the synoptic situation reveals that 
most of the extratropical cloudstreets occur 
in vigorous outbreaks of polar air which is 
heated from below during its progress towards 
warmer regions. An intense cold High is 
almost invariably found on the surface map. 
Along the low level trajectories the flow is 
from cold to warm with the cold center to 
the right and the warm center to the left. 
The temperature gradient across the flow calls 
for a pressure gradient decreasing with height. 
By the time the 700 mb level is reached, not 
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during strong cloud street 


development over Boston 
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profile near Boston, see Fig. 
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Fig. 13. Upper air map of 
same day as in Fig. 12. Solid 700 MB CHART 
lines: contours of 700 mb IOOOE 

surface. Dashed lines: tem- 7 JANUARY 1954 


peratures. 
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Fig. 14. Seasonal march of cloudstreet development, as 

observed in Boston (USA) during 1953 and 1954, indi- 

cating that this type of cloudstreet (formed in polar 
outbreaks) is primarily a winter phenomenon. 


much of the surface High is left and, were it 
not for friction, the (geostrophic) wind would 
be a maximum at the surface. 

The temperature gradient along the flow in- 
dicates cold air advection and hence calls for 
backing rather than veering of the wind with 
height. This vertical variation of the wind 
direction tends to line up the vector of the 
upper (geostrophic) flow with that of the lower 
(frictional) flow giving a uniform direction 
throughout the layer. Fig. 12 and 13 give an 
example. 

It can be stated then that the curved wind 
profile owes its existence to a frictional wind 
increase with height in its lower portion and to a 
thermal wind decrease in the upper portion, while 
the almost two-dimensional uniformity of the 
wind direction follows from the cold air 
advection of the low-level flow which gains 
heat during its progress. This situation will be 
encountered again in our discussion of the 
tropical cloudstreets. Occasionally the strong 
flow in the lower troposphere, between the 
deep low and the cold high, is straight over a 
horizontal distance of the order. of 1,000 km. 
The conditions favorable for cloudstreets ex- 
tend into the area where the flow is anticy- 
clonically curved and where subsidence in- 
versions exist with relatively clear skies. This 
confirms earlier conclusions by MALETZkE 
(1953) who considered the Middle West and 
“Texas Northerns” as optimal situations for 
cloudstreets. 
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As a phenomenon of cyclonic activity extra- 
tropical cloudstreets occur preferably in winter 
except in high latitudes. This is demonstrated 
by the (smoothed) frequency curve of Fig. 14 
which, however, should be considered with 
reservations in view of the small sample. In 
the year Oct. 1953 to Sept. 1954, cloud- 
streets over the Boston area were observed on 
35 days. As complete coverage was not ensured, 
a figure of 40 to so days per year may have 
been closer to reality. In other latitudes cloud- 
streets are far more common. 

As to the “low-level jet” profiles of Fig. 10 
mentioned earlier, it is doubtful if they are 
of the same type as the other profiles. In 
some cases these data represent a shallow sea 
breeze with easterly winds. Again we have 
heating of a cold air current flowing over 
warm ground. In other cases the “low-level 
jets” are most pronounced during the night 
and resemble the strongly supergeostrophic 
type discovered during the Great Plains Project. 
(BLACKADAR 1957, LETTAU and DAVIDSON 
1957.) The wind conditions encountered during 
the day may, in this case, represent just what 
is left over from the “night jet”. 

In the latitude of Central Europe early 
spring offers the most favorable conditions for 
the development of cloudstreets when the 
circulation is still strong and solar heating 
already intense. Especially the cold north or 
northeasterly flows, known as “Bise” in the 
Northern Alps resemble the observations in 
the Boston area and have been utilized success- 
fully by numerous glider pilots. For example, 
seventeen French pilots on April 16to 18, 1954, 
flew motorless 500 to 650 km between Paris 
and Bordeaux using cloudstreets. They re- 
peated the same feat on precisely the same 
dates in 1955 when cloudstreets covered 
considerable portions of German and French 
territory. Under similar conditions, the Russian 
Soaring Champion Iltschenko flew his sail- 
plane from Moscow to Stalingrad on May 26, 
1953 over the astonishing distance of 830 km. 
(ILTSCHENKO, 1954.) A good description of a 
500 km glider flight under long cloudstreets 
was given by COMPTON (1957). 


4. Dry Streets 


It is not surprising that the banded structure 
of atmospheric convection is also in existence 
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Fig. 15. Soaring techniques of sea gulls 
over the ocean inresponse to surface wind 
velocity and temperature difference be- 
tween water and air, as observed by A. 
Woodcock (1942). 


WIND SPEED IN METERS PER SECOND 


when the humidity is too low for cloud for- 
mation. The writer well remembers how an 
invisible “Highway of the air” once saved 
him in a desperate situation when the sail- 
plane he was flying had been carried out many 
miles over the open sea by an easterly gale 
over the Swedish west coast. 

Of interest to our problem are Woopcocx’s 
(1942) observations over the oceans. In connec- 
tion with studies on cellular motion Wood- 
cock, in a fascinating piece of work, explored 
the conditions under which the sea gulls 
change their soaring technique from circling 
to straight flight. Fig. 15 gives his results. 
As a first condition for straight flight the ocean 
surface has to be more than 4° C warmer than 
the air above it. Furthermore, Woodoock 
states, “Below a wind velocity (at the surface) 
of 7 meters per second the birds circle about in 
columnar thermals drifting along with the 
wind as they rise... At wind velocities of 7 to 
13 meters per second, birds soar straight to 
windward in a fast flight in which they gain 
altitude and over-the-water distance much 
more rapidly than in circling flight at lower 
wind velocities. While moving against a wind 
of 28 mph birds have been seen to disappear 
high to windward in a few minutes... Both 
the columnar and strip thermals (as described 
by the gulls) start within 100 to 300 feet of 
the sea surface and extend to an unknown 
height. Soaring altitudes of 2,000 feet have 
been measured with a small rangefinder. Gulls 
are seen at much greater heights.” Again we 
seem to have the conditions of “wind ther- 
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mals”. Above 13 meters per second surface 
winds, no soaring was observed by Wood- 
cock. 

In further studies on the band structure of 
the ocean surface Woopcock and WYMAN 
(1947) found the spacing of longitudinal rolls 
about twice the height of the atmospheric 
layer in which mixing takes place. This layer 
is generally isentropic and covered by an 
inversion. The pattern on the ocean surface is 
thought to be due to a concurrent structure of 
the surface wind. The similarity between these 
and our results suggests that the same phen- 
omenon was observed. 


5. Precipitation Bands 


Radar Meteorology has become increasingly 
interested in the band structure of radar echoes 
seen on the PPI scope. (CLEM and Moxon 
1950, AUSTIN and BLACKMER 1956, LEACH 1957,- 
many more interesting observations are hidden 
in scientific progress reports.) 

In identifying these bands with the phen- 
omenon here under discussion one has to exer- 
cise caution for two reasons: 

1. Since radar reveals precipitation rather 
than clouds, important modifications may 
have been introduced into the dynamics of 
normal convection by the action of the precipi- 
tation. 

2. As precipitation falls relative to the air 
and as the radar beam generally intercepts a 
layer of considerable vertical depth, banded 
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Fig. 16. PPI radar scope photo of banded echoes southwest of Boston, 
200 miles. 


the afternoon of 19 Dec. 1955. Total range: 
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as observed from Blue Hill Observatory in 
For wind conditions, see Fig. 18. (Courtesy: D. Atlas 


and A. Chmela, Air Force Cambridge Research Center.) 


radar echoes may be produced by a number 
of different things. 

For example, they may represent the defor- 
mation pattern of virga in a non-uniform 
velocity field or they may depict the “seeding 
path” of an overrunning “generating cell”. 
(Gunn and MARSHALL, 1955.) 

Althot gh only a minority of these bands is 
likely to b: of the type discussed above, radar 
can, in certain cases, disclose meso-scale features 
of convection which otherwise would escape 
observation. This is particularly true when, on 
cold winter days, low-level convection creates 
enough ice crystals (Fig. 2) to produce radar 
Shoes Two cases of this type are illustrated in 
the radar scope pictures of Fig. 16 and 17 which, 
together with the following information were 


kindly supplied by Dr. Atlas and Mr. A. C. 
Chmela of the Air Force Cambridge Research 
Center. On two consecutive days (roth and 
20th Dec. 1955) the northeastern United States 
were under the influence of a severe outbreak of 
Arctic air from Canada where surface tempera- 
tures averaged —20° to — 30°C. North- 
westerly winds carried this very cold air 
through New England and out over the warm 
ocean waters of thé Atlantic where vigorous 
convection developed. Radar bands running 
in the general wind direction (300°) formed 
solithwest of Boston (Fig. 16, 19 Dec. at 17E) 
and on the next day over the ocean waters 
100 miles east of Cape Cod (Fig. 17,20 Dec. at 
14E). Ship reports the morning of Dec. 20th 
indicate surface water temperatures to be 
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Fig. 17. Radar echoes of 
convective bands over the 
Atlantic Ocean, east of 
Boston, in the afternoon 
of 20 Dec. 1955. PPI radar 
scope photo from Blue Hill 
Observatory. Range mark- 
ers: every 25 miles. Total 
range: 200 miles. B = 
Boston. N = Nantucket. 
S = Ship. For wind and 
temperature soundings, see 
Fig. 18. (Courtesy: D. 
Atlas and A. Chmela. Air 
Force Cambridge Research 
Center.) 


more than 9° C warmer than the surface air. 
From considerations of the beam elevation 
angle the tops of the furthest radar echoes can 
be estimated to be not higher than 3 km m. s. 1. 
The Nantucket (ACK) temperature sounding 
of 4 hours earlier and the ship temperatures 
of 7 hours earlier (see Fig. 18) suggestthe same 
value. The length of the bands may well exceed 
the portion (about 100 km) intercepted by the 
radar beam. The large bands in Fig. 16 exceed 


150 km in length. In both pictures the mini- 
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mum spacing between bands (axis to axis) is of 
the order of 8 to to km, the spacing-depth 
ratio being about 3 or a multiple thereof. 

me wind soundings from Nantucket, Fig. 

8 (the position of ACK and the ship is marked 
on the radar photo, Fig. 17) show the familiar 
profile with a negative mean curvature of 
about 10-7 cm} sec-1 in the layer between 
surface and 8,000 feet. 

Here as in all other aerological studies of 
convection over land we are handicapped by 
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Fig. 18. Aerological conditions observed near Nantucket (Ack) during radar band observations of Fig. 16 and 17. 


Left: 
curves denote wind directions in tenths of degrees. Note: 


temperatures (°C 
ship report (marked “S’’ in Fig. 17) gives air temperature 


winds (knots, running scale) vs. height (ft). Right: 


Small numbers along 


) vs. height (ft). 


9°C colder than water surface. 
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the fact that the 12 hourly schedule of inter- 
national upper air soundings provides no data 
from midday for the U.S.A. 

Finally, we may mention the earlier synoptic 
studies of “hail-streaks” and shower bands 
(PROHASKA, discussed in WEICKMANN 1953, 
SCHIRMER 1951). The fact that they run in the 
general direction of the wind over a distance 
of about 100 km suggests possible relations to 
the phenomenon under study here, even allow- 
ing for reservations expressed in the beginning 
of this chapter. The repetition of hail showers 
along such line (apparently from distinctive 
cells within a shower band) is a well-known 
fact. A new case study by Yasupa (1957) of 
precipitation stripes forming in cold front 
showers confirms these observations and estab- 
lishes the connection with a low level jet. 


6. Cloud Streets and “Spiral Bands” 
over the Tropical Oceans 


Cumulus lines stretching in the general wind 
direction are a characteristic feature of convec- 
tion over the tropical oceans occurring speci- 
fically in the trades and in tropical cyclones. 


Tradewinds 


Little has been written about this striking 
cloud formation which is familiar to most 
tropical meteorologists. If the trade wind velo- 
city 1s not too weak the general wind direction 
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Fig. 19. Tradewind cumuli 

organized in bands parallel to 

tradewind direction, as seen 

from 20,000 ft halfway be- 

tween Puerto Rico and Ber- 

muda, 24 July 1955. (Photo: 
J. Kuettner.) 


can again be estimated by a flying observer 
simply from looking down upon the cloud- 
streets. (Fig. 19.) In spite of the somewhat 
stable lapse rates in the upper sub-cloud layer! 
(Markus 1957) the tradewind cumuli must be 
considered as a convective phenomenon which 
accomplishes the distribution of heat from the 
ocean surface through the tradewind layer. 
Accordingly, the mean trajectories of the 
trades cross from low to high potential tem- 
peratures (RIEHL and MALKUS, 1957) in a 
similar manner as the polar outbreaks dis- 
cussed earlier. 

The vertical wind variation in the trades 
(Fig. 20) shows the same characteristics as in 
the extratropical cloudstreets, only on a more 
modest scale. The profile is negatively curved 
with a vertical shear gradient of the order of 
—5 x 10% (cm? sec). “Winds decrease 
downwards and upwards from the cloud base” 
(RIEHL, 1954). Average level of cloud base 
and tradewind inversion are 900 mb and 800 
mb respectively as indicated in Fig. 20. In con- 
trast to the situation in higher latitudes, the 
inversion does not mark a discontinuity in 
wind speed and direction. The nature of this 


1 This sub-cloud layer is also found in convection over 
land (LUDLAM 1953) and is most pronounced in the 
surroundings of cumulus clouds, not directly under 
them (probably as a product of mixing between cloud air 
and environment). It does not keep glider pilots from 
climbing right into the cloud base. 
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Fig. 20. Vertical profiles of tradewinds. Solid curves: Mean profiles after Riehl, 1954 (Pacific Ocean, July to 
October 1945) and Palmen, 1955 (along half of latitude circle 13° N, Pacific-Atlantic, Jan—Feb. 1952). Dashed 


curve: Doppler radar 


measurement by B-29 research aircraft of wind profile connected with cloudstreet, 


Fig. 21, on 24 July 1955. Small numbers along curv.s denote wind direction (in degrees). 


wind profile has been discussed by RIEHL and 
Markus (1957) who also noted the almost two- 
dimensional character of this type of heat 
driven flow with little vertical change in wind 
direction. From a viewpoint of balance of 
forces they postulate “‘a vertical wind profile 
with a curvature in the turbulent layer such 
that turbulence must act to retard the flow... 
this kind of profile is typical of the portions of 
the tradewind belts where the meridional tem- 
perature gradient is directed poleward. The 
geostrophic wind is strongest at the ground and 
decreases throughout the tradewind layer. On 
account of friction, however, the maximum 
wind is situated some distance above the 
ground. . . In the trade section (along a ground 
streamline) it is warmer at the downstream 
than at the upstream end. Since the pressure 
decreases from right to left across the section 
looking downstream, the isobars will rotate 
counterclockwise with height under the in- 
fluence of the temperature field and thus assume 
more and more the direction of the surface 
wind.” This is a very similar situation as the 
one governing the extratropical cloudstreets 
(Section 3). It appears that the uniformity of 
the flow direction and the developing negative 
shear of the geostrophic wind (resulting in the 
curvature of the actual wind profile) are gener- 
al characteristics of heated flows. We will 
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encounter a similar situation during the discus- 
sion of the squall lines. 

In the summer of 1955, R. Rados and the 
writer had an opportunity to measure charac- 
teristic parameters of a tradewind cloudstreet 
during a flight over the subtropical Atlantic 
Ocean with a specially instrumented B-29 of 
the Air Force Cambridge Research Center 
assigned to Hurricane studies (RADOs, 1955). 
The cloudstreet (Fig. 21) stretched about 25 
miles east-westwards in the general wind direc- 
tion and was arbitrarily selected out of a great 
number of similar streets, 300 miles northeast 
of Puerto Rico. Vertical wind and temperature 
profiles were measured by Doppler radar APN- 
82 and Vortex thermometer in the following 
way: First, the aircraft climbed slowly under 
and through the cloudstreet along its axis, sub- 
sequently it descended on a parallel course in 
clear air. The familiar wind characteristics may 
be recognized in Fig. 20. To compare further- 
more the mean vertical motions beneath the 
cloud base and in the clear air between streets 
two “precision speed runs” were made at1,000 
feet height in the direction of the cloudstreet 
axis holding pressure altitude (and power)con- 
stant by means of the autopilot. This technique 
developed during the “Jet Stream-Mountain 
Wave Project” allows one to evaluate differ- 
ential air speeds in terms of approximate ver- 
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tical velocities of the atmosphere. The mean 
air speed beneath the cloudstreet exceeded 
that on a parallel course outside the street by 4 
knots. According to calibrations, this corre- 
sponds toa relative vertical motion of approxi- 
mately 1 meter per second averaged over the 
whole length of the cloudstreet, the updrafts 
being beneath the cloudstreet while the down- 
drafts are outside in clear air. The 1,000 foot 
flight level war about half way between the 
ocean surface and the cloud base. On fur- 
ther runs across the street at 3,000 ft a similar 
technique indicated maximum local updrafts 
of the order of 4 meters per second inside 
portions of the cloud. 

From this crude investigation it appears 
that the vertical circulation in and around this 
tradewind cloud street may have been of 
the type suggested by the seeding experiment 
mentioned earlier (AUF'M KAMPE et al, 1953), 
namely: “An updraft under the ridge, diverg- 
ing horizontal flow near the top of the cloud- 
line and downdrafts at some distance left and 
right of the cloud edges”’. 


Tropical Cyclones 


It is now generally known that spiral bands 
are an essential part of the hurricane circula- 
tion. (Fig. 22.) Probably the first practical 
evaluation and description of these spiral bands 
from radar data was given by H. WEXLER 
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Fig. 21. Tradewind cloud- 
street, about 300 miles north- 
east of Puerto Rico, probed 
by-B-29 research aircraft on 24 
July 1955, see Fig. 20. Cloud- 
base: 2,000 ft, tops: 6,000 ft, 
length 25 miles. 
(Photo: J. Kuettner.) 


(1947). He reported that the spacing of the 
radar bands varies from $ to 25 miles, averaging 
about 10 to 15 miles. These values have been 
confirmed in later observations. (Fig. 23 suggests 
a characteristic spacing of the order of 10 miles 
or a multiple thereof.) The depth of the layer 
affected by shower activity in the banded area 
can be estimated to be of the order of 5 miles 
giving a spacing-depth ratio between 2 and 3 
(which is a familiar value). The length of the 
radar bands usually exceeds 50 miles, occasion- 
ally 200 miles. In between the precipitation 
bands multiple lower cumulus streets can be 
observed, see Fig. 24, which are rarely visible 
on the radar scope. According to verbal infor- 
mation by Mr. R. Simpson who took the 
picture, the cumulus cloudstreets in the lower- 
left of that figure extended through the wind 
maximum below 10,000 feet. The high, tower- 
ing cumulo-nimbus line in the background 
illustrates the negative shear in high levels 
where the cloud tops lean backwards. (The 
picture is taken from the eye outwards.) WEx- 
LER (1947) has already suggested that the radar 
bands form from tropical cloudstreets of this 
type. A confirmation may be found in the 
recent investigations by SENN, HISER and 
Bourret (1957) who find: “The mean crossing 
angle of the radar rain bands with the isobars 
is the same as that of the wind flow at about 
2,500 ft. This suggests strongly that it is the 
wind flow around the 2,500 foot level with 
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Fig. 22. Spiral bands of Hur- 
ricane Connie, 12 August 
1955, as seen On IO cm radar 
at Cape Hatieras. Range 
marks every 20 miles. (Cour- 
tesy: US Weather Bureau.) 
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Fig. 23. Spiral bands of Hurricane Ione 
as observed by airborne radar of B-29 
research aircraft on 17 September 1955. 
Dashed and crossed lines give path of 
two flight traverses. Hatched areas are 
observed radar echoes, solid areas are 
quasi-permanent echoes. Analysis by 
N. La Seur, 1957. (Reference system 
moved with hurricane center.) 
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Fig. 24. Cloudstreets surrounding hurricane as viewed radially 
band in background with tops falling backwards. Cloud systems are circulating from right to left. 
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from center. Note cumulo-nimbus 
(Photo: R. 


outwards 


Simpson.) 


its associated intense convective activity which 
is responsible for both the generation and the 
shape of the spiral bands.” 

In the later stages of hurricanes which have 
moved into higher latitudes the variation of 
wind direction with height complicates the 
picture, but the orientation of the bands still 
agrees with the direction of the winds near 
the surface. (KESSLER and ATLAS, 1956.) 

Very few vertical wind soundings are avail- 
able from the neighborhood of the mature 
hurricane center, although recent horizontal 
soundings by aircraft have supplied a rather 
complete picture of the circulation at a given 
pressure level (Rapos 1955, LA SEUR 1957). 
Until data from the National Hurricane Pro- 
ject become available conclusions will have 
to be based upon present evidence. Fig. 25 
shows wind soundings in close proximity to 


a Pacific Typhoon near Guam (RIEHL 1954) 
and an Atlantic Hurricane near Bermuda 
(JORDAN 1954). Apparently the hurricane pro- 
file is an enlarged image of the wind profile in 
the extratr opical cloudstreets, the shear gradient 
being, again, of the order of —10-7 a sec-1 
(or slightly more) but the wind maximum 
occurring at a somewhat higher level, close to 
10,000 feet. 


The Hurricane Eye 


Some radar bands spiral directly into the 
rain wall around the hurricane eye (Fig. 22; 
RIEHL 1954, Fig. 11.24). At the risk of stretching 
the point one may venture to suggest that 
the rain wall itself is “the ultimate convection 
band” whose diameter is given by the spacing 
between bands (2 or 3 times its height). As 
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Fig. 25. Vertical wind profiles (Us Air 6 
Force) near center of Pacific typhoon 5 
and Atlantic Hurricane. Numbers along 
curves give wind directions in tenths of 
degrees. For comparison, mean wind 3 
profile of extratropical cloudstreets is à 
shown in lower left. Note: curvature 
values are in all cases of similar order 

(ios? cm=lsec-i). © 
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the rain wall extends in the vertical to as much 
as 7 or 8 miles, a typical diameter of the eye 
of 20 miles would be in accord with expec- 
tations. Following this line of thought the eye 
itself would represent the downward leg in 
the vertical circulation. As “the pressure gra- 
dient-force decreases upward... air particles 
moving upward... will be accelerated out- 
ward” (RIEHL 1954) and there should be a 
systematic outward displacement of precipita- 
tion falling from the tilted band. Air descending 
in the center branch of the circulation will 
quickly evaporate the remainder of the liquid 
water entrained from the cloud wall, and sub- 
sequently, will follow the dry rather than the 
moist adiabate. Downdraft and humidity ob- 
servations along the inner cloud wall by La 
SEUR (1957) support this assumption.! 

The situation, if true, has a striking similarity 
with the classical foehn (FICKER and DE RUDDER 
1943). Here, precipitation also falls outside 
the descending branch of the circulation, name- 
ly, on the upwind side of the mountain range. 
In both cases, hot and dry air is produced 


1 MALKUS (1958) in a paper published after completion 
of this manuscript deals more rigorously with this prob- 
lem. 
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throughout a large portion of the atmosphere 
with excellent visibility and very low hydro- 
static pressure (“lee-low’’) and in both cases 
the horizontal pressure gradient between the 
precipitation and the “foehn area” reaches 
excessive values. 


7. Some Related Phenomena 


At this time it may be permissible to extend 
the speculation started in the last section to 
circulations of the atmosphere which contain 
some of the characteristics of band convection, 
although, at first glance they would appear to 
be of an entirely different nature. The two 
phenomena to be considered are the squall 
lines and the “polar bands”. 


Squall Lines 


Recently it has become known that most 
squall lines are connected with a low-level jet 
(PORTER et al 1955, STAFF MEMBERS SELS CEN- 
TER USWB 1956, BAIRD 1956). The fact that 
intense convective activity organizes itself in a 
line parallel to the low-level flow and that the 
jet, by definition, possesses a curvature of the 
wind profile is of significance in connection 
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with our earlier discussions. The flow usually 
appears in the 850 mb map as a band in the 
horizontal plane. Porter et al (1955) found it 
at this level in 75 % of 171 cases studied, while 
BarrD (1956) did so in all of his 12 cases. Interest 
in this jet has generally centered around its 
ability to transport moisture from the gulf 
coast en in a ‘tongue’. If viewed in a 
vertical plane the jet characteristic of the wind 
profile is even more striking. The flow generally 
peaks below 5,000 ft and the curvature reaches 
values in excess of 10-7? cm-! sec-!, occasionally 
approaching the enormous value of 107$ (cor- 
responding to a change of wind shear of 100 
m/sec per km?). 

Recently this type of vertical wind profile 
has attracted wider attention in the problem of 
the atmospheric boundary layer (BLACKADAR 
1957, LETTAU and Davipson 1957). It has been 
interpreted as an inertial oscillation induced 
when the frictional constraint imposed by 
daytime mixing is released at nighttime, an 
effect most pronounced if the geostrophic wind 
is a maximum at the ground. We have seen 
earlier that the latter is a characteristic of 
heated flows. The geography of the South 
Central U.S. favors the entrance of tradewind 
air from the Gulf of Mexico. This air carries 
its (band favoring) properties along (see Sec. 
6, viz. high moisture content, curved win 
profile, tradewind inversion). The low-level 
geostrophic wind maximum may then be 
expected during the time of maximumheating, 
ie. in the afternoon, while the peak wind of 
the inertia oscillation is reached in the early 
morning hours. Full development of squall 
lines and tornado activities generally occurs in 
between these hours, that is, in the late evening 
(TEPPER 1954). 

An example which is exceptionally well 
documented by synoptic radar data (Licpa 
1956) is presented in Fig. 26. At the time 
selected (2100 CST) a squall line of almost 
1,000 miles length reached its full strength 
with the severest weather concentrated in its 
north center section. To give an idea of the 
nature of the three-dimensional wind field 
the figure shows the vertical wind profiles fold- 
ed over onto the horizontal plane with the axis 
of rotation pointing into the direction of the 
low-level jet. 

The solid black arrows in Fig. 26 giving 
speed and direction at the lowest maximum- 
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wind-level, indicate that the orientation of 
the squall line parallels roughly the direction 
of the low-level jet. The vertical wind profiles 
show the strong curvatures connected with 
this type of flow. In the horizontal plane, the 
low-level jet seems to be most pronounced 
along the leading edge of the squall line and 
just ahead of it. In drawing an analogy to the 
cloudstreets discussed earlier we may suspect 
that here, too, the convective activity organizes 
itself in a line parallel to the wind, in accord- 
ance with the curvature of its vertical profile.1 

Of course, there is enough indication that 
other factors, such as the horizontal deforma- 
tion field of the low-level jet stream, the 
moisture distribution, etc., contribute to the 
formation of a squall line, but inspection of 
numerous squall line cases reveals that the 
features displayed in the example of Fig. 26 
are typical enough to be generalized. 


Polar Bands and the Jet Stream 


It has been known for a few years that 
characteristic cirrus clouds develop in the jet 
stream (SCHAEFER 1953). Prominent among 
the forms seen are bands and lenticular clouds. 
Methods have been described for using these 
clouds as a navigational tool on transatlantic 
flights (Frost 1953, 1954). The cirrus bands 
which run in the direction of the maximum 
wind (Fig. 27, 28) are probably identical with 
the “polar bands’ widely discussed in the 
early cloud literature. It appears (Tycos 1927) 
that Alexander von Humboldt first described 
these clouds which seem to radiate from one 
point at the horizon, stretch across the entire 
sky and converge at the opposite side of the 
horizon, actually covering more than 400 
visible miles. Humboldt believed that they 
formed under the influence of the geomagnetic 
field pointing toward the magnetic poles and, 
therefore, called them “polar bands”. At the 
end of the century CLAYTON (1896), in some 
monumental nephanalysis work, studied 1,392 
cases of these bands and found significantly 
“that at all velocities the majority of bands 
lies in the general direction of their motion and 
but few at right angles to it; but the proportion 


1 It should be noted that a secondary wind maximum 
exists near the 10,000 ft level. There are also indications 
of inversions (I) or shallow stable layers (S) developing 
in the late evening near the ground, particularly in the 
rear of the squall line. 
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Fig. 26. Aerological analysis of squall line about 1,000 miles in length, over Central USA on 30 September 

1954. Position of squall line according to radar analysis (LIGDA 1956). Vertical wind profiles foided over into 

horizontal plane give wind velocity (knots) vs. height (feet). Small numbers along curves give wind direction 

in tenths of degrees. For further explanations, see text. Note: orientation of line follows approximately low 
level jet (black arrows). 
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of the bands which move in the direction of 
their length increases rapidly with the veloci- 
ty”. Between 1889 and 1894 nearly 80% of 
the rapidly moving bands (more than 100 
knots) “varied less than 23° from the direction 
of motion”. He also concluded “that cirrus 
bands have the greatest mean velocity when 
moving longitudinally, and the least when 
extending at right angles to their direction of 
motion”. 


Fig. 28. Cirrus band in close 
proximity of and parallel to 
jet stream core of 180 knots 
over Boston, 3 April 1957, 
as seen from Blue Hill Ob- 
servatory, looking northeast. 
(Photo: Conover). For upper 
winds, see Fig. 30. 
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Fig. 27. Typical jet stream 
cirrus, 18 Dec. 1956, looking 
southeast from Blue Hill Ob- 
servatory. (Photo: Conover.) 
Westerly jet stream overhead. 
Note: Northern edge of cirrus 
shield and cirrus bands paral- 
leling jet core on its right 
side (looking downwind). 


He cites earlier studies by Ley in England 
“which give a much higher velocity for the 
bands than for ordinary cirrus”. This and the 
following discussion suggest that the type of 
cloud, now generally called “cirrus radiatus” 
indicates the existence of a jet stream overhead. 

New information on these clouds has been 
obtained in “Project Jet Stream” (McLEAN 
1957) and “Project Cloud Photo” of the Blue 
Hill Observatory (CONOVER 1958). Tentative 
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assessment of this knowledge may be summa- 
rized as follows: 

In a typical westerly jet stream the cloudiness 
consists of a cirrus shield extending to the south 
of the core and frequently subdivided in a 
band-like fashion. The striking feature is the 
sharp cut-off on its northern edge, parallelling 
the jet stream core over distances of many 
hundreds of miles. The space between this 
straight edge and the core just north of it may 
be as small as a few miles or as wide as 30 
miles. It is either clear or filled with cirrus 
bands extending over distances of the order of 
hundreds of miles in wind direction and spaced 
several miles apart. The core itself and the 
area north of it is generally free of clouds. 
Occasionally there may be just one long cirrus 
band parallelling the jet core over tremendous 
distances. On the record breaking coast-to- 
coast flight of the Jet Stream Project’s B-47 (25 
January 1957) the northern edge of the cirrus 
bands parallelled the zone of maximum wind 
(asmeasured by Doppler-radar) more than 1,000 
miles over the Western Plains of the U.S.A. 
Capt. Frost (1953) of BOAC describes how he 
followed with his transatlantic stratocruiser an 
individual cirrus band over a distance of 1,700 
miles closely delineating the jet stream core of 
about 150 knots. Obviously the polar bands are 
not dependent on the terrain below. 

Flignts across the jet stream reveal that the 
cirrus bands form in the area of the tropopause 
“break” (Enprich and MCLEAN, 1957). The 
structure of this “break” is most interesting 
and by no means discontinuous as the name 
may imply. If viewed from an aircraft flying in 
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the jet stream core, the “break” looks as 
follows (see Fig. 29): 

To the north of the jet stream at the height 
of the core, skies are very clear, corresponding 
to the polar stratosphere. The haze layer mark- 
ing the lower side of the tropopause can be 
seen (visually) to rise rapidly but continuously 
from north to south through the core area over 
a distance of about 30 miles. On the south 
side the haze tops the cirrus by several thousand 
feet depicting a transition layer which is sand- 
wiched between the cirrus deck and the sub- 
tropical stratosphere. (This transition layer 
may be connected with the somewhat diffuse 
tropopause generally found on temperature 
soundings near the jet stream core; Type B 
tropopause, ENDLICH, SOLOT, and THUR 1955.) 
Between the core and the edge of the cirrus 
shield to the south the cirrus bands are found. 
Recent observations of “Project Jet Stream” 
reveal that a wide zone of dynamic instability 
spreads south of the jet stream core which 
may be responsible for the development of 
the cirrus shield. (ENDLICH, KUETTNER, MCLEAN, 
and Rapos 1958.) Imbedded in this cloud mass 
are layers of convective instability. 

Turbulence encountered almost invariably 
in jet stream cirrus is indicative of unstable 
zones developing in the upper troposphere 
just to the south of the jet stream core. The 
level of maximum wind falls in this layer only 
on the southern side of the core while on the 
northern side it continues into the stratosphere. 
The whole phenomenon is perceptable only on 
the meso-scale. 

Fig. 28 (3 April 1957) gives a photo of a jet 
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stream cirrus band taken from the Blue Hill 
Observatory near Boston (courtesy of Mr. 
John Conover) about one hour before the 
GMD-2 rawin sounding of Fig. 30 was taken 
at Bedford, Mass., 20 miles away (courtesy of 
Mr. Leviton). The supplementary temperature 
sounding made at the same time, revealing a 
superadiabatic lapse rate below the tropopause 
is not considered reliable. (The height of the 
tropopause was indicated above the peak wind.) 
Within the 5,000 ft thick layer directly below 


GMD - 2 - Sounding 
3 April 1957 
Bedford, Mass 


Fig. 30. Vertical wind sounding through 
jet stream of 3 April 1957 over Boston 
(see Fig. 28). GM D-2 test sounding. 
(Courtesy: Mr. R.Leviton, Air Force 
Cambridge Research Center.) 
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the maximum wind of 180 knots the shear 
gradient reaches 3 x 10°”(cm! sec-1). 

A similar value is given in the jet stream case 
of Fig. 31 (8 March 1954). Here the Jet Stream 
Project’s B-47 reported cirrus streaks between 
33 and 34 thousand feet close to and slightly 
south of the jet stream core which was located 
at 34,000 ft with 116 knots. Temperatures and 
winds were recorded by vortex thermometer 
and Doppler radar at the levels of 30, 34 and 
36 thousand feet. Since all jet streams display 
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thousand feet) by B-47 aircraft of “Project 
Jet Stream”. Florida, 8 March 1954. 
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Fig. 32. Cirrus bands forming 
in the lee of Sierra Nevada, 
California, and trailing down- 
wind under jet stream con- 
ditions during afternoon of 
24 April 1955. 
(Photo: B. Woodward.) 


spectacular values of vertical shear the vertical 
variations of this shear are always remarkable 
and exceed considerably the value of 10” (cm-1 
sec~) encountered in our other studies on cloud- 
streets. 

According to the criteria found earlier the 
optimal location for a possible developmen of 
convective bands would be directly south of 
the core, where the wind maximum lies in 
the upper troposphere thus locating the maxi- 
mum curvature of the wind profile in the layer 
of maximum vertical temperature lapse rate. 

Preferred but not exclusive places of cirrus 
band development are mountain ranges crossed 
by a jet stream. The cirrus bands shown in 
Fig. 32 formed in the lee of the Sierra Nevada 
under a combined jet stream and mountain 
wave condition. The point of origin of the 
bands (which run, again, parallel to the jet 
stream) is located just downwind of the moun- 
tain crest. It is known that this is not only the 
area of maximal updraft but also of maximal 
streamline divergence in the vertical plane 
which tends to steepen the vertical temperature 
lapse rate thus favoring cloud convection. (The 
bottom streamline falls several thousand feet 
along the contour of the mountain slope, while 
the tropopause streamline is lifted upwards as 
much; HOLMBOE and KLIEFORTH 1957.) It 
has been observed during the “Sierra Wave 
Project” that a strong jet stream tends to break 
up the wave cloud (which runs parallel to the 
mountain range) into bands oriented normal to 
the mountain range. These bands may trail 
over hundreds of miles downstream. 

The idea that the jet stream clouds at the 
cirrus level constitute “helical vortices” of 
Banard type convective cells has been advanced 
by CLEM (1955) and SCHAEFER (1957). The 
similarity between the aerological condition of 
these cirrus-bands and the low-level cloudstreets 
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discussed earlier are striking enough to justify 
such analogy although the term “helical” may be 
somewhat misleading. (There is no evidence 
that individual air parcels describe several com- 
plete circulations.) In addition, observations by 
Pan American transatlantic flights and by 
Project Jet Stream’s B-29 and B-47 show that 
the indicated airspeed deviates systematically 
during level flight in the jet stream by amounts 
corresponding to vertical motions of the order 
of one meter per second (ENDLICH, KUETTNER, 
MCLEAN, and Rapos 1958). These motions 
seem to be organized in bands along the jet 
stream axis. The existence of such circulations 
is of particular interest in the problem of high 
level turbulence and for the art of high altitude 
soaring (KUETTNER 1955). 


8. Preliminary Note on the Underlying 
Mechanism 


The observations presented in the foregoing 
sections indicate that most cloud bands line up 
in the flow direction and originate in convec- 
tive layers. The characteristic wind profile in 
the vertical has a jet-like curvature. Not the 
shear but the shear gradient is the striking fea- 
ture. 

This suggests that inertia forces enter the pro- 
blem, similar in type to those responsible for 
the Tollmien-Schlichting waves in the bound- 
ary layer or the Rossby waves in the westerlies. 

In part II, to be published separately, the 
theory of this type of organized convection 
will be discussed. Here only a preliminary note 
on the mechanism thought to be operative in 
the convection bands shall be given. 

As in any normal convection, unstable den- 
sity gradients due to heating from below or 
cooling from above create buoyancy forces. 
In the vertical circulations set up by these 
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forces adjacent parcels will rise and fall. If, 
as in our case, the heated layer is moving 
horizontally and its velocity as well as its ver- 
tical shear (vorticity) is changing with height’, 
the rising and falling parcels will carry their 
vorticity (shear) along. In this way vorticity 
differences will arise between horizontally ad- 
jacent parcels. These differential vorticities exert 
restoring forces (inertia forces) which tend to 
return the parcels to their original levels (Lin 
1955). Hence, the curved velocity profile acts 
similarly to a stable density stratification and 
tends to counteract the circulation set up by 
the existing unstable density stratification. This 
counteraction will be most pronounced in the 
plane in which the curvature is a maximum, 
i.e. in vertical planes along the flow direction. 
Convective circulations in these planes will be 
inhibited. In contrast, cross-sections perpen- 
dicular to the flow do not show any curvature 
and circulations normal to the flow do not 
meet such opposing forces. 

It is clear, then, that convective circulations 
will tend to restrict themselves to planes nor- 
mal to the flow forming uniform lines with 
axis parallel to the flow direction. Thus, the 
curved velocity profile may be said to exert a 
“lineating effect” on convection. This calls 
for “longitudinal rolls” of the type observed 
in our cloud studies and in many laboratory 
experiments with heated boundary layer flow 


1 Without changing the sign of the curvature. 
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Fig. 33. Convection bands 
forming in air layer moving 
from left to right and heated 
from below (Avsec, 1939 
Fig. 43). In this laboratory 
experiment tobacco smoke is 
injected at left and progresses 
towards right in longitudi- 
nal rolls. 


(BENARD 1900, 1927, TERADA 1928, MAL 1931, 
AVSEC 1939). See Fig. 33. 

One may derive this effect quite rigorously 
in terms of the amplification constant a of 
the convective motion. If exponential growth 
with time f of the type exp (of) is permitted, 
the small perturbation theory! to be presented 
later in more detail leads to the following 
approximate expression 


where o = time constant 
d = [24 m?4n?2 with 
I, m, n=wavenumbers in x, 
y, z direction (z vertical) 
g = gravity 
B = o/o with o =density 


u = horizontal velocity (in x- 
direction) of basic flow 


mean value of partial deri- 
vative with respect to z 


mean value of second deri- 
vative with respect to z 
eddy viscosity 

The first term in the bracket is the gravity 
(buoyancy) term, the second the vorticity 
(inertia) term containing the mean curvature 


ee 


* Developed for the case that the Prandtl number equals 
unity and that the medium is incompressible. 
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u of the wind profile and the last term is the 
viscosity term. Maximal amplification o de- 
mands that the wavenumber / vanish, thus 
eliminating the stabilizing effect of the vortic- 
ity term. Now, /=o means an infinite wave- 
length in the flow direction. This calls for 
convection rolls parallel to the wind. With 
the second term in the bracket removed, our 
problem reduces to the classical form of the 
three-dimensional theory of cellular convection 
(RAYIEIGH 1916), the only difference being 
that the well known uncertainty of this theory 
regarding the ratio of cell sides (ic. of the 
wavenumbers / and m) has now been elimi- 
nated. The observed values of spacing between 
the centers of adjacent cloudstreets (2 to 3 
times the layer height) agree with the theoret- 
ical optimum! and with laboratory experi- 
ments. (JEFFREYS 1926, AVSEC 1939). 

This dynamical interpretation appears more 
attractive than purely kinematical considera- 
tions such as the deformation of pre-existing 
convective cells in shear flows. 

The required curvature of the wind profile 
can also be estimated: Band convection will 
be excited only if, in normal (symmetrical) 
convection, the vorticity forces arising from 
the curvature would be strong enough to 
counteract the buoyancy forces to such a 
degree that the convective growth rate were 
significantly reduced. In atmospheric convec- 
tion the typical growth time is of the order of 
thousand seconds, i.e. the combined buoyancy 
and viscosity terms of equation (1) call for an 
amplification constant o of the order of 10°? 
(sec-1). If the order of magnitude of the vortic- 
ity term is comparable to © ice. if 


O(lu’/2d?) ~ O(a) (2) 


three dimensional (symmetrical) convection 
would be suppressed. In this case marginal sta- 
bility requires optimum cell dimensions (Ray- 
LEIGH 1916) such that //2d? + H/6n, where H 
is the height of the convective layer (about 2 
km in our observations). Using this expression 
in (2) the necessary shear gradient 


lu] » 6no/H ~ 10-7 (cm! sec) (3) 


in good agreement with the observed curva- 
tures of the vertical windprofiles. 


1 Depending on the boundary conditions the theoretical 
values range from 2.00 to 2.83. 


Tellus XI (1959), 3 


9. Synoptic Cloud Analysis 


Our study indicates that the banded structure 
of the atmosphere is more common than has 
been generally realized. Especially in the cumu- 
lus and cirrus level cloud bands oriented in 
the direction of the flow are quite frequent. 
If these cloud types can be identified from 
above the flow direction of considerable por- 
tions of the atmosphere may be estimated if 
viewed from high levels such as a satellite 
orbit (Wipcer and Touarr 1957, GLASER 
1958). Some practice based on a more syste- 
matic study of the interrelation between cloud 
structure and synoptic weather (BROOKS 1951) 
may lead to the development of techniques 
permitting the construction of diagnostic 
weather maps from panoramic cloud photos 
taken over vast areas such as the Pacific Ocean 
from which synoptic observations are lacking. 

A very crude estimate leads the writer to 
believe that at least 30 % of the cloud coverage 
existing over the world has a structure suitable 
for this kind of synoptic analysis. 

Initial efforts based on a network of “whole- 
sky” time lapse cameras are now in progress at 
the Harvard Blue Hill Observatory (CONOVER 
1958). Already now the cloud bands discussed 
in this paper can be expected to depict outbreaks 
of polar air, tropical circulations and the 
position of fully developed jet streams. Synop- 
tic cloud analysis may very well develop to a 
diagnostic and prognostic tool as significant 
as synoptic radar analysis (LIGDA 1956). 
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Appendix 


Excerpt from flight report by Dr. Ernst 
STEINHOFF on the first glider flight in cloud- 
streets (30 July 1935). 

“I was approximately 20 miles to the south- 
east of the city of Bayreuth soaring cross- 
country from Western Germany eastward 
towards Czechoslovakia using thermals and 
cumulus updrafts in a strong WN W flow. Sud- 
denly my climb indicator, in a very smooth 
flight, registered 400 to 600 feet per minute 
so that I approached the base ofacloud in status 
nascendi very fast. I decided to go into this 
cloud since the turbulence did not amount to 
much and to stay there for a while circling in- 
side the cloud. After about 10 minutes, I left the 
cloud towards the east to continue my flight. 

“] saw that there was one cloud further south 
of me much larger than the one I had just left. 
I could see this not only by the cloud itself 
but also from the shadows on the ground. 
Therefore I decided to head south (crosswind) 
to this cloud and try to get underneath its 
base. As soon as I came closer to this 
cloud the intensity of the down drafts increased 
rapidly and I lost all the altitude which I had 
gained in the preceding cloud, and much more. 
I must have been 1,500 feet below the base of 
the cloud when I reached the north rim and 
the air became extremely turbulent. In spite 
of this, I decided to go through this down- 
draft field in the hope that I would find rev- 
erse currents soon. This was not the case until 
I was almost to the south rim when I hit as- 
cending air, again extremely smooth. The 
climbing rates were between 200 and 400 feet 
per minute. It took me less than 10 minutes 
to reach the base of this cloud. 

“This time I paid more attention to the clouds 
south and east of me and could observe the 
next cloud to be even larger than the one 
underneath which I was just flying and to 
stretch in its major extension from WNW to 
ESE. The shadow of this cloud was very 
accentuated on the ground and I could observe 
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that the next cloud, almost invisible to me, 
extended even twice as far as the cloud directly 
south of me. The third one, seen by the shadow, 
seemed to extend from horizon to horizon, 
so I decided that if the same pattern of air 
circulation were present in the clouds south of 
me I should try to reach the third one. The 
distances between the clouds were in the order 
of magnitude: of 2 to 5 miles, the altitude 
above sea level between 1.5 and 2 miles. Each 
cloud was covering approximately half of 
shaded and half of the bright sunlit surface 
beneath it, so I felt that the direction of wind 
and elevation of sun must have something to do 
with this cloud pattern. I flew eastward around 
the next cloud and tried to reach the second 
next. The same pattern was repeated. First I 
met increased downdraft, as reported before; 
then, I met a terrible turbulence with in- 
creasingly strong downdraft until just beyond 
the center of the cloud; then, after a short 
brief turbulence, I reached the area of smooth 
ascending air again. 

“Again I had lost approximately 3,000 feet 
just in crossing over from one cloud to the 
other, I repeated my ascent to the cloud base 
again and took off for the next cloud which 
stretched from horizon to horizon. I now 
could see a number of cloud rows south of 
me leaving about as wide a shadow on the 
ground as the sunlit area in between... This 
time I did not circle any more but pointed my 
nose towards 135° and kept on climbing. With 
increasing climbing rate I increased my speed 
until I reached the cloud base. I flew close to 
the south rim of it, finally with a speed of 
200 kilometers per hour and no climbing rate 
indicated any more. Imust have flown approxi- 
mately 30 minutes at that speed when these 
clouds discontinued. With wind velocities up 
to 120 km/h at 12,000 feet, I must have 
covered 160—180 km in 30 minutes. 

“No one of the southern clouds extended 
further east than the one I was flying beneath. 
I was close to 13,000 feet altitude when I left 
this cloud. Within the next 10 to 15 miles 
there was no cloud ahead of me. The air was 
smooth... When I was nearly 1,000 feet above 
the ground I flew straight over some wheat 
and barley fields and searched again for ther- 
mals. Whenever I found one, I sacrificed some 
of my paper sheets to mark the thermal. A 
number of pieces of paper which were lifted 
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up faster than the other pieces pretty well 
outlined the flow pattern of the thermal. This 
was repeated three times. When I did this 
the last time, I saw an airfield below me with 
the white letters ‘BRNO’.”’ 
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NOTE: The flight ended in Briinn, Czecho- 
slovakia, where Steinhoff met three more 
pilots who had used cloudstreets to cover this 
500 km distance, far beyond the world distance 
record of that time. 
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Authors’ Foreword to the Two Following Papers on the Oceanic Thermocline 


Each of the two following papers attempts to provide a theoretical framework for 
explaining the oceanic thermocline, and the associated thermohaline circulation of the 
ocean. They were developed independently, and then the authors exchanged copies of 
their papers. As they stand, the two theories are not compatible; but the results of each 
are similar in certain respects, and apparently resemble the actual ocean insofar as we 
actually know it. Because of the basic differences in the two theories we have decided 
not to try to combine them into a single paper, but to publish them together in the 
same number of this journal, so that they should both appear in the literature simultane- 
ously, and so that the marked contrasts in the basic formulation should be immediately 
evident to the reader. 


Allan Robinson, Henry Stommel, Pierre Welander 


The Oceanic Thermocline 
and the Associated Thermohaline Circulation’ 


By ALLAN ROBINSON? and HENRY STOMMEL 


Woods Hole Oceanographic Institution, Woods Hole, Massachusetts 


(Manuscript received November 15, 1958) 


Abstract 


A study is made of the thermal structure of an ocean which is dynamically geostrophic, 
bounded by an eastern coast, and driven by an imposed surface temperature distribution 
and wind-stress. The heat equation contains vertical diffusion in the form of a virtual 
eddy mixing parameter and non-linear terms of vertical and north-south advection. Com- 
parison of the results with the North Atlantic shows an agreement with the observed 
temperature distribution in the thermocline region, including the deepening of isotherms 
in mid-latitudes. The value of the vertical component of velocity at the bottom of the 


thermocline is predicted, and a numerical value for the eddy-conductivity obtained. 


1. Introduction 


Ninety years ago the question of whether 
the wind-stresses or thermally produced differ- 
ences of density are the predominant cause 
of oceanic circulation was a subject of public 


1 Contribution No. 1003 from the Woods Hole 
Oceanographic Institution. 

2 Present address, Dept. of Physics, Harvard Univer- 
sity. Part of the material presented here has been included 
in a doctoral thesis presented to the Dept. of Physics, 


Harvard University. 
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debate.? Since that time, this central issue 
of physical oceanography has been treated 
with more caution, and is only indirectly 
alluded to in present-day textbooks. That the 


3 The exchange of letters in NATURE during the 
1870's between W. B. Carpenter and James Croll shows 
how futile simple verbal arguments can be in discussing 
such issues. The reader with a morbid interest in fallacious 
verbal theories may find it entertaining to look over the 
work of the English eccentric William Leighton Jordan, 
“The Ocean, its tides and currents and their causes”, 
Longman’s Green and Co., London, 1873. 


296 


vertical integral of mass transport of the ocean 
can probably be attributed to the wind-stress 
acting upon the surface was established by 
SVERDRUP in 1947. This important concept, and 
its further elaboration in the hands of MuNK 
(1950), Hipaka (1950), and others, are described 
in a survey article by STOMMEL (1957). Compu- 
tation of a vertical integral, however, gives no 
idea of the distribution of a property in the 
vertical, whereas the geostrophic calculation of 
currents from observed hydrographic station 
data gives only the vertical distribution of 
the velocity vector, but not its integral. More- 
over, the vertical integration completely cancels 
the contribution of the thermohaline circula- 
tion. Consequently, there was a tendency, 
starting in 1950, to minimize the possible role 
of the thermohaline circulation. 

The first attempts to formulate a model 
containing an active thermal-convective proc- 
ess were 
(x) an important Russian paper by P. S. LINEY- 

KIN (1955) 
(2) an unpublished doctoral dissertation by 
FOFONOFF (1954 
(3) chapter XI in a book on the Gulf Stream, 
written in 1954—5 by STOMMEL (published 
in 1958). 
Studies (1) and (2) had no meridional boundary 
—the key feature of the Sverdrup wind-driven 
model—and so were not directly applicable to 
actual oceanic basins. Study (3) treated only 
the horizontal patterns of flow that would 
ensue in a meridionally bounded basin on 
the assumption of a fixed vertical mass flux at 
mid-depth, and did not treat the important 
vertical heat flux portion of the problem. 

Later, VERONIS and STOMMEL (1957) and 
LINEYKIN (1957) were able to modify study (1) 
so that the important influence of the variation 
of the Coriolis parameter with latitude could 
be incorporated; and obtained expressions for 
the depth of the thermocline which seemed 
more related to reality. On the basis of these 
models STOMMEL (1958) indicated how they 
might be applied to give a rather detailed 
picture of the patterns of flow in the thermo- 
haline circulation of the world-ocean. Never- 
theless, as explained in the survey article (Srom- 
MEL, 1957) neither of these models imitates 
the actual oceanic density distribution very 
closely; hence it was desirable to construct a 
theoretical model of the oceanic thermohaline 
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circulation with an eastern wall to inhibit zonal 
flow, driven by a meridional temperature 
gradient fixed at the surface, and without 
using the totally unrealistic large vertical 
temperature gradient imposed by both Liney- 
KIN and VERONIS and STOMMEL on the models 
mentioned above for the sake of lineariza- 
tion of the convection equation. 

The present paper is an attempt to construct 
such a model, and we believe it is a substantial 
improvement upon previous studies. However, 
it contains—as they also do—a parametric 
treatment of the mixing processes embodied 
in an “eddy thermometric conductivity” param- 
eter x. This parameter is assumed to be uni- 
form and constant over the entire ocean basin 
—in contradistinction to other studies which 
attribute the oceanic thermocline essentially to 
variations in x (for example, DEFANT, 1936; 
Munk and ANDERSEN, 1948). Thus we envisage 
the ocean as being slowly and evenly stirred 
by some physical process which we cannot 
specify, and the necessity for so doing is a 
measure of how far short we are of actually 
being able at present of giving a complete 
physical account of the thermal convective 
circulation of the ocean. 


2. Formulation 


We consider an ocean thermally driven by 
an imposed north-south temperature distri- 
bution on the surface. The constraint of rotation 
is represented by a Coriolis parameter varying 
linearly with latitude. The ocean is taken to 
be infinitely deep and bounded only by an 
east coast. We inquire into the non-linearly 
coupled temperature and velocity fields in 
the steady state. 

Let (x, y, z) represent the east-west, north- 
south, and vertical respectively, with (u, v, w) 
the corresponding velocity components. The 
origin is fixed in the east coast at mid-latitudes. 
The region of interest is x<o, y<y<y, 
z<o. The fluid system to be considered is 
defined explicitly by the following assump- 
tions. 

The equation of state of the fluid varies 
linearly with temperature only. Furthermore, 
the coefficient of thermal expansion, «, is con- 
sidered zero except when coupled with g, 
the acceleration of gravity. In the equations of 
motion we neglect all viscous and inertial 
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terms, balancing by the pressure gradient the 
gravitational force in the vertical and the 
Coriolis accelerations in the horizontal. In 
summary we use purely geostrophic dynamics 
together with the Boussinesq approximation. 
The equations of motion and of mass con- 
tinuity take the form 


op 


Jar - =0 (1) 
à 
feou +5 = 0 (2) 
g (00) (1 -aT) +2 = 0 (3) 


ou i Ov x Ow 
de ay? ae ae 4) 


where f=f)+ fy, the Coriolis parameter, 0, 
is the mean density, p the pressure, and T 
the temperature. 

In the heat equation we assume a balance 
between diffusion in the vertical by a constant 
eddy-conductivity x, and vertical and north- 
south advections, 


ORT oT 
v 


to : dz a 


v ay fo) (5) 


The choice of x will be discussed later. The 


2 can be justified only in terms 
of boundary conditions and anticipated results. 
The z boundary conditions on the temperature 
will be to specify T as a function of y only at 
z=0, and to require that T approach a constant, 
zero, asymptotically as z > — oo. The horizontal 
velocities, both of the same order of magnitude, 
will be largest at the surface and go to zero 
asymptotically. Thus, near the surface, u, v, and 


negiect of u 


TE 
er have their largest values while TL À speci- 


fied as zero. Clearly v rn will dominate over 


CT : ; 
u—. The inclusion of a western boundary 
x 
would, of course, change the above considera- 
tions locally. 
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The fields of primary interest are T and w; 
two equations in these alone can be obtained. 
By cross-differentiation of (1) and (2) and use 
of (4), we obtain 


ie ahem 
92 


py (6) 


Cross-differentiation of (1) and (3) yields 


aw OT. 
dz ag Ei (7) 


Inserting v from (6) into (5) and (7) yields 
the desired equations, 


22T. is oT fowoT A 
dz? dz B dz dy 5 (8) 


ae Pr 


The boundary conditions imposed at z=o 
are T= T,+ T, y, and w=0;as z>-~, T>0, 
and w remains finite, i.e. w — Woo (x, y). The 
boundary condition at x=o is taken, for 
mathematical convenience, as T=o. This 
assures no flow through the plane x =o, but 
is otherwise physically unrealistic. However, 
in no real fluid or ocean would the above 
equations describe the region near such a 
boundary; therefore our detailed results very 
near the coast will not be valid. Anticipat- 
ing carrying the y-dependence parametrically, 
we reserve discussion of this boundary con- 
dition. 


and 


3. The Similarity Transformation 


T and w are specified at z=0 and ap- 
proach arbitrarily close to their asymptotic 
values along some curve z=f(x). In related 
problems it is often found that at a given x, 
the approach of the fields to their asymp- 
totic values depends only on the ratio of the 
actual depth to the value f(x), of the curve 
at this point (GOLDSTEIN, 1938). The formal 
analogy to these boundary layer problems 
occurs because only vertical mixing is allowed, 
so that only higher order derivatives with 
respect to z are present. Formally, we seck 
separated solutions of the form T=G(x)0(é), 
w=H(x)o(£), where £=zF(x), and all func- 
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tions are tacitly functions of y. This transfor- 
mation is also convenient for satisfying the z 
boundary conditions, e.g., at z=o, &=o and 
T becomes a function of x. In the following a 
prime denotes differentiation with respect to €, 
and a subscript y, differentiation with respect 
to y 

Taking x and z derivatives in terms of x 
and é, and inserting these in (8) and (9), we 
obtain 


xF0" — Hl oy — (4) wd, = 0 (10) 


HF?" - (*) | cos Fe] = (x 1) 


The condition for separability of (10) in x and & 
is merely that H be a constant times F; this 
constant can be absorbed into w, so we require 
F=H. For equation (11) to separate all three 
terms must have the same x-dependence. The 
second and third terms require that G’/G be 
a constant times F’/F, or that G be an arbi- 
trary power of F. The first and second terms 
then require that F be some power of x+c 
where c is a constant. 

In summary, under the general transforma- 
tion 


E=z(x+c}), T=(x +c) +1 H(é), 
w = (x +c)kw(é) 


(8) and (9) take the form 


#9" — wh’ — (5 Jo's, =o (12) 


a” — Ge) [(3k+1)9+ REG] =0 (15) 


The constants k and c are determined by the 
particular boundary conditions. To make 
the surface temperature independent of x, 
we choose k= -4. The choice of c=o gives 
E=zxh, Eco for both z->— and x30, 
and the cast coast temperature will automati- 
cally correspond to the bottom temperature. 


4. Order of Magnitude Considerations 


Equations (12) and (13) are still too compli- 
cated to be solved exactly. However, from 
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them we can determine certain gross properties 
of the temperature and velocity fields: their 
general structure, the asymptotic value of w, 
and the scale of € within which all asymptotic 
values are essentially obtained. Moreover, it 
would not be much more meaningful to 
have the details of the exact solutions because 
the mathematical model which we are using is 
itself not detailed, and also because the results 
are to be compared with the averaged observa- 
tions of the real ocean. 

In seeking unknown orders of magnitude, 
it is convenient to introduce known magni- 
tudes in terms of slightly different parameters. 

J0) ag 

Der I per i 
Let 7 = 108y, y(y) 3f2(y) 
n and y are of order unity while e is of order 
1076. We expect x also to be order unity. To 
determine how « affects the amplitude of w 
and scale of £ we write €= e@£ and W=eo, 
We assume that the scale of £ and amplitude of 
W are of order unity, and further that # and 
W are smooth functions of £, i.e. the functions 
and all their derivatives are also of order unity. 

Under these assumptions, (12) and (13) 
become, with a prime now denoting differ- 
entiation with respect to & 


-10, and e(y) 


e2x9" — er WH — et+>yW'G,=0 (14) 
eth W" + el =0 (15) 


Requiring that the dynamical equation (15) 
always retain both terms gives 2a+b=1. In 
the heat equation (14), we notice that both 
advections are necessarily the same order of 
magnitude. Insisting that vertical eddy-diffusion 


be comparable to the advections yields a=b = = 


The parameter ¢ now disappears from (14) 
and (15), which now contain only terms of 
roughly unit magnitude. From the similarity 
transformation we have obtained the depend- 
ence of the asymptotic value of w and the 
scale depth on x, from the above considera- 
tions the order of magnitude and major part 
of the y dependence of these quantities. 


5. Parametric Linearization and Solution 


The mathematical difficulties remaining are 
inherent in the two nonlinear terms of (14). 
To obtain approximate results, we shall replace 
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these actual terms by linear terms, a technique 
which has found frequent use in convection 
problems (Lewis and CARRIER, 1949). In the 
last term of (14) we allow for convection 
with an average temperature gradient instead 
of the actual #,. In the next to the last term 
we allow for convection by an average ver- 
tical velocity instead of the actual W. The 
averages are not to be taken over the full range 
of Z, but only from the surface to a depth L at 
which asymptotic values are obtained arbitra- 
rily closely. For convenience we define L 
explicitly as the depth at which the fields have 
approached to within e-? of their asymptotic 
values. 

At this point the difficulty arises that the 
average values of #, and W are unknown. 
However, we do know that#, ranges from T, 
at the surface to essentially zero at the depth L. 
We therefore crudely take T,/2 as the average 
©, although the true average will vary in y 
with the shape of #. However, we have abso- 
lutely no information about W. Therefore, 
we shall allow the unknown average value of 
W to enter the equations and parameterize 
the solutions. Then, consistently taking the 
average of the parameterized solution actually 
obtained will yield an algebraic equation in 
the average W. Finally, the solution of the 
algebraic equations, in terms of the external 
parameters, is inserted back into the solutions 
for} and W. 

Replacing (14) and rewriting (15), the final 
linear set of equations is 


ge Le W'=o 


2 


W" +09 =0 


(16) 
(17) 
(18) 


where 


I IB: 
Dez W(C) dl 


A defining equation for L in terms of y and 
the external parameters will be immediately 
available from the solutions of (16) and (17). 

Exact solutions of (16) and (17) for # and 
W have been found. These solutions consist 
of exponentials times Airy functions of purely 
complex argument, containing p in an involved 
manner. The solution # and W have the 
desired properties, ic. they can satisfy all the 
boundary conditions. It would be possible but 
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tedious to calculate y from the exact solutions. 
At this point, with the knowledge of the exact 
solutions, and with the lack of real meaning 
of the detailed form of the solutions W and #, 
we again proceed approximately. 

The quantities of primary interest are the 
depth L and the value of W... To find these 
properties of the solutions, we shall use forms 
for W and # which exactly satisfy the boundary 
conditions and exactly satisfy the first integrals 
of the differential equations, i.e. have the exact 
average properties. Most simply, we take 
2 
Ÿ = Bye 5 (19) 


ee) (20) 


Multiplying (16) and (17) by d£, integrating 
from zero to infinity, 


(21) 


= (0) + 985-72 Wo=0 


W'(o)+ [Bde= 0 (22) 


In (21) and (22) the derivatives and integral 
can be substituted from the forms (19) and 
(20). From (20) we also have simply q = ed 

Solving the algebraic equations, we obtain 


L=2 {(2x) [0 - Tıy]}? 
F2 {el 7 + Fc | ai (23) 


B 
which is independent of y, and 
Woo = — By L?/4 (24) 


The above results could of course be improved 
by taking higher moments of eq. (16) and 
(17) and introducing additional constants in- 
to the forms (19) and (20). 

The approximate forms for the temperature 
and vertical velocity are 


af = Vo (y) exp 


le Ee FO (t+ a i; 


300 
apg 2x Ê À 
3 f(y) To + at 


(exp { }) 


6. Finite depth 


w = 0 (y) 


ES 


(26) 


It is apparent that in mid-latitudes the ther- 
mocline is so shallow that it occupies only a 
small fraction of the actual ocean depth, and 
that the assumption of infinite depth is justified 
in making the calculations of 9 and W in 
the thermocline region itself. James Crease 
showed us in private correspondence concern- 
ing an earlier form of the thermocline model 
that this was likely to be the case. But the real 
ocean in fact has a bottom, and if this bottom 
were flat the vertical component of velocity 
W must vanish there. Thus we regard the 
quantity W as being a measure of the vertical 
velocity just beneath the thermocline (that 
is, at infinity so far as the thermal boundary 
layer is concerned). The form of the similarity 
transformation does not permit us to fix the 
bottom at z=constant; instead, it is necessary 
to put the bottom at ¢ = £, where ¢,>L, a 
very mild compromise. Evidently we can 
still use the form of solution for temperature 
given by equation (19) and need modify 
equation (20) only as follows 


W= male) (20°) 


2 

where the quantitiese L à are so small at €=¢, 
that they are truly negligible (~e-!) and we 
can assert that these solutions very closely 
satisfy the boundary conditions W=o0, =o at 
the bottom ¢ = C4. In the region of the thermo- 
cline the quantity 1-¢/¢, is so nearly unity 
that the formal results given in equations (21) 
and (22) and (23) are not changed except that 

VW will vanish alto- 


in (21) the term 


gether unless we restrict this ‘approximate 
form of dy, T}/2 to depths shallower than ©=L, 
and set it equal to zero for &>L. This is quite 
reasonable to do under the circumstances. Thus 
the computation of L and @ is not affected by 
the presence of a bottom at a depth of several 
times the depth of the thermocline and the 
temperature field is essentially unaffected 
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throughout the ocean. The only modification 
of the vertical distribution of the vertical com- 
ponent of velocity W is that instead of ap- 
proaching a constant asymptotic value W at 
great depths it approaches a maximum value 
just beneath the thermocline, which we find 
convenient to call W,, nevertheless because 
this is still infinite depth so far as the ther- 
mocline depth is concerned, and then varies 
linearly from this maximum just under the 
thermocline to zero at the bottom. However, 
so far as horizontal components of velocity 
below the thermocline are concerned, the 
presence of a finite depth has a major effect: 
it gives them a finite amplitude (in the case of 
infinite depth they vanish beneath the ther- 
mocline). They are calculated from equation 
(6) the terms of which no longer tend to zero 
under the thermocline. 

The depth of no meridional motion is at 
the depth at which v=o or at W’=o, hence at 


E == n = Since throughout most of the 


oceans in mid-latitudes 2¢,/L > 5, the depth of 
no meridional motion (or reference level for 
“dynamic calculations”) is > L, or quite 
distinctly at the bottom of the thermocline. 
In the absence of wind-stress at the surface, W 
must vanish at both top and bottom of the 
ocean, so that the vertically integrated value 
of v, by equation (6) must vanish. This is 
consistent with Sverdrup’s earlier wind-driven 
model. The thermohaline circulation is entirely 
an internal mode and does not contribute to 
vertical integrals of horizontal velocity com- 
ponents—but it is important in trying to 
understand the oceanic circulation just the 
same. 


7. Inclusion of wind-stress 


The dynamical equations (1) and (2) are 
non-viscous. For the wind stress to work upon 
the ocean it is necessary to allow viscosity 
to play a role somewhere, and this has generally 
been done by a boundary-layer analysis of the 
form of Ekman. Experience suggests that the 
Ekman viscous boundary layer is confined to 
an upper one hundred meters in the ocean. 
At any point in the ocean the total vertically 
integrated mass transport in the Ekman layer 
is to the right (in Northern hemisphere) of the 
direction of the applied wind-stress, t, and of 
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magnitude r/f. Since both + and f vary as func- 
tions of position over the ocean, the mass 
transport of the Ekman layer also varies from 
place to place. If we form the horizontal 
divergence of the Ekman layer, we find that 
this implies by equation (4) a vertical compo- 
nent of velocity at the bottom of the Ekman 
layer of an amount 


2 (= < dl ty 

ee 
In regions of divergence of the Ekman layer, 
this vertical velocity is directed upward, and 
we need merely adopt this value instead of 
the value W=o in forming the function given 
in equation (20). In regions where the vertical 
component of velocity beneath the Ekman 
layer is downward—as in the centers of the 
mid-latitude oceanic current gyres—the ver- 
tical velocity actually must pass through zero 
as a function of depth before it reaches the 
upward directed value W just beneath the 
thermocline. Since, in forming the parameter 
y we take a vertical average of W it would be 
inappropriate to try to treat this case by 
proceeding as before. 

The two cases are treated separately: 
Case 1: we< 0, W,>o (the Ekman layer is 
divergent; z is positive upward, ¢ is positive 
downward). The forms of the solution (19) 
and (20) are rewritten as follows 


We = 


(20°) 


At €=0, W=W..(1-4)=W. (28) 


hence a < 1. The definition of y is now 
js 
I a " 
o=t [wid -We (: 5) (18") 
0 
The first integrals are now of the form 
— x0 (0) + pÜo + oe (W.- Woo) =0 (21”) 


W’ (0) + Sod Br (22”) 
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and with the indicated substitutions the fol- 


lowing relations obtain 


2% a yT. 
T Pot Woo (1-3) dor Ce WW) =0 
2aW. 
L 


L 
Pond 


The definition of a is given by (28). 

In these equations x, r, T, W., Ö, are 
regarded as known, the quantities a, L, Woo 
are to be found. In this case the value of L is 
given by the cubic 

L3 


= (9, - YTı) - LW.= 2% 


(23”) 


instead of the simple cube root (23). Once 
having obtained this we find 


129, 


Woo = We- (24”) 


Case 2: w,<0,W.>o(the Ekman layer isconverg- 
ent). Here there is a reversal of the vertical 
component of velocity, W, with depth, and 
therefore we must divide the ocean into two 
layers: one from &=o where W=W, to C=h 
where W= 0; and the other from €=h to £ = oo 
where W=W which we treat exactly as be- 
fore. For our purposes, if we define the un- 
known temperature at h as d, we can write 
the temperature and vertical velocity in the 
lower layer as: 


(19”) 
| -7 Ga | 5 
W=W,.\1 -e (20 
The quantity h is also unknown, of course, 
however, we can write the quantities L and 
Wa in terms of Ÿ, by equations (23) and (24). 
The only change here is that we can probably 
drop the advective term yT, in the expression 
for L, because we will incorporate the merid- 
ional advection mostly in the upper layer. 
As we have already seen, this term does not 
play a major quantitative role in the results. 
If the upper layer is not very thick a very 
rough approximation of the vertical velocity 
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is given by a linear relation (see numerical 
computations, STOMMEL (1956)). 
(29) 


W = W(t - Eh) 


The average of Wfrom ¢ =oto¢ =h is obviously 
W./2 and the vertical derivative W,/h; more- 
over the approximate form of 2T/9y in the 
upper layer is Ti, so that we now can insert 
these approximations in the non-linear tem- 
perature equation and obtain (we can obtain a 
more accurate form for W later by solving 
the dynamical equation (22)) 


7 W. r W. 

a ee (30) 
The solution of this equation is of the form 
Wed 


d= Ce” + Cot aT, (7 (31) 


where C, and C, are constants of integration. 
The temperature must be equal to Ÿ, at C=0 
and to Ÿ, at €=h; therefore 
GC ae Ce =, Oo 
kW 
Ce” +C,= 9, - 2yT, 


(32) 
(33) 


Also, the temperature gradients at {&=h must 
equal each other in both layers 
hW. 
oe wen) 2 
Ce 2 Är a : ah 0, T 


(34) 


and so must the vertical derivative of the 
vertical velocity component: 


Wai 
er (35) 


Eliminating the constants of integration, C, 
and C, and also W., we obtain three equa- 
tions in which there are three unknowns h, 
L, and 0: 


(, n oe | 24 f 2yT, 20: " 
W. h L 


+0,+2yT,-0,=0 (36) 
W. Lo, 
Sr (37) 


L=2(2x/9,)* (38) 


One of these is, unfortunately, transcendental 
so that a general solution cannot be obtained, 
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but different methods will be appropriate for 
different values of the parameters, involved. 
Rather than delving into a lot of algebra, we 
would like to point out that when dealing 
with the application of these ideas to the 
ocean, we do not really know x to begin 
with. Instead we form an estimate of h, L, 
and Ÿ:, from the shape of the vertical tempera- 
ture sounding. We obtain a rough estimate of 
W, from charts of wind-transport of the 
surface layers, and then we do not use equations 
(36), (37), (38) in the same way. Equation 
(38) is used to obtain an estimate of x; and 
equations (36) and (37) can be used as a 
consistency check on the applicability of the 
model to the ocean. We know all the quantities 
that go into them, and they ought to be satis- 
fied if the theory is applicable. We can then 
use equation (35) to calculate the vertical 
velocity just beneath the thermocline, Wo. 
A further point: because of the nature of the 
similarity transformation, it is not easy to 
introduce vertical velocity at the surface, w,, 
with an arbitrary x dependence: in fact we 
are restricted to w, distributions of the form 
xls, 


8. Application to a real ocean basin 


The fact that the field of density calculated 
from equation (25) superficially resembles 
that naturally occurring in the ocean is dem- 
onstrated in figures 1 and 2. Figure 1 is a 
dimensionless perspective diagram of the 
solution (25) with an imposed surface tempera- 
ture distribution of the form 


Oo(y) = 20°C - 25°C x y 


where we have introduced the non-dimensional 
lengths x’, y” by the relations x=L,x’, y= 
=L,y, z=Dz‘. L, and Ly, are dimensional 
constants defining the scales of the phenom- 
enon, and chosen so as to make the numbers 
x’, y be in a range near unity. Thus the 
actual formula used in constructing the dia- 
gram was 


z' 
T=(20-2.5n) exp | -H | 


L,= 108 cm 


where 
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| Fig. 2. Schematic temperature field observed in the North Atlantic Ocean. 


The data for the North Atlantic Ocean is 
sketched roughly in figure 2. This diagram is 
based upon (i) the longitudinal profile of 
temperature on the western side of the Atlantic 
presented in Wiist’s (1936) monograph on 
the oceanic stratosphere, (ii) Böhnnecke’s 
(1936) chart of the North Atlantic surface 
temperature for the month of March, and (iii) 
the 1957 CRAWFORD section at 8° North lati- 
tude. The deep temperature of the Atlantic is 
2° C, not o° C, so in comparing isotherms, 
we write the 17°C isotherm in figure 2 as 
Tellus XI (1959), 3 


15 +2, and compare it with the 15° isotherm in 
figure 1. Certain similarities of the two figures 
are, of course, forced, for example we have 
chosen the surface temperature distribution in 
figure I to correspond roughly with that of 
figure 2. Also the vertical scale of figure 1 
was drawn to look like that of figure 2, but 
actually matching determines the parameter 
x. Even allowing for the forced elements of 
similarity in the picture, the remaining por- 
tions correspond surprisingly well. Of partic- 
ular interest is the fact that the isothermal 
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surfaces are deepest in mid-latitudes—a fact 
which hitherto has been attributed to the 
surface convergence of wind-driven surface 
layers. Here of course, in figure 1, there is no 
wind acting; whereas, in the real ocean dia- 
gram, figure 2, we must be on the lookout 
for the effect of the wind which actually does 
act on the North Atlantic. 

Before attempting to account for the effect 
of the wind, first we calculate various quan- 
tities such as the parameter x and the vertical 
component of velocity by fitting the two 
figures ı and 2 as they stand. Thus the real 
Atlantic Ocean is roughly 5,000 km wide 
and 10,000 km long, and taking L,=5,000 
km, L,=10,000 km makes the horizontal 
scales of the two figures agree. Setting 

=10-“sec-1, T,=20°, 1, = =25..10®°C/cm, 
B=2. 10: cm sec) x æ2x 0 LC ge 
= 10% cm/sec. To fit the vertical scales we want 
the value Hz’ equal to 1 to correspond to the 
depth 1,000 meters at y=o and x =-1, 
and we obtain the value of the parameter x: 


*%=10 cm?/sec 


The vertical velocity at the bottom of the 
main thermocline may be written in terms of 
the thermocline thickness z, defined as the 
depth at which the temperature is e-1 of the 
surface temperature locally. This useful relation 


nd) 


is obtained by elimination of (To D 


between equations (25) and (26). 


aBg 
w(o)=(T,+ T. - 2 
( ) ( 0 1Y) Fey)x t 
If the thermocline thickness were 1,000 meters 


(as it appears at first glance to be at y’=0, x= 
-1) the sub-thermocline vertical velocity is 


w(oo)= 53 x 10-5 cm sec-1 


The expression for the vertical velocity in 
terms of thermocline thickness is interesting 
because it does not involve the parameter x 
explicitly. It is very similar to the expression 
used by STOMMEL (1958) in an analysis of the 
abyssal circulation, but using the much more 
primitive model of VERONIS and STOMMEL 
(1957). We believe that both these values, of 
x and w() are an order of magnitude too 
high. The reason is that we are comparing a 
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model in figure 1 with no winds, to a schematic 
representation of the real North Atlantic 
Ocean which does suffer being acted upon by 
winds. In mid-latitudes Montgomery (1936) 
has computed that there ought to be an 
convergent Ekman layer at the surface, directly 
beneath which there is a vertical component of 
velocity downwards, w.= — 5 x 1075 cm sec}. 
We expect therefore that part of the warm 
water layer in the upper portions of figure 2 
really must correspond to the wind-driven 
layer of thickness h, and that below this layer 
the proper thermohaline regime begins. The 
thickness of thermocline which we used before 
was too big. It included both hand z,. A much 
better idea of the rate of exponential decay 
of temperature in the main thermocline a x’ = 
=1, y =o can be obtained from curve NA 
in figure 3. 

Let us proceed by first writing the quantities 
L, h, and W, in terms of the original coordi- 
nates. This is advisable because the dimensions 
of these quantities are rather odd 


= = (e/x)tz, 


h = (elx)!z, 


W.= (xle)tw. 


where 2; is the “thickness” of the thermocline, 
estimated by the shape of the temperature 
curve well below the inflection point at 
d*9/0€2=0 in the temperature curve, and 2, 
is the “depth” of the wind-driven layer which 
we identify in figure 3, for example, as the 
depth of the inflection point. The temperature 
at the inflection point is defined as d,. The in- 
formation we obtain from the North Atlan- 
tic curve in figure 3 is as follows 


= — 3.5 X 108 cm 
— 9.0 X 104 cm 
Oy = 10. © 


Zh 


IR 


These estimates are for a geographical position 
corresponding to 
x= -5x108cm 
KS 
at which 
e 210-8. (°C sec)-1 
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Fig. 3. Temperature vs Depth 

at mid-latitudes in western 

regions of various oceans. 

N = North, S = South, A = 

Atlantic, I = Indian, P = Pa- 
cific. 


and from Montgomery’s 
charts 


wind-transport 


Wea —§ x 10-5 cm sec} 


The quantities for other ocean basins are 
somewhat different. Using these quantities we 
compute 


L 
5 0-44 cr (sec Cr 
W, = 4.0 cmt! sec—2/3 °C1/5 


h=0.72 cm®® (sec C) 18 9, =10 C 


According to the formulation in section 7, 
Case 2, these quantities are not all independent. 
We need only specify two, and the others are 
computable. However, there is a degree of 
uncertainty in identifying any one of these 
quantities from oceanographic data, so it is 
helpful to list them all, and then using the 
relations (36) and (37) inquiring whether they 
are mutually consistent. Thus, from equation 
(37) we have the remarkable relation 
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Inserting the values of W, h and L/2 we 
obtain d, =12.6 °C which is not very different 
from the rough estimate #, = 10° C obtained 
by inspection of figure 3. The above relation 
is indeed a rather fascinating one, because 
the applied surface temperature does not 
appear directly in it. We now calculate x and 
Wa from the equations (38) and (35) and 
obtain 

# 0.85 cm?/sec 

Ws= -2.5 cm: sec"k °C'ls 


Thus the actual upwelling velocity beneath 
the thermocline obtained from by transforma- 
tion to the original coordinates is 


War 3-1-4405 cm. seo 


We believe that the order of magnitude of 
the mixing parameter x and the deep up- 
welling velocity wo. computed from the North 
Atlantic temperature distribution in this man- 
ner is much mote in accord with the estimates 
of the magnitude of mixing processes and of 
the deep-ocean transports, than the numbers 
obtained neglecting the wind. Thus in applying 
these ideas to the actual ocean it is essential to 
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include the wind stress in the model. The two 
causes of oceanic circulation (wind-stress and 
thermohaline process) are connected in a 
basically non-linear fashion and cannot be 
super-posed as additive solutions. Fofonoff 
has already pointed this out for a zonal ocean. 
Of course, in this model, vertical integrals 
of the transports depend on wind-stress alone— 
as has now been known for twelve years. 

The reader can now substitute the various 
constants into the expression for the tempera- 
ture in the upper layer and see that it affords a 
reasonable approximation to the form of that 
actually observed. 

The state of affairs in the North Atlantic 
Ocean is not markedly different from that in 
other oceans. The vertical distribution of 
temperature at 30° latitude near the western 
sides (but east of the western boundary currents) 
of other meridionally bounded oceanic basins 
is shown in figure 3. (NA, North Atlantic; 
SA, South Atlantic; S.I., South Indian; NP, 
North Pacific; SP, South Pacific.) Actually 
the thermocline in the North Atlantic is 
somewhat anomalously deep compared to 
those in other oceans, but the shape ofthemain 
thermoclines themselves are very similar: i.e. 
the value of z; is about the same. The fact that 
the Pacific is twice as wide as the Atlantic is 
of no consequence because width of the ocean 
appears only as a cube root factor. Some of 
the complications in the temperature curves 
shown in figure 3 are due to salinity effects 
on the density field—but these are details 
which at the present stage of the theory do 
not merit immediate consideration. 

The model shown in figure 1 exhibits a 
number of general features which ought to 
be mentioned. First, the warm water in mid 
ocean does not tend to move toward the poles 
—but flows with an equatorward meridional 
component. At low latitudes it turns west 
(since it cannot cross the equator, according to 
equation (6)) where presumably it forms a 
western boundary current when it encounters 
the western boundary of the ocean. The western 
boundary current is a higher-order dynamical 
regime which can occur under fairly general 
conditions of stratification, so there is no reason 
to doubt that the thermal circulation can be 
“closed” on the western side by such a bound- 
ary current—although we have not attempted 
a formal demonstration for this particular mod- 
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el. Everywhere in deep water there is a slow 
upward component of velocity. The way in 
which the solution for a single basin can be 
connected together with solutions for other 
basins to form a scheme descriptive of the 
general circulation of the world ocean has 
already been pointed out by STOMMEL (1958). 


9. Comments on the parametric treatment 

of mixing, the quantity x 

The inclusion of an eddy-conductivity, x, 
is a parametric way of including small scale 
vertical mixing processes into the theoretical 
model. Although this is a very common thing 
to encounter in the earlier oceanographical 
literature, it is a technique which present-day 
oceanographers are very reluctant to employ— 
except as a last resort. It means that we simply 
do not know anything about the physics of 
the mixing process. Despite our ignorance of 
the cause of mixing or the mechanism, we do 
have some idea, as we will show below, of the 
order of magnitude of the mixing. We take x, 
by these other considerations, to be of order 
of magnitude unity (c.g.s.), or five hundred 
times the molecular diftusion coefficient. This 
assumption has important consequences on 
how the model is to be set up, as was shown 
in equation (14). If x were much greater than 
unity, equation (14) would become 9" =o, 
yielding a linear form for # which would have 
no thermocline-like character. If x were much 
less than unity equation (14) would become 
WH —yW'In=0, which describes a purely 
advective process. This equation, together with 
(15), has been solved by the method outlined 
in the beginning of section 4. The solution 
exhibits something like a thermocline, but its 
dependence on latitude does not agree with 
observations. Therefore, the choice for x of 
order unity seems necessary. In the results we 
note that only the cube root of x appears in 
the depth L, so that x may vary by an order 
of magnitude without changing the results 
significantly. 

In the calculations we have taken x to be a 
constant. A more valid x would probably in 
general be a function of depth. Were this 
the case, the first term of the heat equation 
(16) would be replaced by [x(&)9(&)]. In 
the resulting calculation, the first term of (21 
would be replaced by x (0)0(0), and x(o 
would replace x in all the results. This means 
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that within the validity of the above calcula- 
tions the form of x is not critical. If more 
moments of equations (16) were taken, inte- 
gral properties of x(£) would enter. 

e mean distribution of many oceanic 
variables like salinity, oxygen, carbon-14, etc., 
must depend to a large extent on the same 
mean velocity field and parametric mixing as 
those which determine the mean thermal field. 
For the most part these properties have no 
important dynamical effect—except, of course, 
salinity which can play a dominant role, but in 
subtropical areas usually plays a definitely 
secondary role in determining the density 
distribution. Therefore these properties are in a 
sense tracers whose distribution depends upon 
the thermal circulation and the parametric 
mixing, but which does not in turn effect 
them. Thus they provide a means for testing 
our model of the oceanic circulation for inter- 
nal consistency. In order to make such a test 
it does not seem likely that very small scale 
experiments, or transitory phenomena will be 
very convincing because we are so completely 
ignorant about the nature of the mixing process 
that we do not even knowits scale or spectrum, 
or how uniformly it acts in time. Under the 
circumstances we prefer to employ tracer phe- 
nomena whose horizontal and vertical and 
time scales are all comparable to those of the 
oceanic thermocline itself; and we want these 
tracers to be located quite close to or within 
the thermocline. We offer two examples from 
the Atlantic. In both cases salinity is the tracer. 

Case 1: The first case is the layer of Antarctic 
Intermediate Water in the South Atlantic 
Ocean. At Meteor Profile VII in the eastern 
half of the ocean, between Stations 177 and 
180, the salinity minimum lies at a depth of 
700 meters, where the temperature is about 
5° C, thus lying immediately beneath the main 
thermocline. The thickness of the layer is 
between 500 meters and 1,000 meters—depend- 
ing upon how it is defined. The salinity in 
the minimum here is about 0.2 % lower than 
that 200 meters above and below the minimum 
(925/922 > 10° % cm-?). The salinity gradient 
in the east direction is small (ds/dx ~o) as 
compared to that in the northward direction 
(2s/dy ~7 x 10-19 % cm1). Wüst (1956) has 
made dynamical calculations of the velocity 
between stations 177—180 at 700 meters and 
finds a mean northward velocity of as little as 
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v=0.2 cm sec-!. The mixing parameter com- 
puted from this situation is 0.14 cm? sec-1 
[x = (v9s/9x)/(92s/0z?)|. According to our 
theory there ought to be an upward velocity 
of about 2 x 10-4 cm sec! at this level, so that 
we might expect the level of the minimum to 
slope upward toward the north with a slope 
of about 10-4, which indeed it does. 

Case 2: Thesecond case to be discussed is a layer 
of maximum salinity which covers much of the 
southern North Atlantic in the area of the 
North Equatorial Current, and passes westward 
through the entire Caribbean, during which 
passage the intensity of the maximum decreases 
by about 0.3 % (As/dx = + 0.20 x 10-8 %, cm1). 
Upon entering the Caribbean this layer is at a 
level of 125 meters and upon leaving through 
the Yucatan Channel it is 200 meters deep. 
The mean westward velocity of this layer of 
water is about 3 cm sec! (u= - 3.0 cm sec.-1) 
and its salinity maximum is about 0.6 %, 
greater than the salinity 100 meters above or 
below it (925/92?=1.2 x 10°%, cm?) ; thus one 
computes *=0.5 cm? sec}, again a value 
more or less consistent with the requirements 
of our model. This saline layer originates at 
the surface in the subtropical surface salinity 
maximum, and sinks downward as it flows 
southward and westward. The layer is at a 
temperature of 20° C, hence it lies near to the 
top of the thermocline—and the vertical veloc- 
ity should be downward here: this latter is in 
accord with the downward slope of the layer 
as it traverses the Caribbean (slope of 5 x 1075). 


10. The Abyssal Circulation 


As the reader can judge, this model of the 
thermohaline circulation implies many definite 
quantitative things about the circulation of the 
deep waters of the ocean. For example, we 
should be able to derive the total amplitude 
of the deep circulation of the ocean by inte- 
grating the deep upwelling velocity over the 
whole water-covered globe (in spherical coor- 
dinates instead of the beta-plane). In this way 
we should be able to predict the average age of 
the deep waters, a topic of current interest to 
geochemists, and a subject of much present day 
experimental enquiry, and sea-going activity. 
This further extension of the ideas here pre- 
sented is treated in a separate paper by A. B. 
Arons and Henry Stommel which is being 
prepared for publication after this one. 
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Abstract 


A theoretical study is made of the density field and the associated velocity field produced in an 
unlimited ocean by a prescribed density distribution at the surface. It is assumed that all mo- 
tions take place under geostrophic and hydrostatic balance, and that the density is simply advected 
by the motions occurring. The computation is carried out for a spherical earth. The theory gives 
a depth of penetration of the surface disturbances of the order of 1000 m, if one assumes a 
relative density variation of the order ro~* and a characteristic velocity below the boundary 
layer of the order 1 cm. sec~!. The depth of penetration is proportional to the sine of the 
latitude. Assuming a stable ocean with a surface density increasing from the equator to the 
pole the theory gives a meridional distribution of density of the form observed in the real 
oceans. The associated zonal velocities are westerlies at high latitudes, easterlies near the 
equator. 

To permit a more precise check of the theory by laboratory experiments the corresponding 
solution is derived for a rotating “‘dishpan’’. The solution is found to be of the same type as the 
one studied in the spherical case, but it is pointed out that fundamental differences between the 
spherical and parabolic cases are likely to occur in more general solutions than those studied here. 


1. Introduction 


In the last years there have appeared some 
interesting papers discussing theoretically the 
problem of the ocean thermocline (LINEYKIN, 
1955; STOMMEL and VERONIS, 1957). In these 
papers it was assumed that the (turbulent) 
diffusion of density plays an important role, 
and the coefficient of diffusion enters as a key 
parameter in the solution. The density advec- 
tion terms were not considered in their com- 
plete non-linear form but entered the problem 
only as linearized perturbation terms. The 
reason for this simplification is, clearly, the 
mathematical difficulties in handling the non- 
linear effects. 

It should be noticed that the importance of 
diffusion processes in large-scale ocean dynam- 


ics has not yet been proved. It cannot be 
doubted, however, that density advection plays 
a fundamental role so that it seems more 
natural to start out from a purely advective 
model, in which all diffusion effects are neglect- 
ed. In the present paper a study of such an 
advective model is carried out. It appears that 
the model can explain the main. features of the 
ocean density field below a boundary layer of 
thickness 100—200 meters. The present com- 
putation is carried out for a spherical earth. It is 
assumed that the motions take place under 
geostrophic and hydrostatic balance and that 
all motions disappear at sufficiently great 
depths. No vertical boundaries are introduced 
in the model. Accordingly, the model cannot 
predict “Gulf Streams” and similar boundary 
phenomena. 


1The main part of the present study was carried out during the author’s stay at the Johns Hopkins Uni- 


versity, Civil Engineering Department. The work was spondered by the Office of Naval Research and the 
U.S. Weather Bureau under contracts N-onr-248(31) and CWB-9504, respectively. 
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Although the physical model studied here 
is quite simple the mathematical treatment is 
still hampered by the non-linearity introduced 
by the density advection terms. It is possible 
to derive a single differential equation for a 
certain density function, but so far it has not 
been possible to give the general solution of 
this equation. A particular solution is derived 
which seems to yield the essential physical 
results of the model. 

As a next step in a systematic approach to 
the thermocline problem one would include 
density diffusion, at least in the boundary layer. 
If one could solve the appropriate boundary 
layer equations and match the solution to 
that for the deep water region given by the 
advective model, it would be possible to deter- 
mine the density field in the whole ocean in 
terms of the primary forcing functions acting 
at the sea surface, such as heating and cooling, 
evaporation and precipitation. One would 
also like to investigate such factors as vertical 
boundaries, finite depth, wind-produced con- 
vergence in the Ekman layer etc. 

To test the reality of the theoretical solution 
presented here one would like to set up a corres- 
ponding laboratory experiment, using a ro- 
tating ‘‘dish-pan”. One advantage of the labo- 
ratory experiment would be the possibility of 
controlling the effects of the density diffu- 
sion. In the laboratory one could work in the 
regime of laminar motion where only molec- 
ular diffusion is important, and one would 
then avoid the kind of guess-work that enters 
into all estimates of turbulent diffusion effects. 

To permit quantitative comparisons be- 
tween the theory and such laboratory experi- 
ments the present analysis has been extended 
to the case where the equipotential surfaces are 
not spherical but paraboloidal. In this case, 
there occurs in the equation a new term ex- 
pressing the effect of the latitudinal variation 
of effective gravity. In the particular solution 
studied in the present paper this term changes 
nothing essentially, but it seems likely that it 
can lead to new types of solutions in more 
general cases. 


2. Formulation of the problem 


Consider an ocean of incompressible fluid 
covering the spherical earth. The earth is 
assumed to have a radius R and an angular 
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velocity 2. The sum g of the acceleration of 
gravity and the centrifugal force is assumed 
to be constant in magnitude and directed 
radially. The ocean has a depth h and follows 
the earth in its rotation except for small mo- 
tions produced by surface disturbances. To 
arrive at our model, which is characterized 
by a geostrophic, hydrostatic and advective 
balance, the following specific assumptions 
have to be made: 

(a) The horizontal scale of the disturbances 
is much larger than the depth h. 

" (b) The non-linear acceleration terms are 
much smaller than the Coriolis acceleration. 

(c) The boundary layer within which friction 
and diffusion play an important role has a 
depth much smaller than the depth h. 

(d) The Coriolis acceleration is small com- 
pared to the product of gravity and the relative 
density variation. 

If the horizontal scale of the disturbances is of 
the order L, if the density variations are of the 
order Ao, if the characteristic horizontal veloci- 
ties are of the order U, and if the fluid is 
characterized by its mean density, @, and the 
coefficients of (turbulent) kinematic viscosity 
and diffusion » and x, then the above four 
conditions can be written 


te I (1a) 
U 
Lo<’ (1b) 
I iy I (41 
kV. Qala Vall Re 
Uae 
gd 


It should be noticed that our model cannot 
be applied very close to the equator; in this 
region non-linear accelerations and friction 
must come into play. 

With the above assumptions the equations of 
motion, the continuity equation, and the den- 
sity transport equation, valid in the interior 
of the fluid are 


I dp 
rcosO dd . (2) 
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I op 
22 sin Oou = 0 (3) 
ER 
ee (4) 


i a I. * à 
rcos O ad ul r cos Ba nn 
2 


0 2 
+5 (qu) + -qw=o (5) 
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Here r is the radial distance, & is a longitude 
coordinate and © a latitude coordinate, counted 
positive eastwards and northwards, respect- 
ively. u, v, and w are the corresponding velocity 
components, and p is the pressure. 

The above equations will be simplified fur- 
ther by making the realistic assumption that 
the depth h of the ocean is much smaller than 


x 2 ‚ ar 

R. In this case the term — ow in the continuity 
r 

equation can be neglected in comparison with 

the term — 


or 
replace r by R and dr by dz, z being a vertical 


(ow), and one can everywhere 


coordinate which is zero at the surface (strictly 
at a horizontal surface just below the boundary 
layer). After multiplying the density transport 
equation by @ we have then 


; ae, ‚ 

- 29 sin © (9) = - 68 3g 123) 
2Q sin O (ou) = =e (3°) 
te 
som (4) 


Be ee 
Cire peeing hs. - °°) 


Tellus XI (1959), 3 


2m 


Our five dependent variables are ou, ov, ew, 0 
and p. 
The boundary conditions are 


@=00(¢,) at z=0 (7) 
QU=QV=pW=O at z=-» (8) 


where 0, (8, ©) is a prescribed function. The 
condition (8) implies that @ approaches a con- 
stant value, 0, at z= — oo. Since the physically 
interesting solutions have a stable stratification, 
we have @9< 0. 

In the beginning the ocean was assumed to 
have a finite depth h so that in a strict sense the 
boundary condition (8) should be applied at 
the point z= -h. In the following we will, 
however, only deal with the case where the 
depth of penetration of the surface disturbances 
is small compared to h. One may then in- 
troduce the boundary condition in the stated 
form for the sake of analytical convenience. 

It is to be noted that no condition is imposed 
on the vertical velocity at z =o. In fact, a 
condition that w =o at this surface would mean 
an overspecification of our problem, when we 
also require that all velocities should vanish at 
great depths. In a complete solution there 
should enter also a boundary layer on top 
within which friction and diffusion become 
important. With this boundary layer all the 
boundary conditions can be satisfied. However, 
we are here satisfied to study only the deep 
water solution. Accordingly, z=o will not 
represent the real free surface of the ocean 
but rather a horizontal surface just below the 
boundary layer. 


3. Derivation of the M-equation 


From the system of equations (2’)—(6’) one 
can derive a single differential equation for the 
quantity 


M($, 0, 2)=f fe'(dz? () 


where 0'=0-06 is the perturbation density. 
It is assumed that 9’ is given at z=o and that 
it decays to zero for large negative z rapidly 
enough to secure the convergence of the 
integral (9). 
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Eliminating the pressure in equations (2’) and 
(3’) by means of the hydrostatic equation (4’), 
and using the condition (8) one finds 


a £ 90 
eu ZOR 5 | À = 
“= é (10) 
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The continuity equation yields 


nae 8 20 
PT ~ FOR? sin? of {x ae) 


=00' 10,5) 


00 90 À ' 
Replacing 6 by 30 ott and introducing the 


function M one has 


af Le. 
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SOR ne cos © Ms: 132) 
w= DR nr MO 
where Mo = =. etc.; furthermore 
= M¢-z 
2 = Moz: (13) 
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Introducing these expressions into the density 
transport equation (6) we get the desired 
equation for M: | 


M5 z:Mez — MozzMgz= cot O MgMzzz (14) 
This equation can also be written in the form 


(14) 


J (Mzz, Mz) = cot 0My M: 
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where J denotes the Jacobian with respect to 


& and O1. 


The boundary conditions for M take on the 
form 
M.:= 0($, 0) at z=0 | (15) 
M=o at z= -© 


where 0, (¢, ©) is the given surface distribu- 
tion of the perturbation density. 


4. Derivation of a particular solution 


To avoid unnecessary mathematical com- 
plications we confine our attention only to 


the half-sphere 0 <O < = Introducing a new 


latitude coordinate 


(16) 


which runs from O at the pole to — co at the 
equator, equation (14) takes on the simpler 
form 


n = log sin © 


Mg zzMnz — MnzzMgz= Mg Mzzz (17) 
One class of solutions which satisfies equation 
(17) is obviously 


M=M(n, 2) (18) 
In this case, in fact, both sides of equation (17) 
vanish separately. The corresponding velocity 
field is purely zonal, v=w=o. Since there are 
no vertical velocities the density at the bounday 
is not advected into the fluid and the solution 
has little physical interest. For the general case 
where the surface density distribution depends 
both on ¢ and n, one may try a solution of 
the form 


M=P(¢,n) Q(n, 2) (19) 


It should be noted that the present M-equation could 
also have been derived in a so-called ß-plane, using 
rectangular »quasi-coordinates» x, y, 2 (cf. the derivation 
in the appendix). In the ß-plane the equation becomes 


MxzzMyz — MyzzMaz= 5 M,,Mzzz 


Putting dx=RcosOd¢, dy=RdO, f=2Rsin®, = 


ny 


cos © we get again equation (14). 
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It is immediately found that (19) is a solution 


for arbitrary functions P (8, 7), provided Q 
satisfies the equation 


Q2zQnz = COQ. = OQ... (20) 


Suitable solutions to (20) are found by putting 


al) 


Equation (20) changes then into the form 


(21) 


z De 2 7 Inn F(n) 77 
ee ae 2 er — 4 — 
Fo) OO PTT OO Fin) COR 
and it is seen that (21) is a possible solution if 
we choose ae =, OF 


us 

F(n) = const e* (22) 

« being an arbitrary constant. Q is now 

determined by the ordinary differential equa- 
tion 

@ [Q dr = 


(Or Q'Q’ — QQ” =0-(23) 


where & 2 he the argument of Q. 


F 
Particular solutions of (23) of the form Q= 
const 4 are easily found, with possible expo- 


2% 
nent values a=o, a=1, and a > but these 
x 


solutions are not sufficiently well-behaved at 
z=o and z=-o to be of use. For «=1, 
however, the equation has the particular solu- 


tion 
(24) 


where k is an arbitrary constant, and this 
solution represents a possible physical situation 
If k is positive Q will vanish at z= — © and 
will certainly be well-behaved at z =o. 

In the following, the discussion will be 
confined to the solution (24), and we will 
demonstrate that it is capable of describing 
correctly the main features of the ocean 
density field. Nevertheless, it would be of 
great interest to study more complicated 
solutions of the equation (23), which possibly 
could describe finer details in the ocean density 
field. There seems to be little hope to find 
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the general solution to (23), although it actually 
can be reduced to a first order equation}, but 
one can, of course, resort to numerical methods. 

Using the particular solution (24), M takes 
on the form 


M= P(4, ne = Pl, n) ete” (25) 


or, in terms of the original coordinates ¢, ©, z, 


kz 
M=M,(6, 0) ane 


(26) 


M,(¢, ©) is determined by the boundary 
condition at the surface, while the constant k, 
which is simply a scale-factor for z, is at our 
disposal. It can be fixed by prescribing some 
more physical parameter in the problem, as 
example the total energy or the total angular 
momentum of the system. 

To estimate the order of magnitude of k in 
the real oceans we can relate it to the charac- 
teristic horizontal velocity. Using the pre- 
viously derived expressions for the perturbation 
density and the mass-velocities one finds 
easily that the order of magnitude of k is given 
by 


(27) 


where, again, L is the horizontal scale of the 
disturbances, U is the characteristic horizontal 
velocity, and Ao is the characteristic density 
variation. The corresponding depth of penetra- 
tion of the disturbances is of the order 


1 Introducing, in succession, the following new depend- 
LA 


Q { Vb BL Ct= 


ent and independent variables: Z= 


In €, equation (23) transforms into the second order 
equation 
(Z* —a) (2Z* -3Z* + Z*) + (Z* = 30) (- Z*2+ 
+Z*- Z*’) + (Io) Z*’+ (2Z* — Z*’)- 
*(- Z*+ Z)=o 


Introducing V=Z*’ as new dependent and x=Z* as 
new independent variable this equation is further reduced 
to the first order equation 


dV 


a) Va - V?2 + (x? - 30x + 3a)V - x[ax? + 
x 


(& 


(1 -30)x + 20] = 0 
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UDO 
0 


The neglect of the effect of the finite depth 


requires that 7 <I, or 


Ü 2 

Ao 

ZEN 
Vs“ 


This inequality states simply that the square 
of the internal Froude number is much smaller 
than the Rossby number. 


OU 
OL 


< (29) 


5. Numerical estimate of the characteristic 
parameters 


Assuming h=s5-10° cm, L=10° cm, U=1 
A 

em sec, P=10°cm sec, ze 103, 2=10-* 

sec, »=102 cm? sec l, x=102 cm sec-}, one 

arrives at the following values for the five 

characteristic parameters entering in the condi- 


tion (1): 
| a: 
rs 100%; 107197 le 102% 
«L 
- a 104, = 21072 


Furthermore, the internal Froude number en- 
tering in the condition (29) takes on the value 


TJ © 1.4: 1078 
veer 
0 


One can, of course, always argue about the 
precise numerical values to be chosen 


MR A 
for such quantities as L, U, ER y and x, but 


unless the values are definitely wrong our 
assumptions should on the whole be justified. 
The assumption that the diffusion effects are 
limited to a thin boundary layer is certainly 
not justified by the present estimate, since 


Lge 
the parameter wo takes on a value of the 
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order 1, but it should be noted that this value 
is obtained by using a characteristic velocity 
for the deep sea. In the boundary layer itself 
the velocities must be one to two orders of 
magnitude larger and a more satisfactory value 
of the boundary layer thickness is then ob- 
tained. 

With the above numerical values the depth 
of penetration of the surface disturbances 
according to (28) becomes 


D - 1000 m 


This value certainly is of the correct order of 
magnitude (cf. Fig. 1). 


Fig. 1. Example of the vertical density variation in a real 
ocean. (Average value at 22° S in South Atlantic, com- 
puted from DEFANT and WUST, 1936). 


6. The density field 


With the solution (26) for M the density 
perturbation field becomes 


kz 


o'= 06 (¢, 6) en? 


(30) 


where 0 (¢, ©) is the surface distribution. 
According to (30) the perturbation density 
should decay exponentially downwards from 
the boundary layer, and such a law will 
certainly fit the oceanographic observations 
well (Fig. 1). The depth of penetration should 
increase from the equator to the pole, varying 
proportional to the sine of the latitude. This 
result seems also to be in qualitative agreement 
with the observations. 


Tellus XI (1959), 3 


AN ADVECTIVE MODEL OF THE OCEAN THERMOCLINE 315 
(0) 
Equator 
Oo 
0.6 05 0.4 0 
309 0 (+) 
2 40° 50° 60 Isolines En 
Ap 
Fig. 2. Theoretical density distribution in a meridional plane. Surface density distribution: 0’9=—JA@ cos? ©. 


Vertical coordinate: kz. 


From (30) the latitudinal variation of 0’ is 
found to be 
kz 


1 : kz in © 
= 156 Qo cotO— g e 


Ay aU 
According to this formula => has the sign of 


206 E R : ; . 
aA at high latitudes, while at low latitudes it 
has the sign of 06, at points below the sur- 
face. Ina realistic case one would have 9) <o, 
do LA 0 [2 
a >o, and en 
latitudes but negative at low latitudes, at 
points below the surface. Thus the density 
isosurfaces should lie at maximum depth at 


should be positive at high 


middle latitudes and rise towards the surface 
both at the low and the high latitudes. This 
result seems also to be supported by the obser- 
vations. 

As a numerical example the meridional 
density distribution corresponding to the sur- 
face distribution 09= -Aopcos? © has been 
constructed (Fig. 2). In this case 0, runs from a 
negative value —Ao at the equator to a value o 
at the pole, in approximate agreement with 
conditions in the real oceans. (One should 
assume here that Ao shows some longitudinal 
variation so that the degenerated axially sym- 
metric case is avoided.) For comparison we 
show the mean meridional density distribution 
in the South Atlantic, as computed from the 
“Meteor”-data (Fig. 3). 
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Fig. 3. Example of the meridional density distribution in a real ocean. (Average values for South Atlantic, computed 
from DEFANT and WUST, 1936). 
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Fig. 4. Theoretical distribution of the zonal velocities in a meridional plane. Surface density distribution: 
0 o=—A@ cos? ©. Vertical coordinate: kz. 


7. The velocity field 


The velocity field corresponding to the 
solution (26) is given by 


kz 
kz se 
ou = (ue fr - 2 (6. 0). € 
kz 
QUE (ov) oe*” (3 1) 
kz 
ow = (ow),e"° 
where 
1 
ve 
1 000 
I+tg 0120 


and the subscript o denotes surface values. In 
terms of the surface density the surface veloci- 
ties become 


FREU N. 
(eu) = SORE ce A 


RER IE 

(ero 2QRk cos O 90 (32) 
iin Semi ee 

(ew)o 2QR2k?2 Id 


From the above relations one can conclude: 
(a) ov and ew decay exponentially down- 
wards, while ou may show a maximum below 


the surface. Such a maximum is found when 


d , 
a> I, i.e. when —cot 9<— 2% 0. In the 
05 90 


! 


20 : 
— > O a maximum 


210) 


will always occur at low latitudes. 


realistic case where 09 < 0, 


(b) (ov), and (ow), are proportional to 2% 


Qo 
ag’ 
while (ou), depends both on _ 
C 
de 
77] 


the realistic case og <o,—> 0, one would find 


and 0. In 
oO 
(ou)5=0 (westerlies) at high latitudes, and 
(ou). <o (easterlies) at low latitudes. 

A computation of the zonal velocity field 
has been carried out for the specific case 0, = 
= - Aecos*@, and the result is seen in Fig. 4. 

The prediction by the theory of a prevailing 
westerly surface flow at the high latitudes and 
easterly surface flow at the low latitudes is cer- 
tainly in general agreement with observations. 
We cannot, of course, hope to depict the 
detailed current systems in the oceans by the 
present crude model, and we do not give any 
reference here to actual velocity data. 


APPENDIX 
Solution for the dish-pan 


In order to make possible a numerical com- 
parison between the theory and experiments 
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under laboratory conditions we derive here the 
M-equation for a fluid with parabolic equi- 
potentials. In this case the apparent gravity 
varies with “latitude” and some new effects 
can be expected. 

To be able to make a direct comparison 
with the spherical case studied earlier we in- 
troduce similar coordinates: a longitude coordi- 
nate ¢, a latitude coordinate ©, and the vertical 
height z (Fig. 5). In the derivation of the M- 


Fig. 5. Coordinates in the paraboloidal system. 


equation one may, however, work more 
practically with “quasi-coordinates” x, y de- 
fined by the relations 


dx =rdd | ae 
dy = RdO 
where r is the normal distance to the central 
axis and R is the radius of curvature of the 
paraboloid. r and R are both functions of ©. 
Introducing the velocity components u, v, w 
in the x, y and z-directions, the perturbation 
pressure p’ and the perturbation density 0’, 
the Coriolis parameter f=f(y) =2QsinO, and 
the apparent acceleration of gravity g* =g*(y) 


RE a 
=D: and applying the same physical simpli 


fications as are expressed in the spherical case by 
the conditions (1a)—(1d) one derives the sys- 
tem of equations 


(34) 
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fly) (eu) = + (35) 
o=-L (pe (36) 


à 0 
je let) +5, 0”) + em) +o (37) 


(qn) 52 + (v) 32+ (gu) 32 = 0 (38) 


Expressing ou, ov, and ow in terms of 0’ one 


finds 
wey) rao! mat): is 2 
; 9‘ ORK ‘ 


Harrer] 
: To J 2x 


(39) 


where 
1 po 
or =" 


ov and ew are seen to have the same form as 
in the spherical case, where the acceleration of 
gravity is constant, while ou contains an extra: 
term proportional to the vertically integrated. 
perturbation density. Introducing the variable: 
M defined by (9) in Sec. 3 one is led to the: 
equation 


M, M: = MEME == P M; Vie ia 


f 


= m M; Myzz (40), 
It is easily shown that for a TEE = — Le 


Thus equation (40) is reduced to the form 


My22Myz 7 Me M: > 4 (M M: Ir M; Myzz) 


à 


or, after introducing the coordinates and ©, 
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My zz Mo: - Mozz Mgz= (41) 
= cot O(Mg Me + M;M 22) 


This form should be compared with the 
equation (14) in the spherical case. 

Introducing the new variable n=log sin © 
(41) simplifies to 


M zz Mnz rs MyzzM gz = Mg Maz +M; M2 


(42) 


Attempting a solution of the form M= 
P(P, 1) Q (n, 2) gives 


Qxz Qnz - Qnez Q; = QQ-2z ag Q; Qu (4 3) 
and again solutions of the form Q = ann) 


n 
exist provided F(n) = const. ex. Q is then to be 
found from the ordinary differential equation 


£[Q’ Q” - (Q’)?| + (1 - «)Q’Q’-«QQ”=o 
(44) 
z 
where C= En) 


In this case an exponential solution of the form 
D I 

Q=const ek exists if one chooses «=-. 
2 


Thus we get 


M= P(G, 1) e% = P (by n) ete ™ 


PIERRE WELANDER 


or, in terms of the original coordinates, 


The solution for the dish-pan is accordingly of 
essentially the same type as for the sphere, the 
only difference being that the depth of penetra- 
tion is now proportional to the square of the 
sine of the “latitude”. This cannot, however 
be expected in more general cases, due to the 
occurrence of the new term M;Mgzz in the 
equation (41). 
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Abstract 


The very strong effects in the auroral zone of the solar flares of 7 July 1958 as observed at Ki- 
runa Geophysical Observatory by means of magnetometers, an ionospheric sounder, a cosmic 
noise absorption receiver (riometer), oblique auroral reflection receivers, transpolar communi- 
cations receivers, and cosmic ray telescopes are reported and discussed. Several remarkable 
features of the terrestrial disturbances were observed: 

1. Extremely strong absorption became apparent a few hours after the solar flare. In spite 
of a linearly increasing absorption during the first seven hours after the flare no change in 
height or critical frequency of the F2 layer was noted during this period. 

2. The SIDs reported by Pacific Observatories at the time of the flares were not observed at 
Kiruna although Kiruna was on the sunlit side of the earth. 

3. A magnetic storm and a large decrease in the counting rate of the meson component of 
cosmic radiation appeared simultaneously 31 hours after the flare. 

4. The maximum absorption at 27.6 Mc recorded during this period surmounted 20 decibels. 


Introduction 


The ionospheric black-out phenomena fol- 
lowing the solar flares of 7 July 1958 at 0032 
and 0040 U.T. had in Kiruna a longer duration 
(c. 73 hours) than the black-out beginning on 
23 February 1956 (c. 36 hours). On Kiruna 
records it is surpassed with respect to duration 
only by the black-out period of 22 June— 
26 June 1957 (c. 81 hours) during the last 
sunspot maximum. 

The geomagnetic storm, which began 31 
hours after the flares were first observed, is 
with respect to amplitude the strongest storm 
of this maximum period of solar activity. De- 


1 Kiruna Geophysical Observatory. 
2 Electronics Research Directorate, Air Force Cam- 
bridge Research Center. 
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viations from the quiet day conditions in the 
horizontal component amounted to as much as 
2,200 gammas between positive and negative 
peaks. 

On 7 July the sun never dipped below the 
horizon as seen from Kiruna. Its angle of 
elevation was at local midnight between 6 and7 
July only % degree; a refraction correction of 
approximately % degree should be added to 
this to obtain the apparent elevation as seen 
from the ground. At the time of the flares 
the solar altitude was 3 degrees. Thus, Kiruna 
was situated near the boundary of the daylight 
side of the earth. 

During the time of the flare outburst and 
throughout the entire period, equipments re- 
cording six parameters correlated with solar 
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conditions were being operated at the Kiruna 
Observatory. These are magnetometers for all 
three components, panoramic ionospheric re- 
corder (ionosounder), a unit for measuring 
ionospheric absorption of cosmic noise at 
27.6 Mc/s (riometer), a VHF receiver for 
recording 87.9 Mc/s oblique auroral reflections, 
receivers for recording transpolar transmissions 
from College, Alaska, on 11.805 Mc/s and 
17.900 Mc/s and cosmic-ray meson telescopes. 

It might be of some interest to examine in 
detail the effects of the solar flare on the param- 
eters mentioned for a far northern location 
(geographic coord. 67.8° N, 20.4° E; geomag- 
netic coord. 65.3° N, 115.5° E) situated during 
the first hours after the flare at the boundary 
of the daylight side of the earth and at the same 
time on the geomagnetic night side. 


I. Review of observations of solar flares, 
SIDs, and solar radio noise on 6—8 July 
1958 

Table I presents data concerning solar flares 
of importance 2 or higher that occurred on 
7 July (see also Fig. 1). 

It seems likely that the observations No. 2—4 
concern one and the same flare. 

Most of the solar-flare correlated geophysical 
effects observed during the 7th of July and the 
days immediately following might therefore 
be attributed to one or both of the two flares, 
near the center of the sun’s disk, which started 
at 0032 and 0040 respectively. 

On 6th of July the High Altitude Observ- 
atory in Boulder, Col. (USA), gives infor- 
mation of only three flares, all of importance 1. 
For 8 July there was one flare of importance 3 
(outburst 0619 U.T.), one of importance 2, 
and 11 of importance 1 — to 1+. 

In spite of the fact that the sun was almost 
quite undisturbed on 6 July and the important 
flares of 7 July have been reported to have 


started not earlier than 0026 U.T. two observa- 
tories, Okinawa and Hiraisho in Japan, re- 
ported very strong gradual short wave fade-out 
(importance 3+ and 2+ respectively) with 
the starting time 0000 U.T. The duration was 
given as > 99 min and 140 min respectively. 

Hiraisho also reported the beginning of a 
sudden enhancement of atmospherics only a 
few minutes after 0000 U.T. 

Two sudden short wave fade-outs of im- 
portance 2 were observed in Canberra, one 
starting at 0030 U.T. and with a duration of 40 
minutes and the second beginning o110 U.T. 
and with a duration of 55 minutes. Starting 
times given are less certain than two minutes. 

Radio noise bursts, which might have been 
connected with the large solar flares, were 
reported at 0025 and 0027 U.T. fromHollandia 
(New Guinea) and Mitaka (Japan) respectively 
and at 0038 U.T. from Mitaka. The Mitaka 
observation at 0027 indicated a very strong 
but very short flare on 201 Mc/s. From about 
0103 to approximately one hour later, very 
strong radio noise emission on 201, 3,000 and 
9,500 Mc/s was reported from the same ob- 
servatory. At Sydney intense radio emissions 
of spectral types II and III were observed from 
0027 U.T. and some tens of minutes on. Strong 
receptions with starting times 0028 and 0040 
U.T. were also obtained on the single frequen- 
cies 600 and 1,420 Mc/s at that observatory. 

Thus, both optical flares and solar radio 
noise bursts were observed from the centrally 
located flares. The ultra violet radiation effects 
i.e. the short wave fade-outs were also noted. 

Figures 1 A, B, and C graph the solar flares, 
SIDs and major radio noise outbursts for the 
whole period 7—12 July. 


2. Methods of Measurements 


At Kiruna Geophysical Observatory the 
geomagnetic field components X, Y and Z are 


Table ı 


Time of the | Time of first 


Heliographic | Length of 


PI : er 
No. er beg. of the observ. er coordin. of | duration Quality of 
flare (U.T.) (U.T.) the flare (min.) observation 
I Mitaka ..... 0032 0032 2 26°N, 10°E = 45 Ghosts eh 
0026 and 0032 
2 Mitaka ..... 0040 0040 3+ 28°N, 07° W 
a Honolulu ... 0040 0040 3+ 25°N, 07° W 94 Excellent 
4 |Sidney ..... 0045 3 25°N, 10° W > 90 | Very poor 
5 Mitakom 0417 0417 2 ara N Era = 8 
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7 JULY 1958 
A FLARES 
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Figs 1. 
of the six parameters, 


18 24 06 12 


IGN (24 06; 1127187247067 12 71187729 


MET. 


Collocation of information on solar flares, SIDs, and major radio noise bursts and of records 
geomagnetic horizontal component, the virtual height of F2, 27.6 Mc/s cosmic 


noise power, oblique auroral reflection on 87.9 Mc/s, reception of 11.805 Mc/s transmissions from 
Alaska, and the meson component of cosmic radiation. 


recorded by means of three standard LA Cour 
variometers and a LA Cour quick-run magneto- 
graph. At Abisko magnetic observatory,! sit- 
uated 88 km northwest of Kiruna, the same 
types of instruments are used but they are 
arranged for measurements of D, H and Z.? 

The panoramic ionosounder was constructed 


1Some H-records from Abisko have been used as 
illustrations below, due to better quality of the records. 
The difference between the geomagnetic variations ir 
Kiruna and in Abisko is very small. 

? The magnetic recording work in Kiruna and Abisko 
is supervised by Dr. N. Amsott’s group at the Royal 
Board of Shipping and Navigation. 
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by Linpquist (1949).* The equipment sweeps 
over the frequency range 0.7—17 Mc/s in 
30 seconds. The maximum pulse output power 
is 16 kW. On a routine basis it gives three 
sweeps every quarter of an hour but during 
special IGY world intervals it sweeps contin- 
uously (two sweeps per minute). 

Since the early part of 1958 a riometer 
designed by C. G. LrrTLe at the Geophysical 
Institute, College, Alaska (vide LITTLE and 


3 It is operated at the Kiruna Geophysical Observatory 
under the supervision of Dr. R. LINDQUIST by the Re- 
search Institute for National Defense. 
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Fig. 2. Standard record of geomagnetic horizontal component at Abisko 6—9 July 1958. 


LEINBACH, 1958), is being operated at Kiruna 
on a cooperative basis between Geophysical 
Institute, College, and Kiruna Geophysical 
Observatory. The riometer is a self-balancing 
cosmic-noise recorder giving information about 
the vertical total ionospheric absorption for 
the operating frequency (27.6 Mc/s). Its Yagi 
antenna is pointed at the zenith. 

Since February 1958, aurorally propagated 
VHF signals from the Swedish F.M-broadcast 
transmitter in Ostersund (87.9 Mc/s) situated 
600 km SSW of Kiruna, and the Finnish F. M.- 
station in Kemi (92.8 Mc/s), 300 km to- 
wards SE have been recorded by means of a 
90 Mc/s rhombic antenna pointing towards 
magnetic north, VHF receivers and pen re- 
corders. 

Two IGY backscatter-sounder transmitters 
with frequencies 11.805 Mc/s and 17.900 Mc/s 
operated by the Geophysical Institute in Col- 
lege, Alaska, have been monitored in Kiruna 
since the early spring of 1958 for investigations 
of transpolar communication with special ref- 
erence to the effects of aurora. The trans- 
mitters are pulsed with a peak pulse output 
power of 4 kW. The output of the receivers 
modulate the intensity of oscilloscopes, the 
sweeps of which are synchronized with the 
pulse repetition frequency of the transmitters 
and the intensity of the light spots are recorded 
on continuously moving film. 

The meson component of the cosmic ra- 
diation is measured on a routine basis by 
means of two international cubic G.M.-tube 
telescopes of the type recommended by CSAGI 


and by another big meson telescope, built by 
Dr. A. E. SANDSTROM’s group at the University 
of Uppsala.! 


3. Observations in Kiruna during the first 
hour after the start of the flares 


Fig. 1 shows what happened with the 
previously mentioned six interconnected pa- 
rameters measured at Kiruna during the days 
7th to 12th of July. The geomagnetic record 
was very quiet on 6 July (Fig. 2) and the first 
hours of 7 July. The equipment sensitivity was 
such that a crochet in the H-record of 5 y 
amplitude was observable. No crochet could 
be observed which is in accordance with a 
statement of McNisu (1937), that solar flare 
crochets are never observed further from the 
subsolar point than 70°. In this case the angle 
was 87°. 

According to RAWER (1953), Mirra (1952) 
and other standard works a SID is observed 
over the entire sunlit side of the earth. No signif- 
icant absorption that might be associated with 
the solar flare was, however, observed in Kiruna 
(sun’s elevation 3.4°; cf. Figs. 1 F, 3 and 4). 
The SIDs observed in Japan and Australia were 
very strong, so the result seems to indicate 
that SIDs do not usually reach within 4—s 
degrees from the boundary between the day- 
and night-sides of the earth. During the first 
hour following the flares the sun reached about 
6—7° above horizon, so the change in the 
elevation was small. Such observational condi- 


! The cosmic radiation registration at the Kiruna 
Observatory is supervised by Dr. A. E. SANDsTROM. 
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Fig. 3 a—f. Ionograms a—c show the situation before, 5 minutes after and 20 minutes after the 3 + flare out- 
break at or40 MET (no significant absorption visible on b and c); d was obtained almost two hours after the 


flare (Es, but no other important change); 


e and f show the last two returns that were found before the long 


blackout period began. 


RIOMETER RECORD (27,6 Mc/s) 


0600 


0500 0400 0300 


OI00 0000 
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0200 


Fig. 4. Cosmic noise power measured by means of the Riometer during the hours before and after the 3 + flare 
outburst at 0140 MET on 7 July. 


tions are relatively unusual and can of course be 
obtained in the northern hemisphere only at 
points to the north of the polar circle during 
the weeks around the summer solstice within a 
few hours from local midnight. 

The propagation path of the transmissions 


0100 


from College, Alaska, to Kiruna was entirely 
over the sunlit side of the earth at the time 
of the flare outburst. It is therefore interesting 
that unusually good reception was obtained 
from 0035 to 0300 MET on 17.900 Mc/s 
(Fig. 5) and good reception was obtained also 


200 
0300 


D 5° im 


FL 


I" 


u 


N 


Fig. 5. Record showing reception in Kiruna of transmission on 17.900 Mc/s from College, Alaska. The re- 
ceived signal modulates the intensity of the beam of an oscilloscope, the sweep of which is approximately syn- 


chronized with the pulse repetition frequency of the transmitter. 


The periodicity of the record is due to the 


rotation of the transmitter antenna. 
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on 11.805 Mc/s. This indicates that the observed 
SID had little effect even for angular distances 
from the sub-solar point below 60 degrees 
(at the solar flare outbreak the sun’s elevation 
was about 40 degrees at College). 

Two equipments with which strong solar 
radio noise bursts can be recorded were being 
operated at the time of the flares; namely, the 
riometer and the recording receiver for recep- 
tion of VHF aurorally propagated signals. The 
direction to the sun was almost perpendicular 
to the plane of the vertical riometer three 
element Yagi antenna. The angle between 
the sun and the lobe direction of the 90 Mc/s 
thombic (pointing towards magnetic north) 
used for reception of auroral reflections was 
about 55° at 0130 MET. Thus, the sun was far 
out of the beam for both antennas. In spite 
of this, high intensity signals, not obtained 
during four months of operation, were re- 
ceived on both equipments.! 

On the riometer record of Fig. 4 can b: 
seen a strong intensity increase on 27.6 Mc/s 
starting at 0204 MET and lasting for about 12 
minutes. It coincided with an intensive burst 
of long duration observed at the Mitaka 
Observatory on the frequencies 201, 3,000 and 
9,500 Mc/s. The signal was quite different in 
fading character than interference, indicating 
that the band width of the increasing radiation 
was at least of the order of the sweep range of 
the Riometer (100 kc/s).1 

The broadcast transmitter at 89.7 Mc/s 
used for recording VHF-auroral reflections 
halted operations after local midnight. The 
transmitter off times are indicated in Fig. 1 G. 
A strong signal of short duration was recorded 
during the off period of 7 July (starting at 
0128 MET). Since the transmitter was off, 
the signal is assumed to have been of solar 
origin, which is supported by the coincidence 
with the strong burst on 201 Mc/s observed 
in Mitaka starting at 0219 MET. 

The counting rate of the meson telescope 
did not show any marked variation for the 
hours immediately following the flare out- 
bursts. 

Thus the first hour after the start of the 
flares was characterized by the absence of a 
SID, the absence of a magnetic crochet, and a 


1 Similar records have been obtained later in the sum- 
mer of 1958 in connection with strong solar radio noise 
Outbursts. 
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lack of change in the meson counting rate. 
However, the data show solar noise bursts 
on 27.6 Mc/s and 87.9 Mc/s receivers and 
enhanced reception of signal propagation 
across the pole from Alaska to Kiruna. 


4. The behaviour of the observed parameters 
during the next few days after the flare 


a. The Geomagnetic Field 


As mentioned above the geomagnetic field 
was unusually quiet on 6 July and the excep- 
tionally quiet condition continued during the 
first 8—10 hours of 7 July. From the forenoon 
of 7 July to 0845 MET of 8 July the H-curve 
was of the most common type obtained at 
Kiruna, being only somewhat disturbed with 
the fine structure of the curve having somewhat 
higher amplitudes than on a slightly disturbed 
day. The deflections did not surpass 100 gam- 
mas except during the normal night time bay 
(CE IFR) 

At 0845 MET on 8 July the magnetogram 
for the H-component showed a sudden com- 
mencement followed by a rapidly varying 
trace with peak amplitudes of the order of 
only 100 gammas. At 1115 MET H began to 
increase very rapidly. A positive peak deflection 
of 1,200 gammas above normal level occurred 
at 1250 MET. The smoothed curve crossed 
the Zero-variation line again between 1700 and 
1800 MET and decreased to 1,000 gammas 
below zero at 1810 MET. Except for some 
sporadic small positive deflections H stayed 
negative then until noon of 9 July. During 
the following days similar curves with steadil 
decreasing amplitudes appeared (Fig. 1 D). 

The large scale appearance with positive 
deflections during the days and negative de- 
flections during the nights is the normal one 
at Kiruna. The common interpretation of this 
data with a rotating earth in a declining 
auroral zone current system with the current 
direction from west to east on the dayside 
and in the opposite direction on the night 


side is well adaptable. 


b. The Ionospheric Absorption 


The ionosounder records showed a slow 
continuous increase in absorption during the 
morning hours of 7 July giving a decreasing 
intensity of the echoes and an increasing fmin 
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(compare Fig. 3a—d with 3 e and 3 f). After 
0600 MET no extraordinary ray echo appeared 
until 10 July. Fig. 3f shows the last record 
containing any return before the very long 
blackout (73 hours) started. The next echo 
appeared first at o910 MET on 10 July.! Seve- 
ral days followed with ionospheric reflections 
appearing and disappearing producing an 
oscillating blackout condition. A weak sporadic 
E layer was present before the flares started and 
could be found up to 0415 MET with only 
about half an hour of interruption. At the 
disappearance of Es the fminF was approxi- 
mately 3 Mc/s. 

The last receptions from College, Alaska, 
were obtained at 0230 MET on 7 July on 
11.805 Mc/s, and did not reappear again until 
11 July at 0730 MET (Fig. 1 H). BBC London 
interfered usually very strongly on that frequen- 
cy during the whole day and night except for 
the forenoon hours. The interference could be 
seen up to 0630 MET on 7 July and returned 
on 9 July at 1700 MET, which was 38.5 hours 
before the Alaska transmissions returned. On 
17.900 Mc/s the last reception from College 
was recorded at 0330 MET on 7 July. It appears 
that normal propagation conditions returned 
quite a good deal earlier where the midpath 
was south of the auroral zone than across the 
pole. 

On the riometer record an increasing ab- 
sorption began within two hours after the 
start of the flares (Figs. 4 and 6). The decrease 
in signal strength of incoming cosmic noise 
(measured in db) was nearly linear through 
the morning hours and had at 0700 MET on 7 
July reached - 2.7 db below the average value 
for that hour for the previous week. The 
absorption increased all through the day and 
reached a maximum a few hours after midnight 
between the 7th and 8th of July (Fig. 6). For 
about 10 hours the cosmic noise intensity 
level differed from the zero line (obtained 
when calibrating with the noise diode shorted) 
by only a very small amount. During those 10 
hours the cosmic noise level never was higher 
than - 17 db. Due to the uncertainty in the 
readings of the record it can only be said with 


1 As can be seen from Fig. 1 there were some equip- 
ment troubles during the nights, but there are reason to 
believe that no echoes occurred during those periods 
(there was total blackout in Troms6; TONSBERG, 1958, 
personal communication). 
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Fig. 6. Absorption in db as function of time, obtained 
from the Riometer record. 


any certainty that the lower limit of signal 
strength was lower than - 20 db. 

As far as we are aware absorption as great 
as this level has not been reported carlier in 
the literature.! During the solar flare of 23 
February 1956, a maximal absorption of 1.8 
db during night time and 8 db during the day 
was reported from College, Alaska (LITTLE and 
LEINBACH, 1958), which was on the night side 
of the earth at the time of the flare outburst. 

The sudden commencement of the magnetic 
storm occurred at 0845 MET on 8 July but 
not until 1150 MET did the first high positive 
peak value in H appear (Fig. 2). During that 
transition period a large increase in cosmic 
noise signal strength from about - 17 db to 
a peak value of — 5.2 db was observed (Figs. 6 
and 7). After that period the absorption never 
reached its maximum value again. A very slow 
recovery began in the afternoon of 8 July 


1 An absorption of 27 db was measured by the group 
at Sodankylä, Finland, after the flare of 7 July by means 
of a radio astronomy interferometer. (Personal communi- 
cation to one of the authors by Prof. DIEMINGER.) 
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Fig. 7. Riometer record showing the strong decrease of absorption during the first few hours following the sudden 
commencement at 0845 MET. 
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RIOMETER RECORD (27,6 Mc/s ) 
KIRUNA 9-10 July 1958 


Fig. 8. Examples of correlation between changes of the geomagnetic H-component and of the ionospheric ab- 
sorption. 


bringing the signal strength back to its normal 
value in about 3 days. 

The smoothness of the absorption record 
during the whole period 7—9 July is remark- 
able, and indicates a homogeneous absorption 
over a very large area. The contrast to absorp- 
tion caused by auroral forms is evident. 

During the recovery period two strong 
magnetic disturbances occurred, one positive 
and the second negative. On both occasions 
disturbances could also be seen on the riometer 
record (cf. Fig. 8). The positive magnetic dis- 
turbance, however, only caused a small smooth 
increase of absorption, while the disturbance 
crossing from a positive to a negative value of 
the horizontal magnetic component gave a 
sudden strong increase of absorption. This 


might indicate that different mechanisms are 
operating for production-of positive and nega- 
tive magnetic disturbances. The matter will be 
studied further with emphasis on the study of 
periods when the magnetic field variations 
change from positive to negative or vice versa. 

During the period of short term appearance 
and disappearance of blackouts on 10—12 July 
a fairly good correlation existed between 
blackouts on the ionograms and the absorption 
record on the riometer (see Fig. 1 and Fig. 8). 
The riometer records were characterized by 
sudden absorption dips while the ionospheric 
records showed no returns at all for periods 
ranging from several minutes to several hours. 

For the hours following the flare outbreaks 
and up to the blackout no deviations from the 
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normal behaviour could be seen on the f,F, 
curve. A comparison with the average curves 
for June and July of 1956 and 1957 showed a 
small increase from year to year of the critical 
frequency of the F, layer. The variations of 
the curve of midnight to os00 MET, 7 July 
1958, were similar to those of the same hours 
for June and July of 1956 and 1957. 

The above can be said for the h’F, curve 
too, up to 0500, at which time the F, layer 
disappeared completely. From 0500 on only the 
lowest height from which echoes were received 
could be obtained from the records. That 
height is of course also a function of the 
absorption of lower layers. 

Thus it should be emphasized that while 
the absorption effects were obvious during this 
period of 0230 to osoo from the decreased 
cosmic noise background signal, from the 
decreased intensity of the returns on the 
ionosond and from the increase of lowest 
usable frequency, both the height of the 
F, layer and the critical frequency of the same 
layer were unaffected. 

Up to 0908 MET on 10 July there was 
total blackout. During the following days 
blackouts of long and short duration were 
appearing and disappearing. The ionospheric 
returns would be present and then completel 
absorbed. If all—even very short blackouts 
were taken into account—no clear relation 
could be seen between f)F, before or after 
a blackout. However, in investigating only 
blackouts with a duration longer than two 
hours (five occurred between 10 July and 14 
July) it was found that in all cases f,F, was 
lower after the blackout than before. The 
average change was about 0.5 Mc/s. Again 
looking at only these longer duration black- 
outs there is a clear tendency after the end 
of the blackouts for a slow increase of f,F, 
at least during the first hour following the 
blackout. 

Concerning h’F,, it is more difficult to give 
any statement on the basis of the actual data, 
because it could fairly seldom be read from 
the records. The lowest height from which 
echoes were obtained, however, showed a very 
clear tendency to be a good deal higher after 
a long duration blackout than before (average 
change 91 km for the five blackout periods 
mentioned above). When short blackouts were 
included in the statistics the tendency was 
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found to be not quite as definite. Fig. 1 E, 
which gives the value of h’F, or, when not 
observable, the lowest echo height, illustrates 
this. It can also be seen that the echo height 
almost always decreased during the first hour 
after the end of a blackout. 

Thus, in analyzing the data from the longer 
lived but still relatively short duration blackouts, 
we find that the critical frequency of the F, 
layer immediately after the blackout was 
lower than before the blackout. The height of 
the F, layer was higher after the blackout than 
before indicating a possible expansion of the 
layer during the blackout period. During the 
next hour or so following the blackout the 
foF, would gradually increase and h’F, would 
gradually decrease. 


c. 87.9 Mc/s Oblique Auroral Reflections 


On Fig. 1G a few small amplitude bursts 
recorded with the auroral reflection receiver 
can be seen during the latter part of the 7th 
of July. These short duration signals are fairly 
common. Unlike the bursts mentioned pre- 
viously at 0128 MET on 7 July, which occurred 
while the transmitter was off, these signals 
were received when the Ostersund, Sweden, 
transmitter was on the air. The absorption had 
reached a very high value by this time but it 
was still many hours before the magnetic 
storm-causing beam arrived at the earth. 
These bursts received during the broadcasting 
period of 7 July were in general of much 
shorter duration than the normal oblique 
auroral reflections received before 7 July. 

Not until the afternoon of the 8th of July 
(at 1700 MET) did new echoes of the ordinary 
fading type appear. Some hours later (begin- 
ning at 1827) an echo type was recorded that 
had never been observed earlier in Kiruna 
during four months observation. The fading 
spectrum of the signal had an appreciable 
portion within the frequency range of the 
recorder and the fading amplitude was very 
large. The recorder was overloaded during 
the peaks of incoming signal (see Fig. 9a). 

At the time of the auroral reflection the 
absorption curve showed a smooth, slowly 
rising curve (Fig. ıF) with no indication of 
any rapid changes taking place in that part of 
the ionosphere which was in the beam of the 
vertical Riometer antenna. 
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Fig. 9. Records of oblique auroral reflections on 87.9 Mc/s. Transmissions from the Swedish FM-station Östersund 
were recorded in Kiruna by means of a rhombic antenna pointing to magnetic north. 


The magnetic X-component had just reach- 
ed the quiet day level again after a very 
large negative deflection. The Z-component 
had returned to a quiet level after negative 
deflections of great amplitude. The horizontal 
current parallel to the auroral zone was locali- 
zed to the south of Kiruna, if we interpret the 
magnetic deflections as caused by a horizontal 
auroral zone current system. 

During the largest magnetic deviations no 
auroral reflections could be seen, not even 
during the first few hours when the deflection 
in Z indicated that the current was flowing 
north of Kiruna, which is remarkable. As very 
strong reflections occurred a few hours later, 
when the absorption still was extremely strong, 
it is doubtful that the lack of reflections during 
the earlier period was due to absorption of 
the scattered radiation. It seems more likely 
that the aurora in spite of the magnetogram 
indications was to the south of Kiruna, which 
is supported by the information in the High 
Altitude Observatory’s weekly report that the 
aurora was observed down to 39 degrees lati- 
tude. 

Several suggestions might be put forward 
concerning the interpretation of the strong 
fading reception. As no effect could be seen 
on the cosmic noise incoming signal strength 
the aurora must have been very local and 


outside the beam of the vertical Riometer an- 
tenna. The corpuscular beam giving rise to 
the scattering aurora might be strongly modu- 
lated. However, the rapid rise and decay of 
signal strength are too fast to be caused by 
attachment, plus the fact that the decay time 
of the turbulence is supposed to be longer 
than a few seconds. A possible explanation of 
the strong fading amplitude might be that the 
number of scatterers giving rise to the signal 
was small. As the incoming signal strength was 
very high the turbulence and the electron den- 
sity must have been very high. 

What seems to be the most probable is, 
however, that the signal was not due to 
auroral propagation but came in via other 
modes of propagation from the south, the 
signal strength being so high that it could be 
recorded by means of the back lobe of the 
rhombic antenna. The latter assumption is 
supported by the following facts: F.M. trans- 
mitters in central Europe could be listened 
to at Sodankylä, Finland, during 8 July (Ka- 
TAJA, 1958, personal communication), and the 
fading of the record was fairly like that of 
records obtained when receiving the transmis- 
sions from the Finnish F.M.-Station in Kemi 
with the Yagi antenna pointing towards the 
transmitter. If this explanation were correct 
the condition of the lower ionosphere was very 
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Fig. 10. Intensity of the meson component of cosmic radiation, recorded by means of the vertical channels of a 
G.M.-tube telescope. 


unusual to the south of Kiruna rather than to 
north of it. This is in accord with the fact 
that blackout was first observed after the flare 
in Uppsala, Sweden (geomagnetic latitude 
58.4° N) between 1600 and 2300 MET, on 8 
July but not before that time interval. A 
final possibility is that signals came from both 
south and north. 


Further studies will be devoted to this type 
of propagation with special emphasis on finding 
out whether such signals are arriving from the 
south, or if the auroral propagation mechanism 
in special cases can differ fundamentally from 
the common type. 


On 9 July a long duration auroral reflection 
of common fading type was observed (Fig. 9b). 
During the days 10—12 July a slow decrease 
to ordinary reflection activity took place 
(Fig. 1G). 


d. The meson component of cosmic radiation 


The sudden commencement of the geomag- 
netic storm occurred at 0845 MET, 8 July, 
and a decline could be seen in the counting 
rates of the meson telescopes within the first 
hour thereafter. Fig. 10 diagrams this period, 
with the points representing the countrate 
integrated over one hour before the time in 
question. The decrease continued for about 
ten hours, reaching at 1900 MET a value 
some 5 % lower than during the morning 
hours. After that a slow regression began, 
lasting for several days (Fig. 11 and 10). The 
start of this FORBUSH type decrease, coincident 
with the start of the magnetic storm, seems 
to have been the normal behaviour at the time 
in question (SANDSTROM, 1958; personal com- 
munication). 
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Discussion 


The most interesting feature of the disturb- 
ances mentioned above is without doubt 
the strong ionization in the lower ionosphere 
on 7 July, showing up in the very strong 
absorption of HF radio waves, and the fact 
that no magnetic effects could be observed. 
Such an effect was observed the first time in 
connection with the large solar flare of 23 
February 1956 (according to SHAPLEY and RoB- 
ERTS, 1957) and several reports were published 
on it (BAILEY, 1957; LIED, 1957; SHAPLEY and 
ROBERTS, 1957, and others). 

At least two mechanisms have been proposed 
for ionization of the lowest layers of the 
ionosphere. Based on observations by VAN 
ALLEN and his co-workers (cf. MEREDITH et al., 
1955) CHAPMAN and LITTLE (1957) have postu- 
lated that bremsstrahlung is responsible for the 
ionization. That mechanism may be operating 
during auroral displays and geomagneticstorms 
of ordinary type but for flares of the kind 
discussed in this report it seems improbable 
that the necessary large electron flux should 
not cause any currents at all in the ionosphere. 

BAILEY (1957) has assumed that heavy ions 
(chiefly calcium ions) were causing the strong 
absorption effect of 23 Februari 1956. The 
main difficulty with his hypothesis seems to 
be that effects should be observable far below 
60 degrees geomagnetic latitude if calcium 
ions are the ionizing agent. 

Arguments can be given for the hypothesis 
that the absorption is caused chiefly by solar 
ions of mass numbers below that of calcium 
(cf. HULTQvIST 1959a) and it seems even to be 
possible to explain hitherto observed effects 
with only protons as ionizors. The flux neces- 
sary for the observed absorption is extremely 
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small, so the ions can move in STORMER or- 
bits. Such a hypothesis also makes understand- 
able the geographic distribution and the almost 
linear increase of absorption during the first 
8—10 hours (Huttevist, 1959b). 


It seems probable that the corpuscular radia- 
tion producing the effects discussed in this 
report is identical with that one deduced from 
radio observations of the sun by Wırp, Ros- 
ERTS and Murray (1954). 


The rapid decrease of absorption at the 
start of the magnetic storm may possibly be 
due to heating and expansion of the ionosphere, 
bringing the electrons towards higher levels 
where the absorption effect is smaller. From 
the fading characteristics of the record and the 
time displacement, the possible influence on 
the cosmic noise signal strength of the major 
radio noise burst at 0840—0850 MET on 8 
July can be ruled out. 


The proposed interpretation is in accordance 
with the well established fact that the F,-layer 
in polar regions very often is destroyed at the 
beginning of a magnetic storm (RAWER, 1953), 
and also with the dynamo theoretical explana- 
tion of sudden commencement. ECKERSLEY’S 
(1942) theory for the destruction of the F,- 
ayer, mentioned, seems not to be adaptable 
(cf. HULTQVIST, 1959a). 


Summary 


There are many interesting facets to the 
details of this particular polar blackout. The 
solar flares occurring around local midnight at 
Kiruna produced SIDs in Pacific areas but not 
at sunlit Kiruna. Within two hours after the 
commencement of the solar flares, the ioniza- 
tion of the lower atmosphere had started. 
However, during these first few hours of 
increasing absorption, both the height and the 
critical frequency of the F,-layer remained 
unaffected. The absorption kept increasing 
reaching a maximum about 24 hours after the 
flare. 

The effect of the absorption was to produce a 
blackout of vertical incidence ionospheric re- 
turns and oblique trans-polar reflections for 
approximately three days. The first ionospheri- 
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cally reflected signals to return were trans- 
missions (BBC-London) where the propagation 
midpoint area was situated to the south of 
Kiruna. 16 hours later vertical returns and still 
22 hours later transpolar signals were received. 

During flare observations the magnetic re- 
cords were undisturbed. 31 hours after the 
flare a sudden commencement magnetic storm 
started. Simultaneously the cosmic ray equip- 
ment registered a clear decline in meson count- 
ing rates. During the magnetic storm oblique 
signals from distant transmitters at 87.9 Mc/s 
showed signs of very fast and very deep fading, 
similar to sputter on scatter circuits. 

The ionospheric and magnetic storm period 
which followed the long blackout period 
showed many of the well-known effects; i.e., 
appearance and disappearance of ionospheric 
reflections, sudden short dips in the cosmic 
noise intensity, and large intensity variations 
of the horizontal component of the earth’s 
magnetic field. 
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On the Interpretation of Ionization in the Lower Ionosphere 
Occurring on Both Day and Night Side of the Earth 
within a Few Hours after some Solar Flares 


By BENGT HULTQVIST, Kiruna Geophysical Observatory, Sweden 


(Manuscript received November 28, 1958) 


Abstract 


Some strong solar flares cause very strong absorption in high geomagnetic latitudes on 
both day and night sides of the earth. The absorption effect is measurable a few hours after the 
solar flare outbreak. It is shown that the effect can be interpreted as being caused by a high- 
energy-ion beam of very small density, emitted from the sun and moving in Störmer orbits. 


Large ionospheric absorption effect following 
a solar flare outburst within a few hours 
without giving rise to any geomagnetic or 
auroral effects was reported for the first time 
from many ionospheric observatories north 
of 60° geomagnetic latitude after the big solar 
flare of 23 February 1956 (cf. SHapLEY and 
ROBERTS, 1957). 

Since the end of March 1958 a riometer 
(Relative Ionospheric Opacity Meter) designed 
by C. G. Lirtre at Geophysical Institute, 
College, Alaska (concerning construction and 
operation vide LITTLE and LEINBACH, 1958) has 
been in operation at Kiruna Geophysical Ob- 
servatory. That apparatus records the cosmic 
noise vertical flux on 27.6 Mc/s. Relative 
values of the total absorption in the whole 
ionosphere can be obtained if the absorption 
is supposed to be zero on days with very high 
recorded cosmic noise influx (that the absorp- 
tion is very small on such days will be shown 
below). 

During the period of operation three solar 
flares followed by heavy ionospheric absorption 
(peak value > ıs db) have been observed, 


namely on 7 July, 16 August and 26 August 
1958. The effects of the former flare on several 
upper atmosphere parameters as observed at 
Kiruna Geophysical Observatory have been 
reported (Hurrevist, AARONS and ORTNER, 
1959). 

In spite of the small number of observed 
flare effects a preliminary report will be given 
here and an interpretation will be proposed. 
The latter might be considered as a working 
hypothesis until it has been confronted by 
more complete observational material. Due 
to the rarity of the phenomena and the declin- 
ing solar cycle we will perhaps have to wait 
for a long time for such statistical data. 


Description of the absorption observations 


The outbursts of the solar flares of impor- 
tance 2+ and 3 + respectively of 7 July, which 
with all probability were responsible for the 
absorption, occurred at 0030 and 0040 U.T., 
fairly near the center of the sun’s disk. On 
7 July the sun was circumpolar in Kiruna 
(geographical coordinates 67.8°N, 20.4° E; 
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geomagnetic coordinates 65.3° N, 115.5° E). 
At the time of the flare outbursts the elevation 
angle was 3.4 degrees. Kiruna was, however, 
at the same time on the night side geomagneti- 


cally. 
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Figs. 1—3. Collocation of information on solar flares, 
SIDs and major radio noise bursts and records of the five 
parameters geomagnetic horizontal component, blackout 
periods for the ionosounder, 27.6 Mc/s cosmic noise 
power, oblique auroral reflections on 87.9 and 92.8 Mc/s 
respectively, and the meson component of cosmic radia- 
tion for the three periods 7—9 July, 16—18 August and 
26—28 August 1958 respectively. The small arrows above 
the cosmic noise curve indicate the time of sunset and 
sunrise in Kiruna at the earth’s surface. 


On 16th of August an importance 3 flare 
was observed to start at 0430 U.T., the location 
on sun being 13°S, 49° W, at which time 
Kiruna was on the sunlit hemisphere. Finally 
an importance 3 + flare outbreak followed 
by a very strong SID (importance 3 +) was 
observed at 0005 U.T. on 26 August 1958, at 
which time it was night in Kiruna. 

About two hours after the flares on 7 July 
and 26 August and approximately four hours 
after that of 16 August a slow decrease in 
incoming cosmic noise signal began to be 
visible on the riometer record. 

The variation of absorption with time for 
the three periods being considered in this 
report can be seen from Figs. 1—6. Figs. 1—3 
gives a schematic collocation of data, collected 
at Kiruna, for several parameters interconnected 
with the ionospheric absorption. The values 
of Figs. 4—6 have been reduced from the origi- 
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Figs. 4—6. Absorption as a function of time for the three 
periods 7—9 July, 16—18 August and 26—28 August 
1958 respectively. The maximum limits of errors due to 
reading of the record are given for a few points to indicate 
the order of magnitude. The uncertainty in the relative 
position of neighbouring points is far less. 


nal Easterline-Angus records for every hour. 
During the night between 7 and 8 July and in 
the afternoons of 16 August and 26 August 
the absorption reached its maximum, amount- 
ing to more than 20 db (the signal strength 
was so near zero that only a lower limit can 
be stated with any certainty), 16 db and more 
than 20 db respectively.! As far as we are 
aware these peak values are the highest ab- 
sorption values reported for vertical radiation 
around 30 Mc/s. For instance, LITTLE and 
LEINBACH (1958) found a day time absorption 
peak value of about 8 db after the solar flare 
of 23 February 1956, and a corresponding 
night time value of 1.8 db for 30 Mc/s cosmic 
noise. CHAPMAN and LITTLE (1957) report 
that the peak value of the absorption in College 
(geomagn. lat. 64.5° N) “may reach or even 
exceed 10 db by day and 5 db by night”. 
As Kiruna had midnight sun during the 
first absorption period no decrease of absorp- 
tion during the night was observed. During 


! The group at Sodankylä, Finland, measured 27 db 
absorption on 7 July (DIEMINGER 1958, private communi- 
cation). 
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the August periods, however, very marked 
night time decreases (down to about 2 db) 
can be seen from Figs. 2, 3, 5, 6 and 7. That 
effect can be attributed to solar control of an 
electron detachment process in the ionosphere 
(cf. e.g. CHAPMAN and LITTLE, 1957; BAILEY, 
1957; HULTQVIST, 1959). 

During the first hours of the magnetic storms 
on the days after the flares the absorption 
decreased to 5.2, 1.3 and 3.5 db respectively 
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for the three periods in chronological order, 
and then increased again. 

About 24 hours after the flares a slow 
decrease (the diurnal variation not taken into 
account) of absorption began. The decrease 
was, however, often interrupted by periods of 
low increase (Figs. 1 and 6). The increases and 
the short time changes of the curve can in most 
cases be attributed to special events as aurora, 
magnetic storm effects etc. Slow increases 
(such as those during the night between 8 and 
9 July) must be attributed to an ionizing agent 
or possibly drift motion. 

In spite of what was said above about ir- 
regularities in the curve during the strong 
absorption period the smoothness of the rio- 
meter curve is much more remarkable. A 
comparison between the first and last days of 
Fig. I gives an idea of that. This indicates that 
the ionization is homogeneous over a very 
large area, in contrast to the ionization caused 
by aurora. 

The most remarkable feature of the whole 
phenomena is, however, the absense of currents 
(except such causing quiet day magnetic varia- 
tions) in the ionosphere (Figs. 1—3) before 
the sudden commencement (which was ob- 
served also after the flare on 23 February 1956; 
cf. SHAPLEY and ROBERTS, 1957). 

Only observatories above 60 degrees geo- 
magnetic latitude observed absorption effects 
on the night side after the flare of 23 February 
1956 (SHAPLEY and ROBERTS, 1957). For the 
three flares discussed here similar restriction 
to geomagnetic latitudes greater than 60° has 
been found. In all cases the Swedish ionospheric 
observatory at Lycksele (geomagn. lat. 62.5°N) 
observed long duration blackouts but the 
duration was shorter than in Kiruna. At Upp- 
sala Ionospheric Observatory (geomagn. lat. 
58.4° N) no blackout was observed during 
any of the periods of interest. 


Review of hypothesis for the night-time 
ionization of the D-layer 


Based on observations made by VAN ALLEN 
and his coworkers (vide MEREDITH et al., 1955) 
CHapMman and LITTLE (1957) have proposed as 
explanation of the low altitude ionization in 
connection with aurora, that bremsstrahlung 
photons emitted by some thousandth of in- 
coming electrons in the usual auroral heights 
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Fig. 7. Increase of the cosmic noise signal strength during night hours. Arrows indicate sunrise and sunset at the earth’s surface, 
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is the ionizing agent for the lowest ionosphere. 
That mechanism may very well be operating 
in connection with aurora, but there seems to 
be some important difficulties in adopting it 
to the solar effect discussed here, namely that 
the large absorption observed makes it nec- 
essary to assume a fairly large number density 
of electrons in the incoming beam, which 
seems difficult to combine with the total lack 
of geomagnetic and auroral effects. 

That the ionization in the D-layer is caused 
by fairly high energy ions coming from the 
large solar flares, in cases of the kind being 
discussed in this report, has been proposed by 
BAILEY (1957) and EHMERT and REVELLIO (1957). 
In his comprehensive article Baıey has 
discussed the problem in detail. There seem, 
however, to be some difficulties left, the most 
important one being that his adopted ion 
masses and energies would, with all probabil- 
ity, make the effect detectable far below 60° N 
geomagnetic latitude. He has also neglected 
the ionization of the primary ions in his 
calculations, which ionization seems to be very 
important (cf. below). 


Analysis of the phenomena and outline of a 
mechanism for ionization of the D-layer 


Velocity of particles: 


The possibility of determining the velocity 
of the particles giving rise to the absorption is 
fairly poor. From the riometer record we 
have to determine when a decay in signal 
strength that is not of daily variation type 
starts. The daily variation sometimes changes 
significantly within a week or so, as the 
cosmic noise temperature in the antenna beam, 
moving over the sky, changes, which increases 
the difficulty. 

The technique adopted is to construct an 
average daily curve from undisturbed days 
within one or two weeks before the flare out- 
burst of interest. The uncertainty in the time 
interval between the flare outburst and the 
beginning of decline of cosmic noise power 
reaching the antenna is estimated to be not 
more than an hour. The determination of 
such a “travel time” for the solar particles 
does not rule out the possibility of existence 
of much faster particles in the beam, which, 
however, must be too few to cause measurable 
absorption. 


The travel times determined correspond to 
specific rectilinear velocity values. If the trajec- 
tories of the particles were Stérmer orbits in a 
dipole field the measured travel times would 
correspond to slightly higher velocity values 
for the rigidities in question (see below). 
The magnetic field in interplanetary space is, 
however, certainly far from a dipole field. 
The sun’s magnetic field may vary from time 
to time due to distortions caused by solar 
flares (ALFVEN, 1956). This and possible irre- 
gular distribution of interplanetary matter 
having a magnetic field makes it probable 
that the interplanetary magnetic field is irregu- 
lar. That means that the trajectories of particles 
of a low density beam (cf. below) may be 
much longer than the rectilinear distance. This 
is the fact for instance with the interplanetary 
field proposed by ALFVÉN (1956) (cf. Brock, 
1958 a). As far as the author is aware it is, 
however, impossible to give any quantitative 
relations based on observational data. 


For the event of 7 July a travel time of about 
two hours could be extracted from the record. 
On 16 and 26 August the corresponding values 
were four hours and two hours respectively. 


After the solar flare that occurred at 0331 
U.T. on 23 February 1956 effects in the lower 
ionosphere on the earth’s night side were ob- 
served within 20 minutes after the flare out- 
burst (vide e.g. BAILEY, 1957; EHMERT and 
REVELLIO, 1957). 


Conclusions from effects of the earth’s magnetic 
field on the particle orbits : 


The fact that in all cases discussed above the 
effects were restricted to higher geomagnetic 
latitudes than approximately 60 degrees indi- 
cates that if particles originating in the sun 
are assumed to cause the ionization these cor- 
puscles must be charged. 


If the number density of the beam is supposed 
to be very small (arguments for that will be 
given below) then the particle orbits are Stér- 
mer orbits. From the velocity (unknown; 
only “travel time’’ can be measured) and the 
upper geomagnetic colatitude limit (©) reached 
by the particles, a lower limit of the heaviest 
mass number present in the beam can be 
obtained from the relation 


sin © = 0.01202 (Ho)!!4 (1) 
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Table 1. Mass number of particles reaching the latitude limit 
a a ee A ee 2 
Time interval between the flare outburst and observed 


Lower geomagnetic 
latitude limit 


start of absorption effect 


20 min. 2 hours | 4 hours 
60° | 2 (H—He) 13.1 (C—N) | 26.3 (Mg—Al) 
62° | 1.6 (H—He) 10.3 (Be—B) | 20.7 (Ne—Na) 


(where Ho is the rigidity of the particles; STÖR- 
MER, 1955, p. 296) with sufficient accuracy for 
the rigidity values of interest. Reasons why 
only a lower limit of the heaviest mass number 
present in the beam can be obtained for a 
given velocity value will be discussed below. 

All solar flare effects observed in the summer 
of 1958 had about the same geographical ex- 
tent, as mentioned above. The lower geomag- 
netic latitude limit was somewhere between 
62.5° (Lycksele) and 58.4° (Uppsala). For two 
different limit latitudes, namely 62° and 60° 
the mass number of the ions reaching the 
latitude, rectilinear trajectories being presumed, 
is given in Table 1 for the solar flares being 
discussed (time intervals between flare out- 
bursts and observed start of absorption effect 2, 
4 and 2 hours respectively) and also for a “travel 
time” of only 20 minutes, which should 
correspond approximately to the flare effect of 
23 February 1956. 

As can be seen from the table the mass num- 
ber which must, for postulated rectilinear 
orbits, be attributed to those corpuscles reaching 
the lowest latitudes differs very much between 
the four cases of flares, and in all cases it is 
lower than 40 (for calcium), which Batrey 
(1957) has generally claimed to be the ionizing 
agent. Thus calcium ions would reach lower 
latitudes than 60 degrees even if they had the 
lower limit velocity, corresponding to recti- 
linear orbits. 


Depth of penetration of the particles : 


Batey (1957) has worked out a very valuable 
diagram (copied in Fig. 8) that gives a good 
survey of the relations between velocity, mass 
number, lower geomagnetic latitude limit and 
penetration depths. The penetration depths for 
the heavier ions, which are theoretical, might, 
according to BAILEY (1957), be fairly much 
overestimated due to the neglect of stripping 
effects. The diagram might be used for a dis- 
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cussion of the effect, keeping in mind that 
numerical statements are fairly unsure. 

For the travel time range of interest the 
penetration depth (D) can be written 


D=mv- F(v) (2) 


(m is the particle mass and v the velocity) 
where F(v) is a function only of v and not of 
m for the particles. 


Discussion of the reason why strong absorption 
effects are not observable below about 60° geo- 
magnetic latitude : 


Let us for the moment assume that the mass 
spectrum of the beam emitted from the sun is 
fairly broad. The lower geomagnetic latitude 
limit is a function only of the momentum for 
equally charged particles (cf. (1)). If the pene- 
tration depth were also a pure momentum 
function (that means F(v) in (2) was a function 
of mv instead of v, which is equivalent to the 
“equal penetration depth lines” in Fig. 8 being 
parallel to the latitude lines) strong absorption 
could be observed in those latitudes where 
the corresponding penetration depth was 
suitable, and only in those latitudes, irrespective 
of the velocity of the beam. 

That the particles should not penetrate too 
far into the atmosphere is due to the fact that 
for very high electron collision frequencies (y) 
the absorption coefficient for HF radio waves is 
inversely proportional to v. This is discussed 
below. 

If we postulate a broad mass spectrum of 
the beam, the primary condition that must be 
fullfilled is that the penetration depth must be 
suitable. Suppose the limits of the suitable 
range are 60 and 85 km. We see that latitudes 
lower than 60° geomagnetically are reached 
only for very restricted values of velocity and 
mass number. For very fast beams (several 
tenths of the light velocity) no strong absorp- 
tion would be observed, because the penetra- 
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Fig. 8. Diagram showing the connections between rectilinear travel time sun to earth, 
velocity, kinetic energy, atmospheric depths of vertical penetration, lower geomagnetic 
latitude limit for singly charged particles, and kind of charged particle (After BAILEY, 1957). 


tion depth would be too great for the whole 
mass spectrum of the beam. However, the 
velocity spectrum of the beam may perhaps 
be fairly broad and there will also be some 
interaction between the beam particles and 
the atmospheric layers far above the height 
where the particles are stopped. Even parti- 
cles, which do not reach the D-layer but the 
lower E-layer, may produce measurable ab- 
sorption. We should therefore not expect 
very distinctive limits. 

Against the dominating role of heavy ions 
in the beam is the fact that the cosmic abun- 
dance of such ions is very small (less than one 
per mille), 91 % of the atoms being protons 
and 9% helium atoms (ALLEN, 1955). 

There are two difficulties for the hypothesis 
that protons are responsible for the absorption 


in and to the south of the auroral zone. One is 
that a high velocity of 0.5 c is required for them 
to reach even 65 degrees geomagnetic latitude. 
Such a velocity corresponds to a rectilinear 
travel time of approximately 20 minutes. In 
order to reach the earth not earlier than two 
and four hours after the start the length of 
the proton trajectories must be 6 and 12 times 
the rectilinear distance respectively. For 60 
degrees geomagnetic latitude the corresponding 
figures are about 10 and 20 times.! The other 
difficulty is that protons with such velocities 
will penetrate the atmosphere to depths far 


1 Of course permanent ring currents around the earth, 
storage of magnetic clouds in the vicinity of the earth 
etc. may decrease the lower geomagnetic latitude limit 
but that will not be discussed here as no quantitative 
results are available as far as the author is aware. 
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below 40 km altitude, which means that only 
a small amount of their energy will give rise 
to absorption-producing electrons in suitable 
heights. The penetration depths for protons 
given in BAILEYS diagram (Fig. 8) presume 
vertical rectilinear orbits in the atmosphere. 
The curvature of the orbits and scattering in 
the latter part of the penetration path, where 
the velocity has decreased appreciably, and also 
an existing small angle with the vertical direc- 
tion will decrease the penetration depth a little. 
Such a correction of the proton curve of Fig. 8, 
however, will certainly not amount to more 
than some 10 km, which is not enough to 
make the protons stop in the lower D-layer. 

However, protons in large enough amounts 
would give the observed absorption and there 
would also be a gradual decrease in absorption 
with decreasing geomagnetic latitude, due to 
the increasing penetration depths of the pro- 
tons, as has been observed. The needed increase 
in number density compared to an equivalent 
flux stopping just below the D-layer can be 
estimated from the “stopping cross section” 
(SEGRE, 1953). In order to produce the same 
ionization per cm in the upper layers there 
must be some 100 times more protons of 
200 MeV than of 1 MeV and about 10 times 
as many as of 15 MeV energy. For constant 
proton flux (number density of the beam in- 
versely proportional to the velocity) the absorp- 
tion effect would increase rapidly with in- 
creasing latitude. Postulating protons, the lower 
geomagnetic latitude limit would thus be due 
to the fact that even if there are protons of 
energy enough for reaching a lower latitude, 
interaction with the atmosphere would occur 
too deep within. For all flare effects discussed 
in this report the lower latitude limit has been 
the same which means that the upper limit of 
the proton energy spectrum must have been 
a few hundred MeV. It seems probable but is 
not necessary according to the interpretation 
given that flare effects of the type described 
here are more common inside the auroral 
zone. A detailed study of the effect at high 
latitudes may give information that can throw 
some light on the problem if the ionizing 
particles are protons or heavier ions. 

For helium ions all problems discussed 
above for protons will be reduced considerably. 
It seems probable that helium ions are respon- 
sible for a significant fraction of the absorption 
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effect, because of the more favourable geo- 
graphic distribution and penetration depth and 
their relatively high abundance in the solar 
surface layers. This assumption is supported 
by the discovery of helium lines in the auroral 
spectrum (GALPERIN, 1958). 

The fact that the absorption in connection 
with the high energy flare of 23 February 1958 
was much lower than the absorption effects 
observed during the summer of 1958 fits into 
the picture. 


Estimation of necessary number density of the beam: 
The absorption is given by 


0 


where A is the absorption in db; 
E, is the incoming field strength with- 
out attenuation; 
E is the attenuated field strength; 
k is the absorption coefficient, 


(3) 


the integral being taken over the path of the 
radiation. 

For Kiruna’s geomagnetic latitude and 27.6 
Mc/s vertical radiation the propagation is qua- 
silongitudinal. 

For that case k can be written 


p ß 
"a Gp À (4) 


where p is the angular frequency of the radia- 
tion; 


x = — p*/po; 
po =47Ne?/m; 
B =pip; 
YL =ppılp; 
pz = Hıe/me; 


Hr =the component of the earth’s magnetic 


field along the propagation path; 
v =collision frequency for the electrons; 
N =electron number density; 
e and m refer to the electron and 
c is the velocity of light. 
— is for the ordinary ray and 
+ is for the extraordinary one. 


For the case of interest the values of the absorp- 
tion coefficient are for the ordinary ray 
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Table 2. Normal electron density (N cm-?), electron collision frequency (v sec-1), and apparent 
recombination coefficient (y cm? sec!) as functions of height 
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3.20 + 1016 + y? 


(5) 


vN cm1 


0 = 


and for the extraordinary ray 


0.0529 
2.60: 1016 + y2 


(6) 


yN cm! 


e 


Table 2 contains values of N for the normal 
case, y and y (the apparent recombination 
coefficient). The values of N for the normal 
case are those given by WAyNICK (1957). 
The apparent recombination coefficient has 
been ed from Mirra and JONES (1954), 
whose values have arbitrarily been reduced 
somewhat to make them better adaptable to 
both day and night conditions (cf. Waynick, 
1957). The collision frequency values are taken 
from NICOLET (1953). 

With the values of Table 2 the normal 
absorption is found to be 0.5 db. To this should 
be added the absorption in the F-layers, which 
in any case is less than the given value. 

The term »? in (5) and (6) can be neglected 
down to 65 km. For that altitude the error 
introduced by neglecting this term will be 
only 4%. For heights where the values of » 
are greater than 10° the absorption decays 
rapidly with decreasing height. From extra- 
polation of Nicoret’s (1953) curves the »- 
value can be estimated to be attributable to a 
height between so and 60 km. 

If we assume that the ionization causing 
absorption is distributed in a parabolic layer 
with N=o for h=60 km and the peak value 
at a height of 85 km (which coincides with the 
normal peak height, vide Tabie 2)! a numerical 
integration gives, for the peak electron number 


1 For other height distributions the corresponding 
electron number density will be different. As the real 
height of the absorbing layer is unknown, it has been 
presumed that the peak height has the normal D-layer 
value. 


density corresponding to an absorption of 23 
db, the value 8 - 10 cm-3. The number den- 
sity distribution with height is given in Fig. 9 
together with the normal case of Table 2. 
The E-layer has arbitrarily been supposed to 
be unaffected. 

Taking into account the extraordinary ray 
reduces the value by some 7 %, and the effect 
of the broad beam of the antenna is supposed 
to make the number too large by less than 10%. 


ELECTRON DENSITY .CM® 


10 


ne 1 
60 70 80 90 100 110 
ALTITUDE, KM 


Fig. 9. Distribution of electron density (cm-®) with 

altitude (1) in a parabolic layer with a peak height of 

85 km and zero density at 60 km causing a nondeviative 

absorption of 23 db for 27.6 Mc/s vertical radiation (solid 

line) and (2) in an average D layer (according to Way- 
NICK, 1957). 
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Table 3. Electron productions giving rise to the two electron distributions of Fig. 9 as functions 


of height 
Height 6 
lege 5 70 75 80 85 90 95 100 105 
electr./cm? sec 
flare case 4.224702 Trot | Leto 102 17:7, 10% 5.22 102 13:0 - 102 7.8. - 102 77 23 
normal case | Salon Ir.8 . 10-216.4 Ê 1074/3. . 10-2] 18 | 5 | 0.45 [3.0 . 10-214.8-- 1072? 


Using the apparent recombination coefficient 
of Table 2 we find the values presented in 
Table 3 for the electron production per cm* 
and second, q, for the given altitudes. 

The values of Table 3 have been obtained 
from the relation q=yN?, which gives the 
“effective” electron production. This is equal 
to the total electron production only if the 
electron number density, N, is much greater 
than the negative ion number density, N, 
and if dN/dt is much less than q/(1 + À) and 
yN?(A= N-/N), that means in the stationary 
case. During the rapid increase of the ab- 
sorption the electron production was 3.5 times 
as high as the stationary value. An average 
value of N, weighted with the electron col- 
lision frequency (v), is found to correspond 
to a height of 68.5 km for the electron distri- 
bution of Fig. 9. With the average number 
of N, 1.4 10% cm #, the v-value corresponding 
to the same height (2.5 - 107 sec~1), the average 
value for the ion-ion recombination coefficient 
(3: 10-4 cm? sec!) obtained from the rio- 
meter records (cf. HULTQVIST, 1959) and a 
À-value of so (reasonable for the 68 km level; 
vide HULTQVIST, 1959) a stationary state value 
for the electron production (q) of 15 elec- 
trons/cm sec was obtained. During the rapid 
increase of absorption the q-value was found 
to be 52 electrons/cm? sec. That means that 
in a column with one cm? area 2.3: 108 
electrons per second are produced in the flare 
case. 

For mass numbers and velocity values of 
particles causing absorption down towards 60 
degrees geomagnetic latitude under the con- 
ditions postulated above the energy varies 
between roughly 10 MeV and 300 MeV. If 
the energy for production of one electron-ion 
pair is taken as 40 eV (which is certainly not 
too low) we can see that one ion produces 
between 3 : 10° and 5 : 108 electrons. In case 
of the strongest absorption effects an influx of 
50—1,000 ions/cm? sec is required to produce 
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the absorption. The corresponding veloci- 
ties vary between 0.03 c and 0.4 c, approxi- 
mately, as deduced from the rectilinear travel 
time. If only protons and helium ions are 
supposed to be effective the velocity figures, 
deduced from the lower geomagnetic latitude 
limit, vary between approximately 0.2ando.7¢. 

The effective number density of the beam 
moving with the velocities mentioned must 
thus be within 5 - 10-® to 1- 10-8 cm, if heavy 
ions are assumed and between 2 : 10° and 
2. 10°” cm? if protons and helium ions are the 
ionizing agents. 

In the discussion here only the ions effective 
in the atmosphere have been considered 
hitherto. It has already been pointed out that 
protons will penetrate to great depths in the 
atmosphere. As has also been mentioned before 
this fact means that we have to increase the 
density by a factor not greater than 100 and 
probably of the order of 10, to take into 
account the energy release to layers giving 
negligible radio wave absorption. 

The ions of the beam are focused towards 
the poles of the earth which means an increase 
in density compared to the conditions far 
from the earth. On the other hand the major- 
ity of the ion trajectories in the vicinity of 
the earth do not reach the atmosphere, which 
counteracts the focusing effect. Due to lack of 
quantitative information, the counteraction is 
supposed. to be exact. 

If on the other hand moderately heavy ions 
is the effective agent we must assume that they 
are accompanied by protons which reach the 
ionosphere much nearer the poles. If not very 
selective emission mechanisms are operating 
the total beam density must be higher than the 
density of the effective ions by at least a factor 
1,000. 

The densities mentioned above are extremely 
small values compared to those which are 
thought to exist in beams causing aurora and 
magnetic storms. CHAPMAN (1950) and FER- 


342 


RARO have estimated the minimum density 
of a neutral ionized gas sufficient to produce a 
moderate magnetic storm at 100 electron-ion 
pairs and Martyn (1951) has given a corre- 
sponding value of 20. CHAPMAN (1950) has 
stated that a density of far below one ion per 
cm® is necessary for application of STÖRMER’S 
theory. According to the electric field theory 
the number density of a beam giving rise to a 
magnetic storm is roughly o.1 electron-ion 
pair per cm? (Brock, 1958). Whichever of 
the named values of number density in mag- 
netic-storm-producing beams is adopted, we 
see that the values for the beam discussed 
above are less by several powers of ten. 

Thus the existence of the observed intense 
ionization without any currents seems reason- 
able. 

The linear increase of absorption during the 
first hours is a natural feature of the interpre- 
tation given (cf. HULTQVIST, 1959), as is the 
observed maximum in effect during the night 
between 7 and 8 July when the sun was over 
the horizon, and at the same time the higher 
corpuscular flux on the geomagnetic night 
side was acting. 


The decrease of absorption during the first few 
hours after the sudden commencement of the geo- 
magnetic storm : 


In the three observed cases of the effect 
under discussion a very marked decrease of ab- 
sorption occurred after the sudden commence- 
ment (Figs. 1—6). This indicates that the in- 
teractions of the two beams (the first causing 
absorption only, the second causing magnetic 
disturbances) with the ionosphere are different. 

It is a well known experimental fact that the 
F2-layer very often in the polar regions is 
quite destroyed at the start of a magnetic 
storm (RAWER, 1953). That the F-layers were 
strongly affected by the dense beam is sup- 
ported by the fact that the F-echoes returned 
at lower total absorption (as measured by the 
riometer) than they disappeared at. It was 
also observed after the flare of 7 July 1958 
that the F-echo heights were usually about too 
km higher at the echo’s return after blackouts 
than when the blackouts began (cf. Hurr- 
QVIST, AARONS and ORTNER, 1959). 

The most simple explanation for the de- 
crease of the absorption seems to be that the 
ionization is removed from the height region, 
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where its absorption effect is maximal, either 
by a vertical electric field or by thermal 
effects (heating and expansion) caused by the 
magnetic storm producing beam. 

The increase of absorption following the 
transient decrease may be attributed to the 
geomagnetically active beam, and if protons 
are ionizers in the high energy beam, perhaps 
also to a gradual expansion of the lower 
electron-rich regions of the ionosphere up into 
levels of higher absorption. 


Discussion 


The interpretation presented above gives a 
possible explanation of the similarity in effect 
irrespective of the difference in beam velocity, 
the geographical extension, the total absence of 
any magnetic effects, and the gradual increase 
of the absorption during the first hours after 
the flare outburst. 

As given above there may be two slightly 
differing ways of interpreting the effect. One 
postulates a broad mass spectrum of the beam 
and velocities of the ions differing only slightly 
from those corresponding to rectilinear travel 
time. The other is based on the presumption 
that the ionization in the lowest ionosphere is 
caused by protons and helium ions of high 
velocity, in which case either irregular trajec- 
tories in the interplanetary space, of lengths 
about one order of magnitude larger than the 
rectilinear distance sun to earth, which trajec- 
tories also must differ from time to time, must 
be assumed, or a profound influence of ring 
currents and magnetic clouds near the earth 
on the ion orbits must exist. None of these 
alternatives seems to be highly improbable. 

It seems probable that the corpuscular radia- 
tion producing the effects discussed in this re- 
port is identical with that one deduced from 
radio observations of the sun by Wırp, Ro- 
BERTS and Murray (1954). 

The basic feature of the absorption effect, 
namely the lack of magnetic deviations in spite 
of the dense ionization, seems to be under- 
standable with any of the alternatives of the 
hypothesis that solar ions moving in Störmer 
orbits are the agent. 

The experimental bases—three absorption 
periods—is of course very meagre and more 
thorough investigation of a greater number 
is highly desirable. It will be carried out as 
soon as alarger observational material has been 
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obtained. Due to the rarity of the phenomena 
this might take several years. 

Measurement of the absorption by means 
of riometers at as many points as possible 
differing in geomagnetic latitude by only a 
few degrees from 55° towards higher values 
would probably give information that would 
tell us, which of the proposed alternatives of 


mechanism (if any) is operating. 
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Evaporation of Boric Acid from Sea Water’ 
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Abstract 


Previous investigators have shown that the boron-chlorinity ratios of rain waters are many 
times greater than the boron-chlorinity ratio of sea water. The presence of boron in the at- 
mosphere has been attributed to sea spray, volcanic activity, accumulation in dust, evapora- 
tion from plants, and industrial pollution. In this paper data are presented to demonstrate that 
boric acid in sea water has a vapor pressure at ordinary temperatures of the sea and, when sea 
water evaporates, boric acid occurs in the condensate of the water vapor. It is postulated that, 
while some of the boron in the atmosphere can be attributed to the sources mentioned above, 
most of the boric acid results from evaporation from the sea. 


Introduction 


SUGAWARA (1948), in his studies of the 
composition of rain water, noted the presence 
of boric acid and stressed the fact that the boron- 
chlorinity ratio was hundreds of times greater 
than the boron-chlorinity ratio of sea water. 
Muro (1952) collected four different samples 
of water during a continuous three-day rain. 
His analyses showed that the chlorides decrease 
more rapidly than the boric acid content as 
the rain continued, and that the boron-chloride 
ratio was from 50 to 270 times greater than 
the boron-chlorinity ratio of sea water. Muro 
(1953) examined the boron and chloride 
content of snow and obtained results similar 
to those on rain water. 

Several theories have been advanced to ex- 
plain the presence of boric acid in the atmos- 
phere. TRUNINGER (1944) attributes the boron 
to the accumulation of fine dust; KALLE (1945) 
to volcanic activity and cyclic salts supplied 
by sea water; LANDERGREN (1945) to vol- 
canism; and Muro (1956) toindustrial pollution 
of the atmosphere and to evaporation of boron 


* Contribution No. 229 of the University of Washing- 
ton, Department of Oceanography. 


compounds from plants. In support of this 
latter theory Muto showed that boron com- 
pounds are given off by some plants and that 
the rain waters over heavy vegetated areas 
have higher boron content than the rains 
and snows of high altitudes. 

The authors of this paper postulate that 
much of the boron content of the atmosphere, 
and thus of rain and snow, results largely from 
the evaporation of boric acid from the sea. 
Experimental evidence is presented to sub- 
stantiate this hypothesis. 

SCHAFFGOTSCH (1859) demonstrated that some 
boric acid is removed with the water vapor 
when an aqueous solution of the acid is evapor- 
ated. TCHIJEWSKI (1884) stated that the loss of 
boric acid by evaporation of aqueous solutions 
is not proportional to the concentration of the 
boric acid in solution. LESCOEUR (1886, 1890) 
made vapor pressure studies of boric acid and 
its aqueous solutions. At 20° C boric acid had a 
vapor pressure of 2 mm and at 43.5° C, 5 mm. 
HEHNER (1891) stated that orthoboric acid 
would vaporize from its solution at 100° C, 
but WATSON (1893) maintained that it was the 
metaboric acid that vaporized. KoNINGH (1897) 
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Fig. 1. Schematic cross section of the evaporators. 


A. Polyethylene carboy of so liter capacity 
B. Solution of NaOH 

C. Glass wool and NaOH pellets 

D. Cooling condensers 

E. Warming condensers. 


concluded that the evaporation of orthoboric 
acid noted by Schaffgotsch applied only to 
concentrated solutions and that dilute solutions 
could be evaporated without appreciable loss. 
This small loss can now be easily determined 
and is most pertinent to the present problem. 

In the well-known mannitol method for the 
determination of boric acid, the carbonates are 
first destroyed by acidification and the resulting 
carbon dioxide is expelled by boiling. In this 
process the analyst always employs a reflux 
condenser to prevent loss of boric acid. It has 
been recognized for many years that boric acid 
in aqueous solutions has a definite vapor 
pressure depending upon the temperature and 
concentration. As far as the present authors can 
ascertain, this knowledge has not been con- 
sidered in investigations of the evaporation of 
sea water nor offered as an explanation for 
the presence of boric acid in rain and snow at 
atmospheric temperatures. RANKAMA and Sa- 
HAMA (1949) cite referenccs to some geochem- 
ical studies which state that the boric acid is 
volatile at 100° C. 


Experimental 

Methods of Analysis. The improved mannitol 
method of Gast and THompson (1958) was 
used in the determination of boric acid in sea 
water. The colorimetric method of PHILIPSON 
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(1953), because of its sensitivity for deter- 
mining traces of boric acid, was employed in 
the analysis of distillates and condensates, a 
Beckman DU spectrophotometer being used 
for the measurement of optical densities. 
Distillation of Sea Water. 500 ml of sea water 
with a chlorinity of 16.22 % and containing 
3.753 mg/l of boron as boric acid was distilled, 
using a Hopkins head, until 100 ml of distillate 
was obtained. Distillates from samples of the 
same water, which had been acidified before 
distillation, were likewise collected and ana- 
lyzed. Averaged results were as follows: 


Sea water contained 

3.7537 mg/l boric acid-boron 
Distillate of untreated water, 

0.1126 mg/l boric acid-boron 
Distillate of acidified water, 

0.1408 mg/l boric acid-boron 


None of the distillates showed the slightest 
trace of chlorides. 

Evaporation of Sea Water. Both the carly 
literature cited above and the fact that some 
boric acid is removed from solution by distilla- 
tion indicate that boric acid in solution has a 
definite vapor pressure at atmospheric tempera- 
tures. When sea water evaporates the water 
vapor should also include small quantities of 
boric acid. To demonstrate this the apparatus 
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shown in Fig. 1 was assembled. Two carboys 
of polyethylene were used (A). One carboy 
was filled with so liters of sea water having 
a chlorinity of 16 %, and the other with the 
same quantity of distilled water. Air was 
gently blown over the surfaces of the waters. 
No spray or bubbles formed. The air was 
first cleaned by passing through a concentrated 
solution of sodium hydroxide (B), passed over 
pellets of sodium hydroxide (C), to remove 
much of the water vapor, and then warmed 
in condensers (E). A common air supply, 
used simultaneously in both carboys, eliminated 
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the possibility of variable contamination, 
should such a condition exist. The dried air 
entered the carboys at a temperature of 32° C, 
plus or minus 3° C. The waters in the carboy 
were kept at a temperature of 25° C. The air, 
after passing over the water surfaces, entered 
cooling condensers (D), to condense much of 
the water that had evaporated. Portions of 
too ml of condensates were collected and 
analyzed. In one series of experiments the tem- 
perature of the cooling condensers was 12° C, 
and in another series 0°C. The results are 
summarized as follows: 


Air Water Condenser Boric Acid-Boron 
Temp. Temp. Temp. in Condensate 
HE SE Ae ug/L 
Distillede Water nae sure ware secre dem arene sek exe Bee 25 12° 5.0 
Sea a Wiaterie 1e eee emeie eee cet 322 255 127 60.0 
Seas Waterers ee mec once eine olenere 322 25° où 65.0 


The condensate from the distilled water 
yielded the equivalent of 5 ug of boric acid- 
boron per liter. This result is considered as a 
blank determination. 

Evaporation of Dilute Solutions of Boric Acid: 
Dilute solutions of so liters of boric acid were 


prepared by dissolving boric acid in distilled 


Solution, 
Solution, 
Solution, 


1.083 
2.166 


3.249 


water. The concentrations of these solutions 
were 1.083 mg, 2.166 mg, and 3.249 mg 
boric acid-boron per liter, respectively. Dry 
air was passed over these solutions and the 
condensates were collected and analyzed as 
described above with the following results: 


Air Water Condenser Boric Acid-Boron 
Temp. Temp. Temp. in Condensate 

°C ° IC ug/L 

325 25 722 19.0 

322 255 12° 42.0 

324 254 T2> 60.0 

32% 25° o° 66.0 


Solution, 3.249 


Graphical portrayal of the above results in-di 
cates that the concentration of boric acid evapo- 
rating is proportional to the boric acid in 
solution and confirms the observations of 
SCHAFFGOTSCH (1859). 


Discussion of the results 


In both series of experiments with sea water 
and the experiment with dilute solutions of 
boric acid, the boron content of the condensates 
at 0°C showed slightly higher concentrations 
than at 12° C, indicating a vapor pressure of 
boric acid at these lower temperatures. Thus, 
the boron content of the condensates is actually 
less than the boron that had evaporated from 
the sea water and the solutions. 

The vapor pressure of boric acid is a function 
of the temperature and the salinity (boric 
acid concentration) of the sea water. The rate 


at which botic acid may evaporate or be 
‘steamed distilled’ from the sea surface is 
directly proportional to the rate of evapora- 
tion. The concentration of boron in the rains 
or snows resulting from this evaporation is a 
function of the temperature of condensation 
and the temperature at time of precipitation. 

In tropical regions, during periods of little 
rainfall, the evaporation of boric acid into the 
atmosphere would tend to give a boron- 
chlorinity ratio of the surface waters slightly 
lower than normal. In areas where water and 
air temperatures were considerably lower, the 
amount of boric acid entering the atmosphere 
would be less and thus lessen the effect on the 
boron-chlorinity ratios of the waters. These 
ratios might even be slightly higher than 
normal if cold waters were subjected to the 
influence of air and moisture coming from 
tropical or warmer regions. 
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In a recent paper, the authors (1958) cited 
data by the United States Geological Survey 
showing the boron content of several rivers 
of the Pacific Northwest. The waters of these 
rivers on the western slope of the Cascade 
mountains represent the drainage of rains 
falling in the region. These rains had their 
origin over the Pacific Ocean. The boron 
content of the river waters are comparable 
with the boron found in the condensates ob- 
tained from sea water in the above experi- 
ments. 

Morpy (1957), in his report of the 3rd 
Annual Conference on Atmospheric Chemistry, 
cites the data of several investigators who 
reported on the analyses of the boron content 
of rain waters collected over northern Europe. 


347 


These results are much lower than those re- 
ported by any of the Japanese investigators; 
however, all the data tend to support the 
hypothesis that much of the boron in the 
atmosphere and in rains and snows is the result 
of evaporation from sea water. Most of the 
rain falling over Japan had its origin in warmer 
regions of the Pacific Ocean. The origins of 
the rains mentioned by Mordy are more com- 
plicated, having come from the evaporation 
of colder sea water or from the evaporation 
of fresh water on land areas. 
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On the Frequency Response of Some Different Models Describing 


the Transient Exchange of Matter between the Atmosphere 
and the Sea 
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Abstract 


A theoretical study is made of the response characteristics of two “box-models’’ and 
two continous models of the sea which could describe an air-sea exchange of matter. 
It is demonstrated quantitatively that the box-model can depict a transient exchange 
process correctly only when the time of internal mixing of the sea is small compared 
to a certain “transfer time”. Comparing the two continous models studied, one of which 
assumes a purely diffusive sea and the other a purely advective sea, one finds some inter- 
esting differences in the responce characteristics. 


In recent years there has been a great deal of 
interest in transient air-sea exchange processes 
involving different radioactive and nonradio- 
active elements. As one example we may men- 
tion the uptake in the sea of radioactive mate- 
rials produced by atomic bomb explosions. 
Another example is the extra addition to the 
sea of carbon dioxide that has been released in 
the atmosphere during the last half century 
by combustion of fossil fuels. 

In special, this last problem has been sub- 
jected to several theoretical investigations 
(CRAIG 1957, REVELLE-SUESS 1957, BOLIN- 
ERIKSSON 1959). In these investigations use 
has been made of so-called “box-models”, in 
which the atmosphere and the sea are considered 
as separate and internally well-mixed reserv- 
oirs. In some cases the sea has been split up 
in two reservoirs, one comprising the upp- 
ermost layer down to about 100 m depth, 
the “mixed layer”, and the other comprising 
the deep sea. 


The general structure of such a box-model 
is as follows. Each box i is characterized by 
the amount of material under consideration N; 
and by a number of transfer coefficients kj. 
The flux of material from the reservoir i to 
the reservoir j is assumed to be ki;N; (Fig. 1). 
One can then set up the continuity equation 
for the box i as follows: 


Fig. 1. The general box-model. 
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dN; fe % & L 

an oe iy ees wt ft) 
Gj #1) (iJ) 

St is the rate of increase of material, P; is the 
eventual production in the box and - yN; 
represents the decay, if the material is radioac- 
tive. Writing down such continuity equations 
for all the boxes i=1, 2,...... n gives usa 
closed system of equations for the deter- 
mination of the N;:s as a function of time. 

The transient behaviour of the box-model 
is correct under the assumptions 

(i) that the exchange between the reservoirs 

is effectuated by a steadily working 
transport mechanism, 

(ii) that each reservoir is internally well- 

mixed. 

The first assumption asserts that the model 
is linear. Increasing all concentrations to the 
double valuesshould simply double all the fluxes 
between the reservoirs. The second assumption 
asserts that these fluxes are proportional to 
the total amounts of matter in the boxes. 
While the first assumption on the whole 
seems to be justified the second assumption is 
more critical. It requires that the matter is 
uniformly mixed in the reservoir within a 
time tY short compared to the “transfer times” 
I 


=; characterizing the reservoir. I£ there is a 


ij 
production and radioactive decay of matter in 
the reservoir, 7? must also be small compared 


to a characteristic production time 7} and the 


decay-time t? ee It should be stressed 


that this is the requirement that the box-model 
depicts correctly the transient change of the 
reservoir (the deviations from the momentary 
equilibrium state). When studying the over- 
all transient change of a system it is not abso- 
lutely necessary to consider the transient 
changes of every reservoir. Those reservoirs 
which have mixing- and transfer-times very 
short compared to the period of interest can 
often be treated as equilibrium reservoirs, and 
even if they are not well-mixed in the above 
sense a boxmodel representation may again 
be justified. 

The purpose of the present note is to compare 
the behaviour of two different box-models and 
two continous models of the sea in a simple 
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air-sea exchange problem. In all the cases the 
atmosphere is treated as an equilibrium reser- 
voir, which seems to be justified in most cases 
of practical interest. For the sea we assume: 


(i) a one-layer box-model 
(ii) a two-layer box-model 
(iii) a diffusive sea 

(iv) an advective sea 


The investigation is restricted to the case of 
a non-radioactive material, it may become 
actual to study later also the case of radio- 
active matters. 

Contrasting the results obtained for the 
box-models and the continous models should 
yield some precise information of the defects 
of the box-model in such cases where the 
above-mentioned assumptions are not satisfied. 
The comparison of the results valid for the two 
continous models should also be of interest. 
The question whether matter in the sea is trans- 
ported mainly by advection or diffusion seems 
so far not to be settled, and it is possible that 
the application of theoretical results such as 
are derived here could cast some further light 
on this. 

In the comparison between the models use 
is being made of the so-called frequency res- 
ponse diagram, which has found much use 
in engineering sciences. It should certainly pro- 
ve useful in discussions on transient exchange 
and transport problems in geophysics, which 
are governed by linear equations. An advantage 
of the technique that should bepointed out 
is that it can also be applied to observed data. 
This should be of special value when de aling 
with such transport links as, say, the biologic 
transport, for which it seems very hard to set 
up a specific mathematical model. 


Case (i) 

The simple box-model representation of 
the sea is seen in Fig. 2 a. The governing equa- 
tion for the amount N, of our matter in the 
sea is 


dN. 

he k, N, = RNs (2) 
Considering N, as a “forcing function” and 
N, as a “response function” we like to solve 
N, 
——, where AN, and AN, are 
A N, a s 
the complex amplitudes in a harmonically 
varied state: 


for the ratio r= 
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ka No VA kN, 
Sea 
Fig. 2 a. Simple box-model for the air-sea exchange. 


Atmosphere 


Deep sea 


Fig. 2 b. Two-layer box-model for the air-sea 
exchange. 


Ne ze-h 


OZ 


Fig. 2 c. Diffusive sea. 


N,=N,+ AN,e 
(3) 


N,= N; + AN, ei?” 
One finds easily 


rw)= Gra 
AN, k, tio 


I 


= Nees 


(4) 


I+i— 
Ws 


J k 
with r(o)= =, 
k; 

tude ratio can be represented by a curve in 
the complex plane. In this case the curve is 


w;=k;. This complex ampli- 
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Fig. 2 d. Advective sea. 


the half-circle seen in’ Fig. 3. From this fre- 
quency response curve we can read the am- 
plitude and phase of AN, relative to AN, for 
any frequency. As example, for © =w, we see 
that the amplitude (magnitude) of AN, is 


2.0.71 of the value given by an equilibrium 
2 
theory (w=o), and that the phase difference 


between AN, and AN, is—~. The frequency 


response curve obviously contains all infor- 
mation about the model. Every forcing func- 
tion N,(t) can be split up in harmonic com- 
ponents, the N,-responses of which can be 
read from the response curve, and these har- 
monic responses can be summed up to give 
the complete time-response N,(t) of themodel. 


Case (ii) 

The two-layer box-model of the sea is 
presented in Fig. 2b. The governing equations 
are now 


dN„ E 
dt br ka Na a ka Na = (Rie A Rd) Nin 

dN, 

= — Rina Nm = kaNa (5) 
dt 
N, = Nin Er Na 


O O5 
oo | Men 


05 
-0,5 


Fig. 3. Frequency response curve for the simple box- 
model. The figures along the curve give the values 
œ 
Of —. 
Ws 
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ie) 0,5 1,0 


-0,5-- 
k 
| ee fe) 


kmd 
Fig. 4 a—c. Frequency response curves for the two- 
layer box-model. The figures along the curve give the 


@ k a . 
values of —. The case — = 0 corresponds to Fig. 3. 
Ws md 


k 
The dashed line in Fig. 4 c represents the case ee ES 


md 


k 
In all cases LS er a. 
kmd 50 


N, is again the total amount of our matter in 
the sea. In this case one finds, 


r(@) = N = 
I+ _ Rmaka _ joe 
» (Rind Ar ka)? os 
ent > Le oe .O bs Rigs Rind (2)’ 
er 7) Fo (ma + ka)? \o, 
(6) 


Since the deep sea can store a very large 
amount of matter compared to the mixe 
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layer (Na> Nm) we should have ky<knq and 
(6) simplifies to 


kma\.© Rina f © \? G2) 
+ (ine (>) 


The form of the response curve depends here 


r(w) = r(o) 


À : 
on the two parameters = and —-. As is seen 


Rind Rind 
from (6a), if we let +o for a fixed ratio Fas 
md 
we get again the solution of the simple box- 
model. This result should be anticipated, since 
we have here a case where the internal adjust- 
ment between the mixed layer and the deep 
sea becomes infinitely rapid. For finite values of 


Er 
kind 
ably from the result of the simple box-model. 
Assuming that the depth of the mixed layer 
is about 100 m and the total depth of the 


oceans about 5,000 m we would guess that 


k 
es: and using this value the frequency 


Rind 


response curves have been computed for some 
The 


the result may, however, differ consider- 


different values of the parameter =" 


Rnd 


result is seen in Fig. 4. 


Case (iii) 

In the case of a diffusive sea (Fig. 2c) we 
assume that the matter introduced at the top 
of the sea spreads vertically over a depth h by 
turbulent diffusion. The coefficient of diffusion 
is assumed to be constant. In the horizontal 
directions conditions should be uniform every- 
where. It is of interest to note that this diffu- 
sive model is approached by a multi-layer 
box-model where the division in horizontal 
layers is made more and more refined. 

The governing equation for the diffusion 
process is 
an, 


on, 
Ge K 722 (7) 


n, is the mass concentration of the matter, K 
the coefficient of diffusion, and z a vertical 
coordinate. The boundary conditions are 
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’ y on, 
kana = kin, = Qk 
at z= 0 (the sea surface) (8) 
Au 
dz 


at z= —h (the sea bottom) 


where ki, ki are transfer coefficients valid at 
the sea surface and 9, is the density of water. 
(8) expresses the equality between the flux of 
mass through unit area of the boundary layer 
at the sea surface (including both the atmos- 
pheric and the oceanic boundary layers) and 
the turbulent mass flux in the sea just below 
the boundary layer (note that K is multiplied 
by @; to give an absolute mass transport). The 
transfer coefficients k} and k! are, of course, 
related to the coefficients k, and k, used earlier. 
In equilibrium the integrated flux of matter 
from the atmosphere into the sea is kin,A 
and the equally large return flux kin,A, 
where A is the sea surface area. These fluxes 
may, however, also be written k,N, and k,N,, 
respectively. With N,=n,M,, N,=n,M,, 
where M, and M, are the masses of the atmos- 
phere and the oceans, respectively, we find 
b kiA b kA 
OMR M: 
harmonic variation in time: 


As before we assume a 


Na = fa + Ane 
n,=n; + An,(z)e@ 


The solution to equation (7) which satisfies the 
boundary conditions (8) is then 


ki, An, 
An, ER Er has * 
, iw iw . iw 
ki cosh ve h+K ve sinh \ a h 
DER 
x cosh \/ = (z +h) (9) 


Introducing the actual values of AN, and AN, 
from the relations AN„=M,An, AN,= 


Oo 
= M,An,= My fan, dz and introducing also 


—h 
the transfer coefficients k,, k, given earlier we 
can finally write 

AN: 


mr ec) 


(10) 
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Fig. 5 a—c. Frequency response curves for the dif- 
fusive sea. The figures along the curve give the values 


© hz. Gite ale & 
of —. t* = —— is a characteristic time of internal 
Ws 2K 
er I ic 
mixing, and ts; = —. The case — = corresponds 
s Ts 
to Fig. 3. The dashed line in Fig. sc represents the 
tx 
case — = © 
Ts 


=> té bi) Ve Jo oh (ati) y= V+ 1 
Ts 


Ts ws ws ws 


6 il an 
Here r(0)=, w,=k,, Che t* is a new 
Ay 


k, 
> h2 
parameter with the value t* = SK It represents 
3 
a characteristic mixing time for the sea. If 
this time is very small compared to r,, we find 
again the response curve of the simple box- 
model, as expected. For other values of the 
| à 
T the 
parameter — we get deviations from the box- 
Ss 


model curve, as seen from Fig. 5. 
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Case (iv) 

In this case all internal transports in the sea 
are advective. The specific model is sketched 
in Fig. 2 d. On top we have a thin surface 
layer, which is assumed to be well mixed 
horizontally. From this surface layer water is 
advected into the deep sea, carrying with it 
our matter, and returning to the surface layer 
again after a certain circulation time. This 
time will in general vary from one particle to 
the next. Consider now a small “flux-tube”’ 
having a given circulation time r. If the mass 
of water in the tube is dM and the mass 
flux through it is dF we have 


dM = tdF (11) 
The “Aux-mean” circulation time is 
F, 


T = = ip tdF, where F, is the total mass flux. 
Oo 


From (11) we find 
T=— (12) 


This time represents a kind of “mixing time” 
for our system. 

Let now F(t) be the distribution function 
for the flux after t (F(t) gives the total mass 
flux for which the circulation time is smaller 
then r). The distribution of M after r is then 
also fixed, by equation (11). Consider the 
balance of matter in the entire surface layer. 
There is a net transport from the atmosphere 
A(kin;—kin), as in the diffusive model 
studied earlier, there is an advective transport 


Fn, into the deep sea, and there is a return 
F, 


transport from the deep sea Jn; (t—1) dF (rt). 


oO 
Neglecting the capacity of the thin surface 
layer we have 


F, 
A(kitta — kins) - Fen; + fn;(t — 1) dF (1) =0 
(13) 


As before, we write n,=n,+ Angel", ns = Ns + 
Anse”. Note that n, is the momentary con- 
centration in the surface layer. From An, and 
An, we can come over to AN, and AN, b 

help of the relations 4N,=M,An, AN,= 
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M, 
An, f e—ie dM (t) (from the relation N, = 
M; 
Ing (t —t)dM(t)). Expressing k/, and k/ in terms 


of k, andk, by help of previous formulas and 
replacing dM by tdF according to (11) we 


get now 
u Je 
AN, ” 
r(®) AN = an (14) 


Tm represents here the maximum circulation 
time of the sea. To come further we must 
make some specific assumption regarding the 
function F(t). The simplest reasonable assump- 
tion we can make is that it is linear, running 
from o to F, when + runs from 0 to Tn: 


F(t) 7 
Er (15) 
Ss m 
In this case 7, =2T. We can now write 
r(w) = 
œo du 
li? ler ERP EST 
TW 
=1(0)— Sm 2 5 
iw -2i-— T w 
=: =e “#) HAS (+ 1) (2) 
Ds sNUs Ws 
(16) 


where r(o) = w,=k;. 


Again, making our “mixing time” 7 small 
compared to 7, we will get the result of the 
simple box-model, while for larger values of 


T à h 
— the result will be quite different. Curves 
Ts 


£ T : 
for some different values of — are seen in 
Ss: 
Fig. 6. 


From the discussion of the above four models 
we may verify what was stated in the begin- 
ning, namely that the simple box-model is a 

ood approximation only when the “mixing 
time” is small compared to the “transfer 
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time’. When the mixing time becomes 
comparable to or larger than the transfer time, 
the result deviates from that of the simple box- 
model in the way shown in the diagrams. It is 
of interest to note that the largest deviations 


generally occur at high values of = since such 
high values may be found in cases of practical 
interest. As example, in the problem of the 


carbon dioxide uptake in the sea the period of 


‘ r\n 
interest r(-2) is of the order of a few 
© 


I 
decades while t,=— certainly is expected to 
9) 


Ss 
take on values of the order of several hundred 
years. 

It is also of interest to note the considerable 
differences in the response characteristics of the 
diffusive and the advective sea, when the 
“mixing time” takes on large values. In the 
diffusive model the phase angle will be rela- 
tively small and negative, reaching a minimum 


value of about RTE The amplitude ratio is 


finite and is more or less independent of the 
ratio of the mixing time and the transfer 
time, once this ratio is large enough. All 
changes will then take place in a thin surface 
layer and the actual capacity of the reservoir is 
irrelevant. In the advective model the phase 
angle decreases indefinitely for increasing 
frequencies. Making the ratio of the mixing 
time and the transfer time large the amplitude 
ratio tends to zero. 


Fig. 6 a—c. Frequency responce curves for the advec- 
tive sea. The figures along the curve give the values 


oO =. : : ; I 
of —. T is the mean circulation time, and rt; = —. 
Ws Ws 
T à 
The case — = o corresponds to Fig. 3. 
Ts 
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Permeation of Gases through Ice 


By EDVARD HEMMINGSEN, Institute of Zoophysiology, University of Oslo 


(Manuscript received September 9, 1958) 


Abstract 


The permeation rate of carbon dioxide and oxygen through ice has been determined at 
several temperatures. At —9.5° C the carbon dioxide permeation constant was found to be 
0.7 x 10711 cm? sec~? atm™?, which is about two million times less than in water. The per- 
meation increased hyperbolically towards 0°C. No permeation of oxygen through the ice 
could be detected. It is concluded that the gas penetrates the ice through intercrystalline 


brine films rather than through the ice crystals. 


Knowledge of the permeability of gases 
through ice is useful on several scores. It has 
been shown that many northern plants and 
animals survive the winter frozen solid, some 
of them with 90°% or more of the water 
frozen to ice, and yet slow but measurable 
respiration continues down to — 20° or lower. 
A great many aquatic animals and plants have 
also been found to survive within the ice of 
bottom frozen arctic lakes and ponds. Since, 
demonstrably, freezing of the organism does 
not stop the respiration, the question arises 
whether gas exchange can possibly take place 
through the ice which entombs the organism. 
A few determinations undertaken by ScHo- 
LANDER, FLAGG, Hock and IRVING (1953) 
showed that in general this cannot be so, for 
the ice is extremely gas tight, indeed as tight 
as is glass to helium (NORTON, 1957). 

This opened up the interesting possibility 
that gas bubbles enclosed in ice might under 
ideal conditions remain unchanged for very 
long times, perhaps millenia. The whitish 
colour characteristic of icebergs is caused by 
masses of gas bubbles dispersed throughout the 
ice. If these were formed by trapping of air 
when the snow was compacted to form the 
glacier, analysis of the bubbles might accord- 
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ingly give information of the composition of 
ancient atmosphere. Such a possibility became 
even more challenging when it was demon- 
strated that glaciers can be radio-carbon dated 
through the content of carbon dioxide trapped 
in their bubbles (COACHMAN, HEMMINGSEN, 
SCHOLANDER, ENNS and DE Vries, 1958). The 
validity of such determinations clearly depends 
upon how stable are the gases in the ice. One 
factor which could interfere is gas diffusion 
through the ice. 

In the present investigation gas penetration 
through ice has been studied for carbon dioxide 
and oxygen at various temperatures. The 
method used is a modification on macro- 
scale of the micro method employed by ScHo- 
LANDER et al. (1953). The macro scale has made 
possible greater accuracy in measurements of 
ice thickness and area and a more accurate gas 
analysis. 


Methods 


Principle. A thin ice sheet floats on mercury 
and separates carbon dioxide gas, or oxygen, 
above the ice sheet from a known volume of 
air deposited as a blister underneath the ice 
sheet. After a known period of time at constant 
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Fig. 1. Method for determination of the permeability of carbon dioxide and oxygen through ice. 
A. Preparation of ice sheets. 
B. Assembly for permeability measurements. 
C. Introduction of gas under the ice sheet. 


temperature the air is analyzed, and the amount 
of permeated gas is calculated in relation to ice 
thickness and blister area. 

Experimental procedure. An ice membrane of 
uniform thickness was prepared between glass 
plates, spaced by a ring of waxed paper (Fig. 1, 
A). To prevent the ice from sticking, the glass 
plates were lightly greased and aluminum foil 
was smoothly rolled onto them. Boiled distilled 
water was flowed, while still hot, into the spacer 
ring between the plates, and left to freeze. When 
ready, the assembly was taken apart, and the 
foil was carefully stripped from the ice sheet. 
Only clear ice sheets containing an insignificant 
amount of air bubbles were used. The thickness 
of the ice was measured with a dial micrometer. 

The ice window was floated on a pool of 
mercury, and a known volume of carbon 
dioxide-free air was delivered with a pipette 
under the ice, where it formed a shallow gas 
blister (Fig. 1, C). An inverted Petri dish, 
resting with its rim in mercury, formed an 
air-tight chamber over the ice, through which 
a continuous flow of carbon dioxide or oxygen 


could be maintained (Fig. 1, B). To minimize 
thinning of the ice by evaporation cracked ice 
was placed around the ice window. With this 
precaution the average thinning was reduced to 
s—10% of the total thickness. If the thinning 
was allowed to proceed much further, there 
was risk of local break through, especially at 
the warmer temperatures. The assembly rested 
in a thermostated box. 

After a known period of time the gas 
blister was withdrawn and analyzed for carbon 
dioxide or oxygen, in a % cc gas analyzer, 
with an accuracy of + 0.015 % (SCHOLANDER, 
1947). The increase in carbon dioxide concen- 
tration was usually no more than about 0.03 %, 
i.e. barely above the accuracy of the gas anal- 
ysis. Higher concentrations would have been 
desirable, but would require longer exposure 
with consequent danger of break through. 
The blister area remained clearly visible on 
the ice sheet, and was traced on transparent 
graph paper for measurement. The ice thick- 
ness was measured anew and averaged with 
the first measurement. 
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Determinations of the permeation were made 
at —0.5°, — 1.5°, — 5.0° and -9.5°C. At each 
temperature three to four experiments were 
performed, with a new ice membrane each 
time. 

For oxygen determinations the same proce- 
dure was followed, using outdoor air in the 
blister and flowing pure oxygen through the 
chamber. 

Calculations. The quantity of gas which pene- 
trated the ice was calculated from the change in 
composition of the gas blister, disregarding the 
change in volume, which was too slight to be 
measured. According to Fick’s law, the per- 
meation constant K is given by the following 
equation: 

q 


APP)" 


K 


where q is volume of permeating gas in cm$, 
À is area of gas-ice contact in cm?, f is time in 
seconds, d is thickness in cm, and (P, — Pg) is 
the tension difference in atmospheres. The unit 
is accordingly cm® cm! cm”? sec"! atm-! and 
the dimensions cm? sec-! atm-}. 


Results and Discussion 


The data pertaining to the carbon dioxide 
permeation at various temperatures are given 
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in Table I and Figure 2. It will be seen that 
the carbon dioxide permeation increases very 
slightly from - 9.5° to - 1.5°, but rises steeply 
above - 1.5°. 

Four determinations at - 5.0° and two at 
-0.5° gave no measurable penetration of 
oxygen through ice. The permeation constant 
or oxygen was therefore less than 0.6 x 10-4 
cm? sec-! atm-t, which is the limit of accuracy 
of the method. 

The permeation constant for carbon dioxide 
through water at 10° C is given as 1.7 x 10-8 
cm? sec"! atm? (HANDB. CHEM. PuHys.; 
BRUINS, 1929), i.e. carbon dioxide moves 
about one million times slower through ice 
at — 5° than through water. 

The value for the permeation constant of 
carbon dioxide through ice at - 20°C was 
determined by ScHOLANDER et al. (1953) as 
0.006 to 0.01 mm? mm! cm? hour-! atm-!, 
or transposed to the present units 1.7 to 2.8 x 
10-10, which is some ten times higher than was 
found in the present determinations at — 9.5°. 
This discrepancy is far beyond methodical 
uncertainties, and is most likely caused by 
impurities in the lake water used for the 
determinations (SCHOLANDER, personal com- 
munication). In agreement with the present 
determinations these authors found no pene- 
tration of oxygen through the ice. 


Table 1. Permeation constant of carbon dioxide at various temperatures. 


Permeated Permeation 
Temperature * Time Ice-thickness| Bubble area | Bubble vol. CO; constant x 1074 
(2€) (hours) (cm) (cm?) (cm?) (increase (cm? cm! cm? 
in) vol4.9% sec! atm-!) 
— 0.5 10.0 0.072 + 0.002| 15.4 40.1 2.10 + 0.05 | 0.05 + 0.01 14.4 + 2.7 
» 10.0 0.070 » 14.8 » 2.00 » 0.03 » 7.9 » 
» 10.0 0.073 » 14.3 » 2.00 » 0.03 » 8.5 » 
= 1.5 26.0 0.039 » 15.2 » 2.00 » 0.01 » 0.6 0.6 
» 41.3 0.038 » 15.1 » 2.00 » 0.07 » Ze AR OS 
» 27.0 0.046 » 14.9 » 2.00 » 0.03 » 198 20.0 
— 5:0 38.8 0.081 » 13.6 » 1.80 » 0.01 » 0.8 08 
» 26.8 0.075 » 7.8 » 1.00 » 0.03 » 3.0 » 
» 53.0 0.081 » 9.5 » 1.50 » 0.02 » 1.3 » 
— 9.5 25.3 0.072 » 14.9 » 2.00 » 0.01 » Eon eed 
» 35.5 0.045 » 14.5 » 2.00 » 0.02 » TON 10.0 
» 63.7 0.042* » 14.3 » 2.00 » 0.01 » 0.340.033 
» 72.8 0.042 » 14.5 » 2.00 » 0.02 » 0.4 » 
—20 a — — — — 20* — 
» — — = —- == 27 = 
» — a= — — = DRE 


*) SCHOLANDER, 
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Fig. 2. Permeation of carbon dioxide through ice in relation to temperature. Dimen- 

sions of the permeation constant are cm? sec! atm~!. The data are compared to 

two arbitrary drawn hyperbolas. The oxygen permeation is less than the value 
indicated by the stippled line. 


One may think of two possible pathways 
for the permeation of gas through ice, either 
through pores in a dry crystalline ice structure, 
or through a liquid phase which might exist in 
the ice. In the first case, according to Graham’s 
law, one would expect that the smaller oxygen 
molecules would move through the ice in 
greater amount than the carbon dioxide mole- 
cules, but if permeation takes place through 
liquid pathways carbon dioxide would pass 
through in greater amount because of its 20 
—30 times greater solubility. The latter was 
clearly the case, and this strongly suggests that 
we are dealing with permeation through a 
liquid phase. 

It has been established that impurities in 
the water, such as dissolved salts, concentrate 
as brine films on the surface of the ice crystals, 
when the water freezes (QUINCKE, 1905; Dor- 
SEY, 1940; RENAUD, 1949). The concentration 
of these films will very likely be proportional 
to the ice temperature below zero degrees 
Centigrade, and the volume of unfrozen brine 
will hence decrease hyperbolically with de- 
creasing temperature. Since the brine channels 
are limited in length by the ice thickness, it 
follows that their cross section area must also 
decrease hyperbolically with temperature, and 
consequently also the gas permeation. Increase 


of brine concentration would slow down the 
permeation and would tend to accentuate the 
effect. It will be seen from Figure 2 that the 
data indeed approximate a hyperbolic tem- 
perature dependence of the permeation con- 
stant, and this, plus the fact that it is carbon 
dioxide rather than oxygen that penetrates 
the ice, furnishes independent support to the 
concept that permeation takes place through 
intercrystalline brine films. 

Although the ice has a certain permeability 
for gases, it is justifiable to say that ice consti- 
tutes a rather tight barrier for the gas penetra- 
tion. In a case like the radio-carbon dating of 
ice, this tightness to gases is of importance, 
since there must be no exchange of carbon 
dioxide between the atmosphere and the 
enclosed bubbles. Estimations based on the 
presented data show that this holds true, as 
exemplified by the following situation: A flat, 
cylindrical, horizontal gas bubble, 0.5 cm in 
volume and with a top area of 1 cm? is situated 
in the ice 1 meter below the surface. The 
bubble has a pressure of 1 atm. and contains a 
gas mixture with 0.03 per cent or 1.5 - 10-4 cm$ 
carbon dioxide. Assuming solely a perpendi- 
cular permeation through the ice, and using 
Fick’s law with a permeation constant of 
10-11 cm? sec! atm-!, one will find that ro-® 
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cm? of carbon dioxide permeates from the 
atmosphere into the bubble in 1,000 years. 
This means that the enclosed carbon dioxide 
would be contaminated by about one percent 
of atmospheric carbon dioxide. With a per- 
meation distance of 10 meters and a bubble 
pressure of 10 atm. the contamination would 
be about one hundredth of a per cent. With 
regard to the dating problem, an exchange of 
this order of magnitude is so small that it 
would be a negligible source of error. 
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Note Concerning the Influence of Rotation on Convection 


By VICTOR P. STARR, Massachusetts Institute of Technology 


(Manuscript received January 9, 1959) 


My purpose in writing this note is not to 
present a model of an actual convective process, 
but instead to work with a particularly simple 
hypothetical scheme with the intent of illus- 
trating the action of one single factor, namely, 
rotation. The material presented below may 
therefore be thought of as a scientific parable— 
a vehicle for conveying a thought. More de- 
tailed models embracing the effects of other 
pertinent factors in addition to rotation have 
been presented by many persons (see, e.g., 
SALTZMAN 1958 and references there given). 
There is much to be gotten, however, from 
analyses which focus one’s view on individual 
items singly, to the extent that this can be 
done. 

The work described herein was completed 
many years ago, during an early period of my 
career. Accordingly, the general outlines of it 
follow along the trend of my interests at the 
time, which ran to the study of the readjust- 
ment of simple fluid systems under the con- 
servation of vorticity (see STARR 1939). In 
actuality, the treatment of the present topic 
depends upon three conservation principles— 
conservation of mass, conservation of vorticity, 
and conservation of energy. 

We consider an ideal fluid layer of uniform 
depth Dy and density o constant throughout. 
This layer lies next to the plane ground surface 
beneath another layer of ideal fluid of great 
depth as compared with Dy and of constant 
density 0'>o. The presence of rotation is 
allowed for by assuming that the Coriolis pa- 
rameter, f, is constant and uniform. We thus 


have an unstable discontinuity in density in 
the initial state which, moreover, is a state of 
rest with respect to the ground. The problem 
now is to say something concerning the shape 
and size of the finite disturbances which might 
develop at the internal boundary of the two 
fluids, if the latter is slightly perturbed. The 
entire point of this paper is contained in the 
suggestion that: No given hypothetical disturb- 
ance can grow to finite size spontaneously from 
infinitesimal beginnings, if its growth requires the 
violation of one or more of the three conservation 
principles for mass, vorticity and energy. Our 
technique is to assume that mass and vorticity 
conservation is fulfilled for each finite disturb- 
ance selected, and then to examine the fulfill- 
ment of the energy requirement. Other tech- 
niques could presumably result from other 
choices in this regard. 

Let us suppose that there exists in some 
appropriate sense a size-parameter related to 
the horizontal dimensions of the disturbances, 
which is significant for our problem. Let it 
be supposed. further that we 36 not run afoul 
of our energy criterion for a certain range of 
the size-parameter, but that beyond a given 
critical point energy considerations are violated. 
As this point is approached, presumably a state 
of indifference as to stability is likewise ap- 
proached. Thus if some such behavior is pres- 
ent, then in seeking the critical point we are 
justified in using equilibrium concepts in 
regard to the hydrostatics of the system, and 
also in regard to the geostrophic balance as 
will be seen later. 
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As a simplest trial, we consider a finite 
disturbance of the form 


D=Dy+ À cos Ty (1) 


where use is made of a cartesian coordinate 
system with the x- and y-coordinates hori- 
zontal and z vertical, while D is the deformed 
depth of the lower layer, A is an amplitude 
factor and L is the wave length which now 
serves as the horizontal size-parameter. It is 
clear that (1) represents sinusoidal waves inde- 
pendent of the x-direction. 

For a unit column it is easy to show that the 
change in potential energy occasioned by the 
presence of the disturbance is 


= g(e' - 0)(D* - DD (2) 


g being the acceleration of gravity considered 
uniform and constant. This is the amount of 
potential energy released by the disturbance. By 
use of (1) it is possible next to find the area 
average of this energy release, namely, 


Z 


I I ñ 
APE= > [32 -0(D*- DD dy= 


= 3(¢'- 0) A? (3) 


We observe that due to the sinking of heavy 
fluid and the rising of an equal volume of 
light fluid, there is a lowering of the center of 
gravity, resulting in a decrease of potential 
energy, which is independent of the wave length 
L 


The equation of motion for the x-direction 
may in our case be simplified to the statement 
that ' 

du y 

dt = fi v =. ae (4) 
the pressure term being zero due to the uni- 
formity of the system along x. Here u and v 
are the x- and y-components of particle velocity 
and f is time. An individual integration with 
respect to time between the undisturbed (ini- 
tial) and disturbed (final) state may now be 
made so that 


u-uo=f (Y-Yo) (5) 
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The significance of u, is the value of the x- 
velocity component in the initial state which 
in our model is zero. Likewise y, is the initial 
and y the final coordinate of the same individual 
particle. Equation (5) is in fact tantamount to a 
statement of the vorticity principle, namely 
that the potential vorticity remains unchanged. 
Its equivalence for our problem with this 
latter principle may easily be shown by differ- 
entiating it and combining it with the conti- 
nuity equation in the form given later (see 
equation (14) below). 

The equation of motion for the y-direction 
in the lower fluid is 


eae ES @-@ aD 
ep ee ine 2 (6) 


where p is pressure as will be discussed pres- 
ently, all motions and hence the disturbances 
of the pressure in the upper fluid are neglected 
for reasons of convenience. According to the 
borderline or marginal stability arguments 
already made, we are justified in replacing the 
pressure by its hydrostatic equivalent as shown 
in the second equality. It is to be noted that once 
this is done, the y-component of the pressure 
gradient is no longer dependent upon z. When 
this fact is combined with notion that in 
accordance with marginal stability concepts 
the left-hand side of (6), namely, dv/dt may be 
made as small as we please, if follows that u 
and therefore in its turn y —y9 is independent 
of z, according to (5). 

We may thus at once write the equation of 
continuity under the several conditions of 
our problem in the simple form 


do 
since columns initially vertical remain so in 
the disturbed state. It should be kept in mind 
that here again it is the initial and final states 
which are compared in quasi-Lagrangian fash- 
ion (see STARR 1946). 

We may now proceed to examine the con- 
cepts concerning kinetic energy implied by the 
foregoing material. Although it is by no means 
a necessity but merely a convenience, we shall 
neglect the kinetic energy of the upper fluid in 
toto. This is justified as follows. Marginal 
stability considerations point to the smallness 
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of the y- and z-components of velocity not 
only in the upper but also the lower fluid as 
well. As far as the x-component goes, the 
analogue of equation (5) for the upper layer 
would involve only small values of y-yo 
because the divergences due to the disturbance 
would be spread over a great depth. In the end 
it would result that the u kinetic energy is 
small. 

Giving due consideration to what has been 
said, we now write the expression for the 
kinetic energy of the x-component of motion 
in the lower layer as follows, with the aid 


of (5). 


I I 

5 ew D== of Dy — yo) (8) 

Just as in the case of the potential energy 

change, this quantity is for a given unit column. 
The continuity equation (7) becomes, upon 

elimination of D by (1), an ordinary differen- 

tial equation for y, in terms of y, i.e., 


dy, A 27 
eat (9) 
whose solution is 
al ae 2G 
Y-y- Dir À (10) 


where from symmetry considerations we drop 
the constant of integration. 

We now insert the expression (1) for D 
and the expression (10) for y-y, into (8) 
and average over one wave length, so that 


L 
I I 
AKES+ | ZefDiy- por - 


2,074 L? (11) 


Upon appeal to the corresponding expression 
for APE in (3), it is seen that whereas this 
latter quantity is independent of L, the ex- 
pression just found shows that A KE increases 
as the square of L. Since, as one should expect, 
the amplitude A occurs squared in both quan- 
tities, and since all other factors are constants, 
the principle of conservation of energy carries 
the implication that values of L exceeding a 
certain unique limiting value, L,, are not 


PESITARER: 


admissible. The critical wave length thus 
turns out to be 


27 a) 
ree ey gee 
PV ae ee 


If a diagram were to be made using energy 
per unit area as the ordinate and wave length as 
the abscissa, a graph of APE would appear as a 
horizontal straight line, while that of AKE 
would constitute an arc of a parabola in the 
vicinity of L. where the two curves intersect. 

We have thus achieved a demonstration, 
without resort to linearization, of the fact 
that rotation produces an inhibition of large 
convective cell sizes in our system, and it is 
rather obvious why this is so. The kinetic 
energy developed is uniquely determined by 
the kinematics of the deformations. So also is 
the potential energy released. However, the 
kinetic energy increases rapidly with cell size, 
while the potential energy does not, depending 
upon the amplitude merely. 

Actually, for cell sizes slightly smaller than 
L, there is an excess of potential energy 
released, over and above what is required to 
set up the prescribed kinetic energy of the 
x-component of motion. It is to be presumed 
that this excess is used to generate significant 
amounts of kinetic energy of the y- and z- 
components of motion, thus implying the 
presence of stronger instability for smaller cell 
sizes. Probably the strongest instability would 
manifest itself for wave lengths tending toward 
zero, although in real fluids viscosity no doubt 
exerts strong damping effects in this vicinity. 

Many variations of this demonstration can 
be, and in a number of particulars have been, 
made by myself and by my students. The 
kinetic energy of the upper fluid can be taken 
into account. The parallel rolls represented by 
(x) may be replaced by a square array. The 
concentration of the instability in a disconti- 
nuity may be replaced by a continuous grada- 
tion of density. As may easily be guessed, 
most of such alterations increase vastly the 
complexity of the manipulations, without 
adding much that is vital to the philosophical 
import of the results. 

It is not without a certain utility to return 
to equation (6), and to verify in a posteriori 
fashion some of the assumptions made pre- 
viously. Let us impress, in explicit fashion, 
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now, the stipulation that the motion corre- 
sponding to the limiting case when L ap- 
proches L,, should be geostrophic, by dropping 
dv/dt. After differentiation we have 


(13) 


Through differentiation of (5) and substitution 
from (7) it follows that 


(14) 


Equation (13) may now be written in the form 
gDy @ - (4 ED 
fi 0 dy? 


whereupon it is seen that we may actually 
derive an expression for D of the form (1) 
provided, of course, that L=L.. 

From a more general viewpoint the material 
presented suggests a number of points for dis- 
cussion. Some, at least, of these topics are the 
following ones. 

1. The kind of size criterion developed here 
probably has little significance for cumulus- 
type convection in the atmosphere. Such con- 
vection represents much too small cell sizes to 
be influenced by the earth’s rotation in the 
manner suggested above. 

2. The larger convective units in the atmos- 
phere (cyclones and the like) are more nearly 
of a scale for which the rotation is apt to be a 
governing factor from the energy standpoint. 
However, due to various altered and more 
complicated circumstances in the actual atmos- 
phere, any computation of actual magnitudes 
from the theory here given, for purposes of 
comparison with observed data or with results 
of other theories which include other effects, 
is not justified and could be misleading. 

3. Basically, it is the development of hori- 
zontal circulations due to divergence and 
convergence, in conformity with the Bjerknes 
circulation theorem which energywise brings 
about a profound difference between the 
rotating versus the nonrotating case of con- 
vection. These additional motions require 


DD; (15) 
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energy fot their generation, which energy fur- 
thermore increases rapidly with cell size on an 
area-average basis, and leads to the results 
obtained. (Compare STARR 1955.) 

4. A number of more sophisticated theories 
of convection deal with the steady motion of a 
fluid which is maintained against friction. In 
cases where rotation is present the horizontal 
circulations are again present and now it is 
their maintenance against friction which con- 
sumes an extra amount of energy. 

s. The notion seems to persist in some 
quarters that convection necessarily depends 
upon fair-sized departures from hydrostatic 
equilibrium, and that therefore one is not at 
liberty to use the hydrostatic principle in 
studying convective processes. This vis not à 
justifiable view. Actually any system in which 
the relative amount of kinetic energy of vertical 
motion developed is small, is sufficiently close 
to hydrostatic equilibrium to permit one to 
use this principle. In such cases when growth 
of kinetic energy occurs, it is the horizontal 
and not the vertical unbalance which plays an 
essential part. The large motions of the atmos- 
phere possess this character. 

6. The gross significance of what has been 
said does not reside cruicially in the precise 
fulfillment of the several approximations made 
for dealing with an idealized model. For con- 
ditions which depart somewhat from those 
postulated the corresponding expressions for 
energy, etc., would no longer possess the same 
simplicity, perhaps, but their fundamental 
nature would still remain about the same and 
lead to comparable physical conclusions. 


REFERENCES 


SALTZMAN, B., 1958: The general circulation as a problem 
in thermal convection. Scientific Rep. No. 1, M.I.T. 
General Circulation Project, Contract AF19 (604)-2242, 
416 pp. 

STARR, V. P., 1939: The readjustment of certain unstable 
atmospheric systems under conservation of vorticity. 
Monthly Weather Review, 67, pp. 125—134. 

— 1946: A quasi-Lagrangian system of hydrodynamical 
equations. J. Meteor., 2, pp. 227—237. 

— 1955: Proceedings of a study group on the general 
circulation held at che Institute for advanced study, 
Princeton, N.J. (in press). 


BOOK REVIEW 


MORSKOI ATLAS (MARINE ATLAS), Vol. Il 
Fiziko-Geograficheskiy (Physical Geography). Editor 
for Vol. U, L. A. Demin, Izdanie Glavnogo Shtaba 
Voenno-Morskikh Sil (Publication of Head Staff of 
Naval Forces), Moscow, 1953; 4 Sections, xviii pages, 
76 Plates, no price. (In Russian.) 


In 1956 copies of the “Morskoi Atlas,” Volumes I 
and II, were purchased by the U.S. Hydrographic 
Office from a New York bookstore. Volume II 
of the atlas is of particular interest because it contains 
a comprehensive treatment of the marine environ- 
ment of the earth. Although the atlas is not con- 
sidered an effort in original research, it does provide, 
in one volume, one of the latest compilations and 
analyses of many oceanographic and marine climatic 
elements on a worldwide basis. 

Other countries have also compiled and published 
similar atlases which contribute to the world knowl- 
edge of marine areas. In addition to the United 
States, other contributors are: the Netherlands, 
Great Britain, Germany, and Japan. 

Volume II of the Morskoi Atlas contains 76 plates, 
each composed of two or more charts or chartlets. 
The first 6 plates present the tracks for some of the 
historical cruises and marine expeditions. The few 
charts that are included showing oceanographic 
cruises are far from being a complete presentation 
of major oceanographic cruises. The next 9 plates 
are devoted to a treatment of the ocean bottom, 
its topography and geology. Plates 16 to 37 cover 
the physical and chemical oceanography. The tem- 
perature, salinity, and density of the oceans are 
shown for the surface and selected depths on s plates; 
4 additional plates are used to show cross sections 
for selected ocean areas and straits. One plate, 
composed of several chartlets, depicts ice formation 
for the Arctic and Antarctic regions. Twelve plates 
are required for current (7 plates), tidal (1 plate) 
and wave (4 plates) presentations. Marine biological 
information, including a treatment of marine birds, 
is contained in 3 plates, numbers 38 to 40. 

A total of 29 plates, Plates 40 to 69, are used to 
present various marine climatic elements; these 
inclutle presentation of the radistion-energy balance, 
air temperature and air-water temperature differ- 
ence, atmospheric pressure, wind, precipitation, 


fog, cloudiness, tracks of Lows and Highs, and the 
climatic zones of the earth. 

The last 7 plates cover related miscellaneous or 
general topics. The first two plates are devoted to 
terrestrial magnetism; the next two give a summary 
of techniques used in preparing cartographic pro- 
jections and a comparison of the dimensions of 
selected seas, gulfs, and straits. The last 3 plates are 
related to astronomy; two of these give informa- 
tion on the time of sunrise and sunset, and the 
duration of twilight, and the final plate presents 
charts and diagrams of the astral sky and the solar 
system accompanied by explanatory text. 

The entire volume is prefaced by five pages of 
general introductory textual material; 14 plates are 
also similarly preceded by explanatory text; however, 
these introductions appear to add little to the chart 
presentations. Usually the text is devoted to some 
rather fundamental or elementary explanations of 
the phenomenon charted, but, in some instances, 
exceptional individual occurrences (or observations) 
are cited in the text. All of the text as well as all 
chart captions and legends are in Russian; no 
abstract or summary in any other language is given 
anywhere in the publication. Since relatively few 
individuals in this country (the United States) have 
developed reading proficiency in Russian, the fact 
that all written material is presented only in Russian 
limits its general use. 

Recognizing that the atlas would be a valuable 
reference tool not only for the relatively few pro- 
fessional oceanographers and climatologists but also 
for investigators in the fields of geography, geology, 
biology, hydrography, and geophysics, the Hydro- 
graphic Office undertook a detailed evaluation of 
Volume IL. The technical content of each plate 
and its insets is evaluated separately in H. ©. Tech- 
nical Report 38; a comparison of each chart is made 
with existing similar charts from other sources.t 


1 (H.O. Technical Report 38 also includes a complete 
translation of the textual material introducing selected 
plates in the Morskoi Atlas, Volume II (Physical Geog- 
raphy). Translations are also provided for the Introduc- 
tion, Table of Contents, and other material in Volume I 
(Navigation and Geography) which may be of interest 
to the environmentalist.) 
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From the very first look at the Morskoi Atlas 
(Vol. IT) one is impressed by its size and the quality 
of material of which it is made. It is obvious that 
the U.S.S.R. has made a concentrated effort to 
prepare a first class publication, with costs appa- 
rently a secondary consideration. With a closer 
look at the individual charts, one can readily see 
the excellent quality of the lithography and the 
drafting and illustrating work. 

Probably the most serious shortcoming of this 
atlas is the complete avoidance of any mention of 
the observational data used, any comment regarding 
its quantity and quality, and the omission of a list 
of the major references utilized. Inclusion of such 
basic information is considered normal standard 
practice in scientific or technical publications such 
as this. 

Another deficiency appears in a number of the 
plates throughout the atlas: the omission of in- 
formation for the Arctic Ocean area adjoining 
the U.S.S.R. In most instances, this area is covered 
by an inset chart showing environmental condi- 
tions which may be somewhat related to the larger 
chart(s). These inset charts are usually similar to 
presentations appearing elsewhere in the atlas at a 
larger scale. The Eurasian landmass, which usually 
contains no data, seems to offer a more appropriate 
location for such small scale inset charts. Information 
is also frequently not shown for the Baltic, Black, 
Caspian and Okhotsk Sea areas. In some cases such 
omission may be unavoidable because of sparsity 
of observational data or because the chart scale 
may be too small to depict adequately the average 
gradients that occur. However, these voids impair 
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somewhat the overall comprehensive value of the 
atlas. 

The last general criticism, while of lesser signif- 
icance than those previously mentioned, is that a 
few of the plates contain inset charts that appear to 
be unrelated to the basic larger scale chart with 
which it is presented. Also, in one or two instances, 
the order of the chart presentation might be more 
appropriate if the sequence were changed. Some of 
the smaller charts might be much more useful if 
given at larger scales and smaller analysis intervals 
were used. 

It is recognized that no atlas could meet every 
requirement when designed to meet a universal 
need. Aside from the three general comments men- 
tioned above, this atlas is considered a rather ex- 
cellent summary of most marine environmental 
conditions. 

It will be nöted in the use of this volume that 
several of the charts resemble closely charts included 
in publications by “western” countries. 

In conclusion, it is important to recognize this 
publication as an index of the U.S.S.R. interest and 
effort in their study of the ocean environment. 
From the quality of material included in this Atlas it 
becomes obvious that U.S.S.R. meteorologists and 
oceanographers are well aware of work done in the 
“western” world. 

Haroıp W. Dusacu.! 

23 December 1958. 


1 The opinions or assertations contained in this article 
are the author’s and are not to be construed as being offi- 
cially endorsed by the U.S. Navy or Department of 
Defense. 


Current Data on the Chemical Composition of Air and Precipitation XVIII 


(For further information see Egnér, H., Eriksson, E., Tellus 7, pp. 134—139, 8, p. 285 and 517, 
10, p. 171.) 
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Hô 143| 140] 370] 54| 81| 210) 35) 31| 50| 4.6 fe) 

Sm 1910250225 6 AT 7 4|— 12] —4.6 o 

ÔT 54| 62] 94] 17] 8) 57| 30] 19] 85| 5.9] 30 

Sy CAN 30 AI Ame SIN 34 5:3 o 

BH 24| 35| 63. 12) 14|, 34| 5) 7, 24) 44 oO 

On 47| 44| 110] 13] 41] 56] 8| o| 25] 48 o 

Sk 38] 66) 220] 19| 28] 130] 15] 16| 34] 4.8 o 

Al 57| 90] 340] 26] 36) 210) 21] 36] 61] 5.4 o 

Hi 35] 24) 140] 13) 14] 72] 9| 11] 43) 4.8) 0 

Ta 33| 30] 200 5 Oro 2172| 2715|2.20105%2 o 

An NM ON EN EN — 

NE 45| 23] 200 3 4| 120| 18} ro} 35] 5.5| 18 

Gj 87| 39} 630] 10] 17] 350] 28} 48] 160] 5.9] 38 

Fn 7212.13] 64 5| ıı 38| 23 4| 79] 6.2} 93 

Fa 81 9| 84 2 6 x x x x 06.4047 

va TOP Tl 27T 5| ıı 6 6 2| 2236.31 eae ss 

ET 12| 15 10 3 5 5 4 77 red) Cee] 32 

Ke 45| 24| 41 7 PM eh 812 4| 46] 5.8] 61 

Sd 107| 24| 470 9] 26| 260| 20| 35| 43| 5.3 o 

Da 18 oOo} 12 I I 2 4 Ties 0587532 

As 59| 19] 85] o| 12] 65] 38] 8! 37] 5.71 24 

Li 35 8|1 850 7 2|1 010) 56| 120] 73] 5.8| 19 

So 24 Et a 3 I 4 I o 9| 4.9 fo) 

Ka 831733 4 5 3 2 I} 18} 6.0 4 

Ku 15| 45 15 I x 4 8 1| II 4.9 o 

UY. 52] On] 58 3] 10! 9 9 92 10 48% 

Pu 55| 37| 28] 4) 11 St 5 Ol 8 54 o 

Tv 51| 75| 214| 18 4 89 15 9| 26| 4.4 ° 

Rj AZ| TAC 32380 TS Ar 27 o 

Vn 72| 108| 313} 18) 38| 190| 21| 48 8} 4.8 fo) 

Gr 38] 97| 1690] 15] 31) 86| 15) 22] 35] 4.3 o 

Öd 69] 121| 470| 17| 56| 262| 26) 48| 6| 4.9 o 
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mg/m? Sl 2 ug/m?® (=kg/km?) 
S| 10 

Sales lc ( Ê PH lil. i 

oO |i | Na! K|Mg|Ca © [RIG] S | Cl | À | Na| K |Mg|Ca 

Zia q Zz, 

Precipitation January 1959 (D oor) | Air January (L 901) 

Bs 73| 158) 1451) 18) 33] 599] 53] 93] 64! 4.8 o| 84 
Fj 78) 147| 2285| 16} 30lı 214] 70| 207 4| 4.9 ol 172 
NL 87| 22914 024| 33] 40/2 230| 186] 230| 192] 4.9 ol 172 
Ly 43| 120| 188] 20] 31] 104| 16] 28| 29) 4.2 o} 58} —| — — 
As 46} 80] 360! 15} 40| 220] 24) 41 3| 4.6 ol 51 
Vd 79| 188| 540| 28| 154| 290] 31] 59 8| 5.7 ol 46 
Bl 35 Ori) 228) 73] 271 125] EO] 36|. 18 4.7 o| 48 
Ty 20) 50, 2177| 1712| 200 226] 22| Tol, 22]| 5:1 o| 63 
Hö 54| 127|1 137] 16) 38| 583] 43| 73) 68) 4.9] 0| 91 
Ad 231 551, 188| £2) 26| Ero] T6| 21 Ton 4:7 ol 55 
Lw 26| 182/2524| Io 411 554] 64| 247| 85| 6.3] 126| 418] 5.9| 35. | 0.6] 18. | 12.5] 2.4|0.7 
Sw 128| 237|2 387 3 glı 341| 54| 163| 90] 5.4 ol 71| 5-9| 8.6| 0.8! 4.3] 1.3! 0.7/0 
Ab 30| 61] 419 4 7| 241] 15) 32] 30| 5.4 Ol) 1631202 | 17. | 03| 20.7| 28 2213 
Ed 35| 46| 166 Ay 271.961, 7785| 20), 251) 5-9 9| 30| 31 16 a SO rl | 
Es 901571220. 287], Tri) 241 284| ı6| 23| 27) 6.1| TS 2317 5-1] 1.0] 2.3| 1.2] 0.4|0.8 
Ag 28] 57| 250 2) Unie ra0 5] 22] 14| 5-4 o| 48 
Ike AEG 230, TO, STE 89| Iol 27| 33] 4.2 ol 56| 29 73113311 3:2| 0.8| 2.409 
Ro 64| 107| 248| 14] 18| 123] 14| 38| 64| 4.9 Or Be >= 
NA 140| 221|1 099| 15| 39) 939| 25| 74| 98| 6.6| 86| 46) 24. | 12 3:31 3:21 0.3) 3.5 2-8 
Ca 189] 674|7 180 8] 3614 032] 237] 464] 227] 5.6 6] 139} 9-1] 9.3] 1.2| 2-5] 0.5] 0.2/0 
Ma 70| 210|4 284 6 5 6.2 al 220 | — 
Bt 71| 156|2 691 6 6 6.2 22\l 1447| =} — | =| =| | —| 
Cl Gri Onn S05), 65h 3S 5.2 2) 3) TE 
DA AHIRTOS sro AIS 5.2 ol 49 
Bi 53} 69] 398) 74| 5 6.4| 20) 41 
Rl 93] 130|1 349 7 7 4.9 o 61, — —| —| —| —| —| — 
Va 119| 219|2 404] 18 8 62} 24| 85 
We 66) 132/4 606 5| 13/2 895| 112] 133] 178] 4.2 o| 247] 39.41191-4| 7-0|101.3| 4.9] 9.2] 8.5 
SC 61| 192| 632 4) 541 322] 35| 38| 67| 4-3 ol 56| 31.3| 30.5| 9.6| 18.3] 7-0] 3:116.5 
BV 251, 80) 173|, z2|) Aol, 50 8| 20] 56| 4.0 ol 60] 49.4] 10.7| 4-4] 3-0} 2.0] 0.3] 2.2 
Bn 49] 115| 146] 22| 23} 65] 1o| 48] 335] 5.9] 80] 60) 28.1] 35-5) 11-1] 10.3) 5:5) 0-9) 4-4 
Au 53| IO1 90| Io] 143] 25 7| 13) 58| 4-8 ol 21[26.0| 3:5| 4:3) Tt] 1.8) 3:1) 2a 
Fe 81] 101 o 4 ol 40 5 8| 32| 4.4 O7, 
Ba I20] I02| 24 o ol 36 7 8| 48| 4.4 o| 1578.9] 43| 2-51 1-5) %4| 0.2) 2 
Ho 83} 83 2 60.7310.301.39|, 17, 9710 5.0 o| 18[21.6| 5.3| 6.2| 1.4| 0.7] 0.2/1.6 
Rm Szlerrol, 82 fe) 2 201 10 gl 260] 5.9] ı30| 38] 8.5| 1-5] 2.5| 3.1] 1.2] 0.2/3.8 
Et 40| 72 8 2 18 9 55 3210425 o| 21| 26.6 o| 8.1] 3.1] 3.3] 0.6) 3.1 
He S|) en | 
Rz 6 3 5 Aly 23 2 2 ZONE ol 43[15.5| 4-9] 6.3] 1.8) 1.3] 0.2|1.6 
Wi 19| 56 o 3 I 4 3 | GSM 520 O|) 3514.01 27 278 3:2| 1.0 0.425 
167; 35| 61 o 8 7 6 3 9.42 5.2 © 191720). EO) 5.010.7.01,.7.41 0448 
Kl 8 8 o al 5 2 2 Ol 205.1 ol 27| 291 2.7| 4.0) 1.3] 0.9] 0.4| 1.0 
BL HA SOE | 209 8} 121 25 ZI 2710. or] 40 ol —| 5-9] 4-2] 0.5] I-4] 0.9] 4.4|9.6 
DH 99] 167|7 029] 26] 26/3 663] 132] 337] 216] 4.0 o| 255] 24-4] 29.7| 2.7] 17.9] 2-9] 4.1] 1.4 
DB 79| 134] 580] 23] 20} 308] 17] 40] 39] 3.6 ol 67} 40.3] 2.9] 4.8] 2.5] 0.5] 1.2/0.8 
W 74| 193| 368| 23] 30] 185| 14] 28) 20] 3.8 o| 43] 20-4] 12.9] 3.4] 2-5] 0.5] 0.7/0.1 
SA 93| 306] 924] 24] 35} 502] 25] 69) 44| 3-8 ol 64] 28.6] 7.7| 5-1) 2.1] 0.6] 0.6/0.7 
U 78023462460 AAA 882512231 5515 339 ol 44|85.9| 11.5] 7.6| 2.1] 1.2] 2.6|0.9 
B 1ı15| 136] 214] 30] 16] 105 9| 15) 23] 3.8 ol 27|15.0| 12.3] 1.8] 1.9] 0.6| 0.4| 1.0 
D SA ao 17070 7710 740295 Wl ee 27103:9 o| 31] 32-8] 4.5] 4.9] 1.8] 0.9] 1.3] 0.5 
MH 77| 182] 206] 16] 11| 101 8517528842 ol 24] 38.5] 11-1] 3.6) 1.9] 0.4| I.2| 1.0 
Rs 135| 32/1 199] 12] 14| 675| 23]. 82| 36) 4.8 002351772010 ONE 0,5) CS) TA 
LM 78| 66| 255| Io 5122 5| 25| 23| 4-9 o| 18 : 
Lx 73518. 65107305 23|) 14] | 82 Al 2x03) 2710 47, ol ıılı5.5| 1-3} 2-4] 2.3] 0.5] 0.0|0.5 
Bg 87| 46| 204 9 6| 94 6| 16] 32| 4.9 ol 17|24.6) 4.7| 3-1] 3-8] 0.7| 0.5|0.2 
Ae 83] 97| 85| 19 7| 59 7 LO OZ 00523 ol 17| 21-4] 0.0 1.6] 0.7] 0.6] 1.1 
Za x x x x x x x x x x x x x x x x x x| x 
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2 mg/m? 4e ug/m? (=kg/km?) 

9 Sl We | 7, 

= Z A I.» lh|o 

mm | | rule aude Fa 

= S| Cl] @] | Na} K| Mg} Ca © |xsic]| S | Cl} A |Na| K | Mg] Ca 
3 2 | 2 je 7 

o- 

7 Precipitation February 1959 (D 902) | Air Fecruary 1959 (L 902) 

) 

Ri 41} 30] 460 2 3| 320), 236]. 371: ki 
Ki CRE (6) 2| 78 4 3 71 FOX 
Ar 5 e el il a ai A 0.7| 0.4| 2.6 
Oj As | Al, al == 1.4] 0.7] 3.2 
Rö 9 TI eT2 o fo) 8 3 3 0.7 (0.3) 24 
Of zu 310 ero fo) fo) 6| 14 4 1.2]: 0.3 x 
ÄF DD GAIN LAN CC GR AINSI IE SEE 
Fö 29| 44| 52 I Ti) 420 8 5 — —ı — 
Rä I Al ah SUN ENS 0.5] 2.5) 5-1 
Am 2 5 — —| — 8 5 — 2.7| 1.1] 9.3 
Sa 5) 14 By SS 5 2 Bf 0.7| 0.9) 7255 
Ul 5] 27 SIE FSS 2 |) = : : 

Er 22| 30} Io 5 6 8 5 2 

St 8 37 8 3 2 6 4 2 

Fo 3 ar a et 4 a 

Kvı 28 16 3 OBERST 7 

Kv 7 6} 22) 15 3 7 8 5 2 

VK 2) 12 US 7 2) — 

La Sr sel a in NC Sp —- 

Bo 5} 30) —| —| —| 120/ 10) — 

Fa 2 3) SS) S| ST CS 

Fl 7 4 6 2| 81 8 6 3 

Am 10) 16] > 31 else male 6 

Pl 9| 36| 1490| 10] 19] 85| ,7| ı0 

Hö 21| 65| 84| 22] 47] 46 6 9 

Sm ol | —| — 

ÔT 51 20] 28] —| —| 15] ıı 

Sy To) 35) 27] 1 7|-413| 5] 4 

BH dl = — 

On 4 of 63 — — 29 5 — 
Sk 5| 32] 160 93 8) — 
Al 7\ 249 1750| WAL 37 = 82 go) 25 

Hi 122 ie 4l — 
Ta 31| 37| 780 2 4| 480] 30] 58 

An = 


Va 30] 13 ZONE ET, 4 33 2 

Tr 271 35 501 fa) FSi" 38] IS 

Ke OW Sale a sit ale ae 

Sd 105} 35} 820] 18) 38) 460} 40] 56] 35] 5.2 
Da 17 6| 23 5 4| ı2| 25 4| 66] 6.7 
Äs or PRET a Te Fe ES En 
Li 29 2|1 330| 23| 13] 600 4| 87778 5.9 
So 6 8 48 I 3 5 I 3 5| 5-1 
Ka 17 co) 70 I I 6 2 ATI 5 
Ku 7 9} 78 I 3 6 3 4 6| 5.0 
Jy 18 7| 45 I I 6 5 2| TO! 5.7 
Pu 2) hl ad I 5 5 I 2 val Set 
Tv ES} EO) 20777 2 X 7 3 4 7| : 4.4 
Rj 125 318 040 3 1/3 580} 26| 600] 184| 6. 
Vn 6| 45] 1o1| 7! 13] 52] 15] 14| 49] 5.9 
Gr 4] 33) 68. 7) 11| 36) 6) mB 24-44 
Od 4) 30! 82) ro) 24) 42] 5] 7] 9 47 
Bs Ol) 26 “real al 99 gt 18! x4! 42 
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2 mg/m? E- ug|m® (=kg/km?) 
3 Te Tat eee ee ee a em 
mm wes poet ey | x 4 

= S | Cl] & | ow] Na} K|Mg| Ca © frig} SE CI cof Na| K | Mg| Ca 

5 yA cn Z 

8 

un Precipitation February 1959 (D 902) | Air February 1959 (L 902) 
Fj 28) 59] 383 17] 52] 1| 4.6] of 74 — 
NL 11| 18|1 503 25| 791 33) 47 of 36444 —| — —| - aH — 
Ly 5| 38 6| 16} 44] 4.0 o| 189 
AS 27 9! Ah Ad 10) 284) D 
Vd 4| 48 Bl a FINS ol 177) | — —| — | | — 
Bl ah) 49 3] 5] 18) 4.3] o| 134 Sal a 
Ty 21. 26 4 ai Tèl 47 ol 243] — — SS — — —| — 
Hö 3| 23 5) 14) 8) 4.1 0| 246, — = 
Ad Bln 22 7 6| 14| 4.8 O2 359 Eh = 
Lw 24| 153 64| 258] 36| 7.2] 936] 407| 5.2| 9.6] 0.8] 3.8] 0.7] 2.9 o 
Sw 28] 140 Sr) 1721, 781.57 21 288 77.110345 7.7838 1272:01102|,5:0 
Ab Al 22 4 ON a oe > 120 4.71 4.3| 2.9] T:0| 2,6 o 
Ed IO} 33 8 9. 1g}, 6.0 8| 52|24 8.710 5:01. 5.0.1782 T7) 22:0 
Es 62| 168 SO Sala su 5-4) 272756176 Re Zul N Bra |) Aik: o 
Ag 21| 68 81.161.351, SAME ME | 22 TO. 7. 35512 3,21 27:00. 23 fo) 
Le ı 45 ea a SS ee 8.2| 8.3] 4.2] 3.3] 7-0] 4.8 
Ro o 2 == I 51 780 CES wesGll er 122]Wersr 
NA 12| 169 26] 1212 7315 0.8| 5377|. 198379] 22 am 0:215.2:31533 o 
Ca 20| 678|4 097 280| 419| 264| 4.0] —| 988| 13 872431232 27:31.70:8 050 
Ma 21| 139|2 367 6.4| 47| 415 a ei — 
Bt 30| 117/2 604 = Ge ES SA NS EN MN 
CI 35| 39 — —| —| 6.3] 13l 34 —| —| —| AH — = 
DA 13] 72 0222) Tar O — — 
Bi 29]. 26 —| —| —| 69 96 400 À — A] — — —| — 
Rl EN Ey Be — —| —| 5-4 ol 282} — —| — — Aa 74 — 
Va 47| 71 SS See) dgl 18cı — ee a 
We 4| 67 30] |, 591 3-9 olı 185] 15.0| 13.7) 3-4] 6.6| 0.9] 1.2] 1.6 
Sc 4| 30 AM Eee Ol) 304 2700| 12.412721 2. 3.0) 2:20 2.7102.2 
BV Ze 8 Sar Oly 3:7 ©| 609] 34.0| 12.2) 9.0) T.7| 1:71 1.3] 1.7 
Bn 2\ = 7310 102751042 o| 925| 58.8] 33.1] 17-6] 4.4] 4-4| 6.0| 26.4 
Au Os — — —| —| —| —| -—-|13.9| 4-3| 12.2| 2.0] 3.1] 7.9| 2.0 
Fe 32| 43 Io 6| 16| 3.8 A1 Age) = 
Ba 7 24 5 P54 3.8 ol %45|23.8| 13.8 55| 74 2:31, 2.7] 74 
Ho 1212 62 13 5| 28] 4.1 ©| 65|22.9| 3.0] 3:0] 0.6] 0:6) 2-4) 1.2 
Rm 4 3 Zu ee > GE INC) EN x x xl Ele xl ees 
Et 4| 13 1602 3 | 400 o| 165] 21.0] 8.5] 6.9] 2.4] 1.2] 0.5] 1.6 
He en ES SS Sa — 
Rz a aks 2\ al ee Leo o| 67| 15.6] 0.9] I0.0| 0.9] 0.6] 0.4) 0.6 
Wi 16| 15 6 9} 89} 5.8] 68! 59] 8.1] 7.3] 5.6] 1-4] 0.7] 0.6) 2.4 
Lz DE 71 12 8| 54| 6.1] 200] 117| 4.1] 5.1| 6.9] 0.7| 1.1] 0.6] 2.5 
Kl == — 2520| 13:5, 177.010 3:31, 0.710 Ley] 550 
so fee — —| 1.7| 1.6| 0.1] 2.7] 0.3] o.2| 0.5 
Ve tes —| — ale aah 
BL 19| 46 ETS) TOM AT o| —|17.1| 7.1| 6.2| 1.4] 1.6| 3-3] 6.0 
DH Ay = 51 —t, 33) 3-9 ol 650] 7.1] 42.9] 7.6] 16.5} 2.0) 5.9] 3.0 
DB 4, — al —| 27| 3.4 o| 304] 3-4] 24.4] 7.3] 2.5] 0.9] 3.9] 0.6 
W 2| — Al == 8| 3.6 0133118301 06 874) 4d 146. O35 | O58 
SA 44 — 6| —-|..18|. 3.9 o| 247| 7.8| 30.9| 11.4] 1.0} I.I| 0.3) 0.7 
Ü 5| — 6| —| 40] 3.8 o| 236| 9.5| 0.9] 10.0] 1.7| 1.7| 9.9} 0.9 
B 171836 5| —| 25| 4.2 ol 67| 0.1] 17.1] 3.0] 0.8} 0.6 o| 0.7 
D 6| 39 57 | 201s 4.7 ©| 125] 1.8] 20,0) O16] 1:2), 1.5] 0:9) 1.0 
MH Il — —| —| —| 3.6 o| 1215| 44.4| 16.6| 9.5] 1.4) 1.5] 1.4) 0.4 
Rs 9] 31 SOIN 10) 5:5 ol 74] —| 22.9] 4.7] 3.9] 0.5] 2-6] 0.4 
LM 4, — 2] — 9| 5.1 ol 65) NES le le 
N 8 23 4A —=| 16) 5.6] =| 52) 2.6) 7-5! 3:6) 0.4 Qs, 01 08 
Bg I 2| —| 20] 5.0 ol 208] 5.8] 21.9] 6.5| 0.4] 0.8) of 0.4 
Ae I} — =| —| —| 6.21. —| 263] 0.5] 13.2] 5-3] 4) £-2) 0) TO 
Za, | 21l 2 7o| 32] 160] 6.7| 340] 89] x| xX| x XI] XI XI xX 


{ pd \ 
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K | Mg] Ca 


4 61 
2 IIO 
4 fo) 
5 fo) 
6 fo) 
4 Oo 
6 fo) 
6 fe) 
19 94 
9 o 
7 fo) 
8 fo} 
20 fo) 
7 fe) 
6 fo) 
7 24 
22 fo) 
Io fo) 
4 fe) 
6 oO 
5 fe) 
9 fe) 
7 fe) 
5 fe) 
5 o 
3 32 
31 fo) 


w 
© 


Ww 


Ny 
Oo OO NY 6 OO 


a NN 
Œ H ND ND OO HJ DN ONAN NWI 


ORONO: CO OS 97 O16: OO. 0.0 


Tellus XI (1959), 3 


CHEMIGA LT DATA Sr 


2 mg/m? FIB pgim® (= kg/km?) 
8 zZ 2 2 | O15 7, 
„ | mm à , 2:9 0° I, 

= S| cn Oo | | Na| K | Mg] Ca Oo |R[G} S | Cl] Æ | Na] K / Mg] Ca 

3 A |A je A 

a Precipitation March 1959 (D 903) Air March 1959 (L 903) 
Fj 2662| 2317201 3417 716 8 Al 27) 4.3 Ol EN ff = = 
NL 54| 114] 137| 32) 48) 72] 19) 18] 40] 4.44 of 43) —| —| — 
Tey: 16} 81} 34] 44] 26| 13 8} ro} 58] 4.0 Oo] 94 = 
As 35| 203] Tog! 25) 58) 38] 14] 20) 34) 4.7 Ol ig, S| 
Vd BrP E20) 3217 25 87 15| Io 8} 59] 4.9 o| 49 —= 
Bl TA atte RES ete ta]. „ON LO) 3704 Ol gO) | Sy = 
Ty 200. 78 231 al 5727| 2816 8 67| 5.3 el ee = 
Hö aa SI aaa a5 org era | à 
Ad 72 RS ONE NTI 2 34 17 21 5.8 SIN 7 un 
Lw 33| 12411514] 26| 16|1181| 60| 165| 43] 7.2] 823| 259| 4.6| 6.7| 0.9 0.5 
Sw 30| 67| 691| 14 8| 413| 15| 58| 24| 4.9 o| 100} 5.0] 2.7| 0.7 = 
Ab 32| ro8| 357| 33) 32] 235] 70 35] 47] 4.7 || 76 zig || 8:818232 = 
Ed 7438| 77 Siero meas 8} 10) 29] 5.9 Bile a4 Sire 2 720 05:9 I 
Es PEO |2r2 7100 20 431.262 2571242) 361 4.7 IN 351 13-)|| 2.2] 7.9 —— 
Ag SA SE 7621 2531 7750. 1785| ro) = rr 241 6.5] 122) 745118. | 9.61, 5.2 o 
Le 24| 95] 173] 19] 36) 64) 11] 26| 78) 4.5) —| 7339. | 5.7| 4.0 cat 
Ro S7ZEr632260 0.271 535187942 23172270T02| 77, 17 248534.) 6:91 3-4 > 
NA 1720787 ,837| 2ı| 45| 497| 228) 65| 112) 6.0| 1838|. 39] 79.7] 6:0| 2.4 1.0 
Ca 107| 283|2 712] 42] 52/1 562| 86| 2ı9| 128] 5.6 ill, TOE) GO|) dah bes) == 
Ma 84| 302|4 998 7 Ol u SS AO) Bas SS E— — 
Bt 90| 24313 105 9 =) ee GA 20 10773 8 INR — 
CI 79| 134] 593] 7| 6 320 Olt UNE SNS 
DA 44| 110| 422 4 AN || NA NO Mere 
Bi 88| 114| 343 8 Same SOIENT ee == 
RI 62| 167|I 041 6 A S| el Al Fa =| =] 
Va 109| 185]2 354 9 9 5.5 o 555 = 
We 51| 204] 920) 33| 38| 427| 23| 65) 69] 3.1 ©} 79] 7-4]14.7| 3-3 1.4 
Sc 307710671028 44) 371 11 9} 98] 3.9 o| 63] 32-7] 25.7] 6.6 1.5 
BV 72082318. 671 2271 361 20 in| 29198) 4.8 o| 169] 52.7| 5.8| 6.8 27 
Bn 44| 22! 226| 49] 64] 20| 26| 72] 753] 6.1] 186| 134|87.2| 80.3] 7.7 2.5 
Au ZS Ae 52102817 AA eons 10) SsAiler7ol) 15.2 0 167|/ 13-0] 0.0/5.8 287] 
Fe 41| 79 ol 24| 34 8 11 518261 3.8 03302 =| sa) == — 
Ba 13040 28/4707 19) 28 T2 0 60379 ol 24| 0.0] 9.9] 4.1 0.5 
Ho 43] 106] 13] 20) 19 Sih ll | le 08.204841 5.3102.7 0.6 
Rm 69| 41 94| 76 oo er ae] all ee) Qe. PL Basis Ss) Oxy x 
Et 5Al227 17727122817 31 TOT AIS? 428 ON ay) ee == 
He oa S| en =F 
Rz 29| 70 o 8 5 2 4 9} 50] 4.6 ol 21) 15.9] 7.2| 4.5 0.4 
Wi 32237034 8 2 SNS zur ur29l085.7102581, 321177.3|| 6.1105.0 2.9 
Lz 27 | ane? SRE 4 4 Wer GRE bya 02 32|(20:2)// 11-4] 6:9 2.8 
Kl 53| 65 o 6 9 3 6| 21| 121| 5.6] 44] 22]20.0|16.6| 4.6 0.9 
SO 46 ol 45 67 22 17 14: TAT! 25.838) 16) 0.0] 0.8034 177 
Ve 18| 44| 873) 12 I] 388] 46) 52| 166| 6.5} 228] 194| 12.8| 13.8| 3.5 27 
BL 2242| 25 3| 45 3 8 Io] 49| 4.3 Ol 1707318723 Oya 7.0 
DH 538 £06) 262) © 4324) 5072) 15) 41 33) 4:2 Oe 7512-91) 2.8) 97-3 1.4 
DB 54| 145| 85] 33| 26| 41] 14| 21] 72| 4.0 0/02.73109:81884 07:2 0.5 
W A0 A Ar 23] 3410 22/0511) 271074 ol 51] 0.4] 6.3] 5.1 0.7 
SA 49] 143] 94] 26] 72] 48 9] 13] 62| 4.8 ol 42] 6.0] 6.7] 8.7 1.0 
U 74| 163} 86; 38] 35] 33] IO} 33] 93] 4-5 o| 36] 19.2] 10.7] 6.9 0.8 
B 83] 73) 62] 26) 13 28 6} 18] 34] 4.8 Ol ||) Be 1:O 
D 52 OO? 5 0.20.7181 20 (ED a3 2 ares 0923183:09 bees IN I.0 
MH 41] 69, 57] 25} 31] 36 9| 32] 18| 4.5 OSSI 2.21 06.7 | er. 5 0.5 
Rs 891 2231 331 34 4| 26| 14) 51| 23] 5.3 Ol 34| 2 174:.010.2:6 1.0 
LM 46| 49| 78| 10 7| 66 6| 22] 26| 5.1 O| 19 == 
Lx oeil EYES INT) IN 771 1.7 9] 24) 37| 5-5 OT =|) Sill So) 0.5 
Bg 46| 86) 33] 17 | sul 20), Al GA SC Gl AN GATE 0.9 
Ae 161| 105) 84| 53 al. 21 al GO) re AI AE NON) 4.2 0.7 
Za 117| 54 14| 401 94] 21| 23| 35| 310] 5.6| 31 39 ake ll 2:51 8.4 5.3| 25. 
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372 Changes in the Networks 


: Lat: Alt:t. ER : 
New Stations Code N Long. Greenw. ee Precip. | Air | CO,| Start. 

Denmark 

Aalarbeaye rc Aa 55°04" 1455’ E == x — | — | Apr. ı 
Store Jyndevad ....| SJ 54°54" 09°10’ E _ x — | — | May ı 
| Germany 

Konissthul mer ar. Kö 49°24’ 08°44’ E 590 Ss x | — — 
einen 7:00 Ln 49°23’ 08°39’ E 140. x X | — — 
Italy + 

Savze di @luxe.. m. SO 45°01’ 06°52’ E — x XL" Feb: 
Vesima te ter FT Ve 44 24° 08°42’ E — X x | — | Feb. 1 
Sweden 

Are Hummeln...... AH MET ARTE 
ABOLISHED 

STATIONS 

Sweden 

Ojeby nies car cick Oj coordinates, see x | Feb. ı 
Erken ve. cena Er Tellus 10 (1958) x | Feb. 
IDEN NG bans né oo Fl Da x7 x, | FER 
PIOMMINGS ep. EEE Pi XLR Pest 
Norway 

Vagamo.wr...ı. Va x |. Heb. 7 
Denmark 

Odum se. eher. Öd < | Feb. a 
Made TR: ern. Vd X May ı 


From May 1 1959 the CO, values are reported as ml/ıol dry air. The earlier chemical method 
gave the values in ml/ror H,O saturated air. 


Correction 


The air values for the Finnish stations in August 1958 (L 808) 
should be ug/m?. 


Code | S | CI | NH3-N | Na | K | Mg | Ca 
So 2.0 0.3 4:7 0.3 I.I 0.8 0.5 
Ka 7.1 1.5 ax. 2.8 5.7 2.6 1.9 
Ku fo) 0.9 23% 1.4 2.2 2.4 o 
Jy 19. 7.0 25. 3-4 4-3 7.1 10. 
Pu 3-8 fo) 10. 0.4 1.4 0.6 o 
Tv 24. 3-5 56. 6.5 13. 5.3 o 
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CO,-values in Scandinavia February— April 1959 373 
(Cf, FonseLius, Kororerr: Tellus 7, pp. 258—265) 


February March April 
2" nn | + 

O 4 . H . 4 
2 lo Le ls elle A ie o a 
8 || °c a = IS SiS] cc q 2 |S8l8| ec 2 -| 4 [se 
n- [m = © [9 ZA = © [Osa = FAN SE 
= = q = E 
17| —3-5| WSW 8 |x° ro] 3.26]|17 4 |NW 5 9| 3.63||14| —5 |E 1.5 o| 3.29 
Ao |18| —2.0| N 0.5|%° ro} 3.11||18 3 |W 9.5|@ 10|3.31||15 1.5INNE 2 73:31 
19] —2.5| WNW 8 |x° 5] 3.26/[r9] —3.2| NNW 5 |x 8 3.32|[16 1.2| ENE. 0.5|=° 0|3.27 
17 rT iow 18 10| 3.14|\17 2115 2 10| 3.73||14 8 |Calm 10| 3.19 
Br |18 5 |W 5 |x 8) 3-15}/18 6 |SW 3 ı0| —|I15| 12 | Calm 0| 3.22 
19 1 |W 10 | 7 IO} 3.13||19 LEN! 5 10| 3.23[|16| 14 |Calm O} 3.55 
140 20 (5 I 8| 3.34 
AH Wall wea |S I o| 3.39 
16| 15 | Calm o| 3.70 
17 2.3| W 2 2| 3.39||17 3.8] SW 1.5 TO} 3.25||14 6.115 2 o| 3.40 
Ka |18 5.11 S 0.5 5| 3.42||18 6.4| SSW 2 2| 3.61|[15 9.51 S 2 IO] 3.54 
19 4.0| W 4 10} 3.50||19 ZN I 2| 3-49||16| 13.6|S 0.5 10| 3.48 
17| —0.5| SW I.5|% 10] 3.74||17 1.8}SSW 3 2| 3-63||14 3.4| WSW 2 10| 3.40 
K |18| —2.2| Calm 0} 3.64||18 6.0] WSW 2 o| 3-41]|15 7.9| W I Io} 3.48 
19 1.4| W I 10| 3.70||19| —3-0] NNW 1.5 o| 3.28||16 9.4| W I IO} 3.38 
[17 2.4| SW I 5| 3-46]|17 3-0] SW 2 IO} 3.52||14 3.5|S 2 IO] 3.49 
Jy |ı8 2.1) Calm 5| 3.381118 4.5| WSW 2 2| 3-29||15| 10.6|S 2 5| 3-40 
a 3.5| SW 2 5| 3.62||19 2.2| N 3 10| 3.78|16| ı11.5|SSW ı 5| 3-48 
17 1.0)NW 2 8| 3.46|[17 4.4| WSW 2 5| 3.64||14 P|] SSW 10| 3.65 
Pu |ı8 1.4| W 2 Oo} 3.37\|18 6.0| W I 5| 3-43||15| 12.6] S I 3| 3-54 
19 1.8] WSW ı —| 3.56||19 2.5|NW 3 2| 3-39]|16| 14.4] SSW o| 3-54 
17] —2.4| W 2 10} 3.61||17 4.6] SW 2 [00° 5] 3.34||14| —0.9| NE T 100° 113.44 
So |18| —7.6|S 2 10} 3.31||18 3.81SW 45 10| 3.35||15| —0.8| E 00® 13.45 
19 t.alSW 3 10} 3.90||19| —4.1| W 3 5| 3.30||16 0.7| NE I |oo® [3.44 
17 40|NW 3 5] 3.66|[17 2.0] SW I = 16[3-.45||14 6.01SSW 2 IO] 3.40 
Tv |18} —4.0|S I o| 3-59)|18 4.0] SW 3 o| 3.27||15 8.0] SW 2 O} 3.53 
19 4:0|NW 3 O} 3-62]|19 6.0] NNW 2 o| 3.27||16 10.0| S 3 O} 3-33 
17| —3-4| NW 7 5| 3-33||17 2.8| SW 3 IO} 3.79||14| —11.0| S I o| 3.22 
Ta |18| —8.8 | SSW 3 5| 3.40||18 3 SW 003 y 10 3.65||15| —12.8]S 4| 3-42 
19] —0.4|SW 3 5} —||rgl —9.0| N 15 | 7 10] 3.36||16| —3.6 —|=° 10| 3.59 

e 

T7 15|SW 8 |# 10] 3.54/17 65ISW. 7 x. TO 3.36|[14 5.6] E I les 13.20 
Vg |18 6.6] SW 7 |= 10| 3.15||18 4.2| SW 7 | ° ro} 3.40||15 5.2| W 3 Se 5| 3-17 
19 1.8] WSW 6 |=° 7| 3.23||19 0.7| W 31309 13:321177 4:0|SW 4 |=°1ol 3.15 
17 4| SE 2 10} 3.06||17 6 |S I 5) 3.424) 13 |sw 5 12° 6! 3.43 
Sd |18 4|SW 13 . 10 3.26)|18 25 I Io} 3.32]|15| 18 |SE 5 O} 3.15 
19 4| W 13 |v 10|3:13||19 TS, 3 0| 3.3416) 13 |SSW 4 O| 3.17 
17 5.0|NW ı |=? 3.10||17 3.2| E 2 7| 3-29||14] 19.2| SE 4 0| 3.73 
As |18 23 NIV 222er) || TS 2.0) NE 2 6| 3.63||15)| 20.6|SSE 3 I} 3.30 
19 6.8| SW 3 251110 6.8|ENE 4 5| 3.39|\16| 14.2] NE 3 o| 3.46 
17 5.8| W 6 |=? 3.73\|17 3.9| ENE 15 US, 3-62||14 13:8 9922 26 8| 3.41 
V1 irs 4.8] SW 8 |= 3.43||18 3-7} ENE 11 |9° 7|3.36||15| 10.5] SW 5 R 9| 3.76 
19 38|SW 6 |=" |3.83||19 5.0| NE 8 o| —||16| 15.4|S cy |i" Al 
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